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1. INTRODUCTION 


In 1974 the Department of Geodetic Science at The Ohio State University (OSU) 
obtained observational data that was acquired during the West European Satellite 
Triangulation (WEST) program and the International Satellite Geodesy Experiment 
(ISAGEX) campaign. 

The purpose of obtaining this observational data was twofold. Primarily it 
was intended to perform a geometric solution to improve the present values of 
coordinates of the European stations in the OSU WN14 solutions [Mueller^ et al. , 
1973]. The secondary aim was to add some new stations and to assess the quality 
of the WN14 solution with the help of the additional data available. 

This report contains the details of the above investigation as follows: 

^'ection 2 the status of the data as received, the preprocessing required and the 
i.:. 3 limmary tests carried out for the initial screening of the data is discussed. 
Sections 3 and 4 describe the details of the adjustment computations carried out. 

The results of the adjustments are discussed in Section 5, which gives some 
concluding remarks on the entire investigation. 
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2. DATA 


2.1 WEST Data 

The imified optical prograirij otherwise known as WEST, was begun in 1966, 

The program was conducted by a subcommittee of the International Association of 
Geodesy (lAG). The program was formally terminated in 1972, as per resolution 
1 of the Sixth Meeting of the Subcommission in 1972 [lAG, 1972], Reduced data of 
approximately 3,500 simultaneous plates formed a major part of the data. How- 
ever, more data may be expected because the plate reduction is still in progress. 

WEST data was received by OSU in two forms. The first form consisted of 
cards in two sets. One set contained the direction cosines, referred to the Green- 
wich Hour An^e/Declination coordinate system, for a single fictitious image per 
plate, for all simultaneous events. The other set contained information about the 
standard errors associated with the sin^e image observational data derived from, 
the polynomial fitting. The second form contained the direction cosines of seven 
fictitious images per plate for all simultaneous events . The information about the 
standard errors of this observational data was not available separately. This report, 
therefore, contains only the results for single image data. 

The Ohio State University Geometric and Orbital (Adjustment) Program (OSUGOP) 
for Satellite Observations [Reilly, et al., 1972] cannot accept optical data in the Green- 
wich Hour Angle/Declination system. A modified version of OSUGOP which had been 
previously used [Reilly, et al. , 1972] to reduce the BC-4 seven image data, was 
further modified to accept single image observations. A separate program was made 
to merge the observational data of the single images on the first set of cards and the 
standard errors on the second set, thus forming an input set for the modified 
program. 

2.2 ISAGEX Data 

The Centre National d' Etudes Spatiales (CNES) provided OSU with data acquired 
during the ISAGEX satellite observation program. The data consisted of 5,186 laser 
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ranges and 3,562 optical observations which were purported to be simultaneous. 

The laser data was examined for simultaneity by computing the observation 
times at the satellite in Ihtemational Atomic Time for all data. With a criterion, 
of 0.2ms discrepancy for simultaneity, no simultaneous observations could be 
detected. Laser data, therefore, was excluded from the present investigation. 

Examination of the optical data was more complicated. The first test per- 
formed was to discover the amount of data which was simultaneous by using the 
above mentioned criterion. A wider definition of simultaneity was used to accom- 
modate any possible variations in observation times which mi^t arise from pre- 
processing. Approximately only ten percent of the data satisfied the test, indicating 
that further examination was warranted. 

The requirement for preprocessing was reconsidered. The ISAGEX Data 
Handling BooMet [Brachet, 1973] contained some of the needed information. Further 
information was obtained [ONES, 1973, 1972 and 1970] which indicated that pre- 
processing by OSU was possible, but not practical due to limited resources. Sub- 
sequently, Wolf Eesearch and Development Corporation offered to imdertake the 
preprocessing of the data using their GEODYN program. Observations made on 
satellites GEOS I and H, DIADEME-C and -D were preprocessed initially because 
their orbital elements were available at Wolf Research and Development Corporation 
for the observational period. Orbital elements for MIDAS and PAGEOS were la,ter 
provided by the Smithsonian Astrophysical Observatory (SAO). However, there were 
difficulties in obtaining the correct input data for preprocessing and it was decided 
not to use the observational data pertaining to either MIDAS or PAGEOS satellites. 

The preprocessed optical data was tested for simultaneity with a criterion of 
0.2 ms for discrepancy. A total of 353 observations proved to be simultaneous, 
involving 13 different stations. The acceptable ISAGEX preprocessed data could be 
input directly to the OSUGOP program for forming the normal equations and the 
subsequent adjustments. 
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3. WEST DATA PROCESSING 


3. 1 Transformation of Variances 

The variances of the observations were given in the form of standard deviations 
along ^d across the satellite trail. In a few cases no statistics were provided, and 
in order not to lose these valuable observations, the standard deviations as given in, 
[Ehmsperger, 1974, Table 4] were substituted. The modified subroutines which were 
to be used to compute the normal equations required as input, the standard deviation 
of the Greenwich Hour Angle multiplied with the cosine of the declination ( CTqha ' cos 6 ), 
the standard deviation of the declination (crs) and the covariance term. The observa- 
tional data also contained information about the length of the satellite trail and the 
declination for the beginning and the end of the satellite trail, which lead to the 
following transformations (Figures 1 and 2) to obtain Oqua ” cos 6 and a 6 for the 
observations. 



Transformation of Variances 
The following notation is used: 

U/ V/ : Topocentric Cartesian coordinate system parallel 

to the conventional geocentric U, V, W coordinate 
system . 

AL : Length of the satellite trail. 

6 i , 62 ' Declinations of the satellite at the beginning and 

end of the trail . 
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ot i Rotation an^e. 

Pi, Pg: Satellite trail. 

The actual rotation is approximated by a rotation around point Pi at the begin- 
ning of the satellite trail, taken as a strai^t line. 



Fig. 2 

Rotation to GHA System 


If X y^ represent the directions along and across the trail and x,y represent 
the required directions parallel to the plane, and along the tangent to the 
meridian of Pi , respectively, the following relation is obtained; 



( cos (- a) 

. ' 

- sm (“ a) 


sin(“ 0 !) 
cos (-or) 



where a. is computed from the spherical relation 

sin (62 “ 6i) 

sin Of = . 

sin AL 


Using the given variance -covariance matrix 



(3-1) 


(3-2) 


(3-3) 


where is the variance along the trail and is the variance across the trail, 
and by using the relation (3-1) and the law of propagation of errors, the transformed 
varianee-co variance matrix is obtained as given below; 
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(3-4) 



cos^( Od)a|+sin^(oi)ao |-sm(cx) cos{ a) (Cc- or^)\ 


~sin(a)cos(a)(ac- a®) isin® (a)a^+ cos®(a)ac 


It should be noted that cr stands for (Tqha • cos 6 and that (jy denotes <T6 . 

3.2 Observation Station Information 

3.21 Station Numbering System 

The whole data set contains observations from 30 different tracking stations. 

They are listed in Table 1 and their relative locations can be seen from Figure 3, 

Since the numbering system for the WEST stations and the WN14 stations are inde- 
pendent, there were some cases where different stations had the same station, numbers. 
In order to avoid confusion, the WEST stations were completely renumbered. The 
station number consists of four digits, the first two digits arbitrarily chosen as 87. 

In those cases where the tracking station was listed in [NASA, 1973], the NASA 
station number was maintained. 

3.22 Stations Common with WN14 Sys tern 

In order to compare and combine the WN14 and WEST systems, the common 
stations were identified, including those which could be established as common 
through relative constraints. These are marked in Table 1 with a superscript .2. 

3.23 Relative Station Constraints 

For the adjustment computations it is important to maintain the exact rela- 
tionship between nearby stations by introducing relative station constraints .In 
most cases the information for these constraints was extracted from the Circular 
Letters [WEST, 1966-1072]. In some cases the relative location of observation 
stations could be established from Cartesian coordinates given in [Ehmsperger, 

1974]. Details of these constraints are given in Table 2. 
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Table 1 

Tracking Stations 


Stafivn 


NAMI 


CQ«rdin«l*S la 

ED- so 


Ntight Ceattrsinft 

osu 

OtigiiMi 


• 

ivlilude 



o 

iangiltfdt 

/ 0 

EIL HI (M) 

Ell. HI. (M) 

a mi 



WEST 



STATIONS 






6006* 

14002 

TROMSO 

69 

39 

44.375 

18 

i 


119. 00 

113. 19 

4.0 

6016* 


CATANIA 

37 

26 

42.345 

15 

j 1 


-7.00 

16.33 

4.0 

6065* 


PEISEN 

47 

48 


11 

1 n 

28.574 

943.00 

960, 09 

2.5 

8004 

6004 

BRNSG 

52 

35 

05.286 


||l 

22.436 

77.17 

81.31 

3.0 

8009* 

9001 

DELFTH 

52 

00 

09.240 

4 



20.70 

41. 11 

4.0 

8010* 

12001 

ZIMLD 

46 

52 

40.318 

7 

27 

58.239 

900.34 

920.58 

2.5 

8011* 

13002 

MALVEN 

52 

08 

39.066 

358 

'nm 

59.567 

108.60 

134.97 

4.0 

8015* 


HAUTE PR 

43 

56 

01.140 

5 

42 

49.280 

651.00 

676.87 

4.0 

8016 

5002 

STRBG 

48 

35 

01.884 

7 

46 

11.135 

151.90 

165.93 

3.0 

8019* 

5004 

NICEFR 

43 

43 

36.496 

7 

18 


369.42 

394.73 

4.0 

8030^ 

5001 


48 

48 

25.354 

2 

13 

51.339 

155.46 

183.23 

2.5 

8031 

13001 


55 

44 


356 

46 

21.114 

285.10 

301.37 

3,0 

8032 

6110 

HOPBG 

47 

48 

08.287 

11 

01 

26.245 

939.30 




6010 

HOPBG 










8033 

6005 

FRNFT 

Bril 

13 

14.257 

8 

43 

51.822 

177,70 

188.07 

3.0 

8034 

9002 

DELFY 

52 

02 

43,850 

4 

21 





8701 

1001 

GRAZA 

47 

04 

03.821 

15 

29 

40.117 

484. 60* 

492.39 

3.0 

8702 

2001 

BRXOR 

Bill 

47 

53.600 

4 

21 

37.750 

100.70 

115.29 

3,0 

8703 

3001 

COPHN 

55 

44 

22.064 

12 

30 

04.101 

52.89 

52.39 

3.0 

8705 

5003 

BRDUX 

44 

50 

06.500 

359 

28 

24.600 


108.82 

3.0 

8706 

5005 

GOULT . 

43 

51 

12.069 

5 

13 

34.032 

200.30 



8710 

6012 

WSNDF 

52 

35 

05.328 


lail 

22.523 

76.41 



8711 

8004 

GATAN 

37 

26 

42.717 

15 

02 

47.326 




8712 

8005 

OPICI 

45 

40 

59.268 

13 

46 


383.28 

395. 12 

3.0 

8713 

8006 

ORIAA 

40 

30 

01.079 

17 

38 

32.862 

179.90 

196.50 

3.0 

8714 

8007 

SRDIN 

39 

13 

20.960 

9 

07 

04.500 

116.97 

146. 53 

3.0 

•8715 

8008 

' TANIA 


41 

39.050 

14 

58 

31.605 

1718. 12 



8716 

10002, 

MADRD 

40 

27 

07.699 

356 

16 

31.979 

642.38 

688.56 

3.0 

8717 

10003 • 

MADRI 

40 


07.687 

356 

16 

31.865 

642. 38 



8718 

6006 

KLSRH 

49 

Bil 

43.050 

8 

24 

43,760 

132.50® 


3.0 

8719 

8009 

CATNA 

37 

41 

39.050 

14 

58 

31.605 

1718.12 



8720 

11002 

LOVOA 

59 

20 

17.576 

17 

49 

48.173 

49.40 

41.48 

3.0 

8722 

15001 

REKVK 

63 

57 

44.980 

337 

24 

40.020 

-26.17® 

0.40 

3.0 



1 S AG E X 



stati ons 




■i 

8609* 

8009 

ST.mCH 

43 

55 

59.186 

5 

42 

48.383 

650.20 


■ 

872l“ 

1147 

ONDREJOV 

49 

55 

19.4 

14 

48 

03.90 

508.61 

512.98 


8723® 

1181 

POTSDAM 

52 

22 

51.4 

13 

03 

58.80 

109.00 

111. 81 


9004* 

9004 

SAN FERN 

36 

27 

51.367 

353 

47 

42.091 

-12.00 

50.44 



^Based on [Ehmsperger, 1974] except where 
specified by superscript 7 and 8. 

^Stations used for tying to the WN14 system. 
^Recomputed in European Datum (ED). 
^Astro Datum, 


®ED values interpolated. 

^Determined by WEST and reintroduced 
as starting coordinates. 

^Based on [NASA, 1973]. 

*For details on soiu-ce see Section 4. 
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Tracking Station Locations 

A WEST Station 
• ISAGEX Station 
A West and ISAGEX Station 
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3.3 Adjustment Computations 


The observational data was given in the form of direction cosines in the Average 
Terrestrial Coordinate System, After conversion to the Greenwich Hour An^e/ 
Declination system, the observations were used in conjunction with the transformed 
variances to form a set of normal equations by the modified OSUGOP subroutines 
as given in Appendix IV., A total of 396 events proved to be satisfactory. An event 
is defined as two or more groxmd stations simultaneously observing one satellite 
position. Table 4 gives the number of such observations on each line. All obser- 
vations turned out to be of good quality as verified by means of the ’’test distance" 
while computing the normals [EeUly, et al. , 1972] . 

3.31 Preliminary Adjustment 

Althou^ the variances are transformed to the Greenwich Hour Angle/ Decli- 
nation system, they still represent only the errors from the smoothing procedure 
[Ehmsperger, 1974], This procedure yields an overly optimistic estimate for the 
variances . Thus preliminary adjustments were carried out using only minimum, 
constraints to find Ihe proper scaling factor for the variances. The variance of 
xmit weight was found to be 36 . 

3. 32 Solution WEST 34 

This adjustment contains only observational data acquired during the WEST 
campaign. The version WEST 33 was submitted to the XVI General Assembly of the 
lUGG held in Grenoble, France on August 18 through September 6, 1975. While 
working on the combination solution for the WEST and ISAGEX data it was decided 
to add to the previous WEST 33 solution, the additional relative constraints between 
station 8015 (HAUTE PKOVENCE) and station 8706 (GOULT). Also, some station 
numbers were modified, giving priority to the NASA station numbering system. 

But above all, it is recognized that the baseline between 6006 (TROMSO) and 6016 
(CATANIA) is not sufficient to transfer scale to the whole network. The WEST 
satellite network is considered as consisting of two blocks: the central European 

block with a large number of observations and the northern block, connected to the 
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central block by relatively few observations— namely, between 6006 (TROMSO) 
and some stations of the central block. An overall scale factor of 10 ppm. was 
computed between the ED50 coordinates and the adjusted values. Comparing 
individual chords in the two systems it became clear that all chords originating 
from 6006 (TROMSO) yield a significantly smaller scale factor. Also, the scale 
factor in the central area is partly inherent in the weighted positional constraints 
of the WN14 stations. It thus became necessary to include more chord constraints, 
especially in the central area. These were taken from [Ehrnsperger, 1974]. 

The present adjustment, WEST 34, contains the following types of constraints 
appropriately weighted: 

a) A priori constraints on coordinates of observation stations which are 
common with the WN14 solution (Table 1) . These coordinates were 
constrained at the values obtained in the WN14 solution and weighted 
as per their variances. 

b) Relative position constraints (Table 2). 

c) Height constraints for stations common to the WN14 solution as per 
[Mueller, etal., 1973, Table 3.3-3]. Ellipsoidal height constraints 
were also applied for the other stations after transforming the Euro- 
pean Datum height information available [Ehrnsperger, 1974] to the 
WN14-system. This was accomplished by means of iterations. For 
the first direct transformation from the ED50 to the WN14 system, 
the mean values of the shifts given in [Anderle, 1974] were taken and 
corrected for the offset of the origin of the WN14 system from the 
geocenter [Mueller, et al. , 1973]. The height constraints for the 
final iteration are shown in Table 1. 

d) The scale was introduced through three base lines (Table 3) . 

The values of the station coordinates as a result of this adjustment are given in 
Table 5. 
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Table 2 

Relative Station Constraints (WEST) 


1 station 

■BH 

■aH 

WMjm 

Source 

from 

to 





8711 

8719 

13329.91 

10072.67 

-22958.92 

3 

8711 

8715 

-4.04 

-8.24 

7.88 

2,3 



1.78 

-1.33 

-0.16 

2,3 

8716 

8717 

-0.06 

2.70 

0.28 

2,3 


■aiMM 

3709.06 

1053.54 

-2925.82 

3,4 

8711 


-1.61 

-0.43 

2.17 

2 

6065 

8032 

20.79 

53.37 

-25.09 

4 

8015 

8706 

-9569.56 

38443.07 

6742.48 

3 , 4 : 



ISAGEX 



8609 

8015 

43.06 

-15.80 

-43.99 

4 


Table 3 

Base Lines (WEST) 


Station > 
from to 

Chord Distance (m) 

m 

Source ** 

6006 

6016 

3 545 871.454 

3.5 

.1 

8011 

8032 

1 046 479 . 89 

0.65 

3 

8032 

8701 

346 945 . 38 

0.27 

3 

6006 

8032 

2 457 733.52* 

1.20 

3 

8032 

8711 

1 194 842.47* 

1.20 

3 


* Base line not used in this adjustment. 

** Sources used in these tables: 

^[Mueller, etal., 1973] 

^[WEST, 1966-72] 

^[Ehrnsperger, 1974] 

“^[NASA, 1973] 
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FRNFT 

DELFY 

CRAZA 

BRXOR 

COPHN 

BRDUX 

GOtJLT 

\VSNDF 

CATAN 

OPICI 

ORIAA 

SRDIN 

TAMA 

MAORD 

MAORI 

KLSRH 

CATNA 

LOVOA 

REKVK 


8033 

8034 

8701 

8702 

8703 

8705 

8706 

8710 

8711 

8712 

8713 

8714 

8715 

8716 

8717 

8718 

8719 
872.0 
8722 
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Table 5 


Cartesian and Geodetic Coordinates 
(Solution WEST 34) 


Sta.No 


Uu 

V 

or. 

w 

Uv 


<0 


\ 


H 

cr« 


a. A. r. 

Sc- . Ae Te 


u,v,w" Cartesian coordinates in meters (Orientation; u = the Greenwich 
meridian as defined by the B.I. H. ? v ^ X - 90" (E); w ~ Conven- 
tional International Origin). 

tp,X Geodetic latitude and longitude in angular units (degrees, minutes 
and seconds of arc) computed from the Cartesian coordinates and 
referred to a rotational ellipsoid of a = 6378155. 00 m and 
b = 6356769.70 m. 

H Geodetic (ellipsoidal) height in meters referred to the same 

ellipsoid. 

Standard deviations of the Cartesian coordinates in meters. 
Standard deviations of the geodetic coordinates in seconds of arc. 
Oh Standard deviations of the geodetic height in mfjters. 

r* Altitude (elevation angle), azimufh. and magnitude of the major 
semi axis of the error ellipsoid, respectively. Angles in degrees, 
, magnitude in meters. Altitude is positive above the horizon. 
Azimuth is positive east reckoned from the north (see section 
4.74). 

at,A)„rb Same as above for the mean axis of the error ellipsoid, 
act Ac, r« Same as above for the minor axis of the error ellipsoid. 
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Table 5 (Continued) 


EEPRODUCIBILITY OF 
OBIGINAL PAGE IS POOH 


6006 

2102928.96 

mm 

HRSBBRI 

2.67 

5958181.83 

2.41 


69 39 45.16 

0.07 

IP 56 26.83 

0.25 

112.55 

2.41 




C .18 

99.14 

2.66 





82.94 

7.70 

2.42 






7,06 

- 170,8 3 

2.22 



6016 

A 8963 P 8.49 

1.74 

13161 70.99 

2.27 

3856668,07 

2.16 


37 26 39.09 

0.07 

15 2 44.56 

0.09 

15.57 

1 . 90 



- 1,68 

109.42 

2.28 





37.63 

20.73 

2.15 





52.12 

“ 162.75 

1.74 



6065 

4213563.93 

1.89 

020824.78 

2.18 

4702785.56 

2.09 


47 48 4,56 

0.07 

11 1 24.47 

0.11 

959.34 

i .75 



15.09 

19.52 

2.24 





- 1.55 

109.10 

2.18 





74.83 

- 166.63 

1.71 



8004 

3818505.85 

4.97 

708050.10 

5.13 

5042639,98 

4.74 


5 ? 35 3.00 

0.20 

10 30 17.33 

0.28 

82.64 

2,99 



- 6.64 

- 156.28 

6.29 





- 0.62 

113.65 

5.04 




• 

83.33 

“ 161 .70 

2 . C -92 



8009 

3923399.47 

4.21 

2998 B 2 .98 

5.19 

5002971.24 

4,05 


52 0 6.35 

0.16 

4 22 15.15 

0.27 

42.67 

3.16 




0.83 

127.83 

5.41 






- 6.47 

-142 .27 

4.66 






8 3.47 

- 134,93 

3.13 



8010 

4331299.33 

3.47 

567503.55 

4.92 

4633113.20 

3.60 


46 52 37.22 

0.14 

; 7 27 52.53 

0.23 

920.42 

2.31 




- 0.95 

113 .37 

4.99 






3.15 

• 23.43 

4,36 






86,71 

-173 .34 

2.30 
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Table 5 (Continued) 


REPRODUCIBILITY OF 
ORIGINAL PAGE IB POOR 


8011 

3920172.38 

3.05 

- 1347A . 3.98 

2.64 

5012722 . 2 a 

3.A7 


52 8 35.59 

0,12 

356 1 53.07 

0.14 

179,31 

2,86 



16.58 

25.61 

3.91 





71.18 

176 . 5 ? 

2.75 





*“8 o 6 ^ 

113.01 

2.25 



8015 

A578317 , 7 «> 

3.36 

4579 A2 .24 

4 .26 

4403192.90 

3.64 


43 55 57.83 

0.13 

5 42 43.95 

0.19 

676.18 

2.74 



2.80 

124.98 

4.33 





7.60 

34.61 

4.08 





81.89 

“ 124.97 

2.71 



8016 

4188648.12 

3.23 

571416,21 

4.09 

4760143.73 

3.77 


48 34 58.73 

0.14 

7 46 6.15 

0.20 

165 . 3 ? 

2 .42 



11.42 

“ 5.69 

4.39 





3.90 

85.10 

4,09 





77.92 

“ 166.35 

2.29 



8019 

4579462.47 

3.64 

586587,06 

6.90 

4386418.66 

3.83 


43 43 33.26 

0.14 

7 17 57.54 

0.31 

394.78 

3.11 



2.98 

98.63 

6.92 





3.39 

8 .4 5 

4,27 





85.48 

“ 130.20 

3.09 



8030 

4205629.67 

4 . 6E 

163702.54 

7.94 

4776538.93 

4.14 


48 48 22.12 

0.19 

2 13 44.73 

0.39 

183 . P9 

2.47 



“ 0.73 

86.18 

7.96 





0.96 

176.16 

5.71 





88.79 

33.46 

2.46 



8031 

3593858.98 

9,44 

“ 202756,98 

6.44 

524SC6C.45 

6 , 9.8 


55 44 0.92 

0.36 

356 46 15.35 

0,38 

290.05 

3.05 



“ 1.03 

161.27 

11.68 





4.16 

71 .34 

5 .82 





85.71 

“ 122.59 

3.02 



8032 

4213543.14 

1 . 9C 

820771,38 

2.19 

4702810.63 

2.10 


47 48 5.84 

0.07 

11 I 22.14 

0.11 

957.33 

1.76 



15.18 

19.1 0 

2.2 5 





“1 .47 

108 ,7<-> 

2.19 





74,74 

“ 166,71 

1.72 
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Table 5 (Continued) ^^EPjRODUClBILIT^’ 

ORIGINAL page IS P * = • 


eo?3 

4041R64.41 

4.P9 

f?06^:»0,3n 

6.00 

487862'’. 41 

4 ,56 


50 13 11.30 

0.20 

8 43 46.61 

0.30 

189.71 

2.00 



1.08 

142.60 

6.39 





4.06 

52.52 

5.61 





85.80 

“112.47 

2.88 



8034 

3910690,42 

4.21 

20 88 29 .43 

5.19 

5005897,06 

4,05 


52 2 40.96 

0.16 

4 21 34.86 

0.27 

23.92 

3,16 



C.80 

127.83 

5.42 





“6.51 

“142.26 

4,60 





83.44 

-135.2?. 

3.14 



8701 

4194425.71 

3.87 

1162696.85 

2.22 

4647199,64 

4.20 


47 4 0.67 

0.16 

15 20 36.24 

0.11 

487.20 

2.89 



8.19 

12.71 

4.07 





81.68 

“157.11 

2,83 





1.45 

102 .02 

2,20 



8702 

4027918.87 

4.78 

307001.98 

7.69 

4910436.33 

5.33 


50 47 50.83 

0.21 

4 21 30.86 

0.39 

115.71 

3,01 



“3.18 

123 .88 

8.26 





10.10 

34.45 

5.83 





79,40 

-163.37 

2.84 



8703 

3513633.88 

7.03 

776936.23 

6.25 

5248194.03 

4.95 


55 44 20.31 

0.25 

12 29 59.01 

0.39 

53.06 

3.08 



0.03 

“144.8 0 

8.39 





“l.Fl 

125.23 

5.77 





87.97 

on .05 

3.07 



8705 

4530506.30 

17.35 

“41732.26 

33.17 

4474382.07 

18 ,05 


44 50 2.77 

0.81 

359 28 20.07 

1.50 

108.90 

3.21 



0.07 

54.73 

39.47 





-2.12 

144.73 

12.66 

. 




87.87 

146.65 

3.18, 



8706 

4587887.34 

3.36 

419519.19 

4,25 

4306450,50 

3.64 


43 51 8.73 

0.13 

5 13 28 .67 

0.19 

226.34 

2.74 



2.51 

125 .20 

4.33 





7.40 

34.8 7 

4.08 





82.18 

“126 .21 

2.71 
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Table 5 (Continued) 


8710 

3818504,07 

4.96 

7080 51 .4 3 

5.12 

5042640,13 

4.74 


52 35 3.04 

C .20 

10 30 17,42 

0.28 

81.35 

2.09 



“ 6.69 

- 1 56 ,2 8 

6.29 





- 0,65 

113 .64 

5.03 





83. 28 

“ U 1 .89 

2.92 



8711 

4896386.88 

1.75 

1316170.55 

2.27 

3856670,24 

>.17 


37 26 39,18 

0.07 

15 ? 44.86 

O.OQ 

15.56 

l.oi 



“ 1.80 

109.61 

2.29 





37.85 

21.02 

2.16 





52 . C 9 

“ 162.70 

1.75 



8712 

4335518.63 

17.61 

1063083 .44 

30.07 

4540932.63 

10 . 4 ? 


45 40 5 5.82 

0.89 

13 46 38.24 

1.32 

395,00 

3.21 



“ 0.08 

132.2 7 

33.42 





0.00 

42.27 

21 .56 





89.92 

133.65 

3.21 



8713 

4628609.90 

19.49 

1471955.74 

39.16 

4120468.50 

27.30 


40 29 58,72 

1 » 16 

17 36 28.73 

1.57 

196.53 

3.22 



“ 0,10 

133.07 

44.42 





“ 0.21 

43.07 

25,98 





89.77 

69.10 

3.21 



87 lA 

4885403.52 

23.58 

784066.50 

29.46 

4011526.62 

3>.33 


39 13 18.48 

1.35 

9 7 3.87 

1.13 

146.37 

3.21 



0.11 

“ 27.01 

45,54 





0.65 

62.99 

19.62 





89.34 

“ 126.53 

3.20 



8715 

4896390.03 

1.76 

1316178 .78 

2.28 

3856662 .36 

2.18 


37 26 38.86 

0,07 

15 2 44.84 

0.09 

15.57 

1.92 



“ 1.88 

109,74 

2.29 


■ 



37.86 

21.20 

2.17 





52.08 

- 162.68 

1 .75 



8716 

4850674,55 

8.07 

- 3150 C 7 , 3 ? 

7.52 

4116626.07 

9.97 


40 27 1,36 

0.41 

356 16 25.64 

0.31 

689,70 

3.15 



0.97 

21 .94 

13.41 





“ 3,20 

111.89 

5.61 





86.65 

128 . 7 P 

3.13 
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Table 5 (Continued) 


8717 

485067A^,f.l 

F.Op 

-31 5910,03 

7.52 

4116626.60 

9.97 


40 27 1.35 

0.41 

356 16 25.53 

0.3J 

689,70 

5.15 



0.97 

21 .94 

13.41 





“3.?l 

111.8B 

5.61 





86.65 

128 .7 5 

3,13 



6718 

4146533.64 

6.12 

613109 .P3 

6.50 

4791487,60 

5,56 


49 0 39,75 

0.25 

e 24 39.09 

0.32 

146.23 

3.13 



C.37 

175.37 

7.71 





3.59 

85.35 

6.45 





86.40 

“88. PI 

3.11 



8719 

4883056. 9R 

1.76 

13060'37.PP 

2.28 

3679620.16 

2.18 


37 41 35,53 

0.07 

14 58 28.82 

0,09 

1741.09 

1.92 



”-l.oR 

109.6 3 

2.29 





38.07 

21.18 

2.17 





51.86 

“162. 9C 

1 ,75 



8720 

3104204.46 

10.54 

9983 50 ,28 

8,55 

5463280.09 

5,62 


59 20 16.53 

0,32 

17 49 42.89 

0.66 

41,o«? 

3.13 



0.91 

-133.00 

12.83 





“0.40 

137 .0 1 

6.45 





89,00 

70.71 

3.12 



8722 

2591994.28 

10.16 

“1078495.16 

18.33 

5707863.39 

5.85 


63 57 42.52 

0.36 

337 24 30.38 

1.36 

0,47 

3.20 



O.OB 

78,45 

18,75 





“1.63 

168.45 

10.56 





88.17 

170.97 

3.18 




NORMAL TERMINATION 
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4. COMBINED WEST-ISAGEX SOLUTION NO. 36 


This adjustment is a combination solution using the optical observations from 
both WEST and ISAGEX. The ISAGEX normal equations were added to the previous 
set of the WEST 34 solution after appropriate scaling as described in [Reilly, et al. , 
1972 and MueUer, et al. , 1973], The scale factor was derived from a separate solu- 
tion with minimum constraints . Since the ISAGEX data set contains no information on 
the statistics of the observations, the latter were treated with equal weights. 

Only the preprocessed data as provided by Wolf Research and Development Cor- 
poration was used. During the formation of the normal equations, several blunders 
in observational data were indicated . A closer Inspection of the data shows that all 
observations which satisfy the simultaneity criterion are GEOS II flashing light ob- 
servations. The observation times show the characteristic four second interval from 
the zero to the twenty -fourth second of a minute. It is obvious that the blunders were 
a result of mis identifying the GEOS II flashes or timing errors , Attempts were made 
to correct this error by simply rearranging the observations. The results, however, 
are unsatisfactory, probably because such misidentification also affected the prepro- 
cessing. Thus, a mere rearranging of data does not rectify the situation. 

Table 6 illustrates the additional simultaneous ISAGEX observations used in 
this adjustment data. Originally, several observations that satis fi ed the 
test distance criteria, were extracted for the stations 8019 (NICEFR) and 8031 
(EDNBG) . Their use in the adjustment caused an extraordinarily large shift of approx- 
imately 20 to 40m in each coordinate for the corresponding stations, even though the 
station identification seems to be correct, especially for station 8031 (EDNBG). 

But since no reasonable explanation could be found for this significant change, these 
observations have been deleted. 

Station identification could be uniquely resolved for the stations 8010 (ZIMLD), 
8031 (EDNBG) and 8034 (DELFY) through [Schuerer, 1972; Mclnnes, 1972 and Aardoom, 
1972]. Since there is no special ISAGEX numbering system (see [Brachet, 1972]), it 
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is assumed that the ISAGEX station 8011 (MALVRN), 8019 (NICEFR) and 9004 
(SAN FERN) are identical to the NASA stations of the same numbers . 

The coordinates as given in [Marsh, et al . , 1975] were taken as approximate 
coordinates for stations 8721 (ONDREJOV) and 8723 (POTSDAM), The height con- 
straints could be derived, after appropriate transformation, from MSL heights 
given in [Karsky, et al. , 1974 and Bracket, 1973]. The undulation of Potsdam was 
taken as zero; while that for Ondrejov was interpolated. (See also Tables 1 
and 2. ) The values of the station coordinates as a result of this adjustment are 
given in Table 7. 


Table 6 


Distribution of ISAGEX Observations per Line 




8011 

8034 

8609 

8721 

8723 

9004 

ZIMUD 

8010 


13 

16 

11 


13 

MALVRN 

8011 


6 



5 

7 

DELFY 

8034 



6 

4 


24 

ST MICHEL 

8609 




8 


12 

POTSDAM 

8721 






5 

ONDREJOV 

8723 






9 
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Table 7 


Cartesian and Geodetic Coordinates 
(Solution WEST - ISAGEX 36) 

REPRODUCIBILITY OP’ T -.: 
ORIGINAL PAGE IS POOi^ 



u,v,w Cartesian coordinates in meters (Orientation: u s the Greenwich 
meridian as defin^bjr the v s X “ 90“ (E); w s Conven- 

tional International Origin). 

Geodetic latitude and longitude in angular units (degrees^ minutes 
and seconds of arc) computed from the Cartesian coordinates and 
referred to a rotational ellipsoid of a = 6378155. 00 m and 
b = 6356763.70 m. 

H Geodetic (ellipsoidal) height in meters referred to the same 

ellipsoid. 

OutOrtOv Standard deviations of the Cartesian coordinates in meters. 

Standard deviations of the geodetic coordinates in seconds of arc, 

<7m Standard deviations of the geodetic height in meters. 

a»»A„r, Altitude (elevation angle), azimuth and magnitude of the major 

semi axis of the error ellipsoid, respectively. Angles in degrees, 
magnitude in meters. Altitude is positive above the horizon. 
Azimuth is positive east reckoned from the north (see section 
4.74). 

a^,A^, Tb Same as above for the mean axis of the error ellipsoid. . 
ac»Ae,rc Same as above for the minor axis of the error ellipsoid. 
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Table 7 (Continued) 


SS£”Sf®io» 


f>00f> 

2107929. 

05 

2.07 

721666.39 

2.61 

5958182.41 

1 .00 



Bi 

0.07 


0,24 

113.13 

0,69 


nnm 


0.01 

99 . 44 . 

2.63 





- 

G.13 

-170,56 

2.17 





89.87 

-167.77 

0,69 



6016 

AB963G8. 

85 

l.?o 

1316171.20 

2.17 

3856668.64 

1.62 


37 26 39. 

10 

0.06 

15 7 94.56 

0.09 

16,23 • 

0,67 



• 

0.03 

107 .57 

2.22 






1.73 

17,57 

1.91 





8G.27 

-163.39 

0.67 



6065 

4213564. 

02 

1.62 

820824.14 

2.10 

4702786.55 

1.47 


47 48 4., 

58 

0.07 

11 1 24.44 

o.in 

960.04 

0,43 




0.37 

-137.34 

2.14 





“0.21 

132,65 

2.11 





89.58 

-166.74 

0.43 



800^ 

3818505. 

?o 

4.89 

70 8049 .15 

5.07 

5042640.20 

3.66 


52 35 3. 

03 

0.19 

10 30 17.29 

0.28 

82.32 

0,53 



- 

0.17 

-155.46 

6.14 


, 



- 

C.02 

114.54 

5.01 





89.8? 

-160.58 

0.53 



8009 

3923397. 

97 

3. IP 

299876,48 

3,68 

5002970.72 

2.50 


52 0 6. 

39 

0.13 

4 22 14.81 

0,19 

41.03 

0,69 




6*21 

-174,44 

3.98 






0.08 

95,56 

3.69 





89.77 

-15.70 

0.69 



8010 

433129E. 

35 

2,36 

567502,32 

3.58 

4633114.47 

2.16 


46 5? 37. 

27 

0.10 

7 27 52.48 

0.17 

920.57 

0.44 




0.03 

81 .00 

3.63 






0,07 

171 .00 

3.12 





89.9? 

-32.26 

0,44 
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Table 7 (Continued) 


8011 

30 ? o 160.55 

2.69 

- 1347 ^ 5.00 

2.52 

5012719.45 

2.19 


52 . e 35,60 

G.ll 

*3 5 f ‘ 1 53.01 

0.13 

135.30 

0,60 



- 0.48 

- 155.77 

3.62 





0.52 

114.23 

2.18 





69,20 

- tl 08.41 

0.69 



8015 

657 S 322.26 

2.10 

457962.32 

3.45 

4403189,27 

2.18 


43 55 57.65 

0.10 

5 42 *, 3.94 

0.15 

676.81 

C .69 



0.52 

94.79 

3.45 





0 , 32 . 

4.79 

2-94 





80,30 

“ 116.8 9 

0.69 



8016 

4198647.81 

3.07 

571417.65 

3,96 

4760 144 , PI 

2.76 


48 34 58.76 

0.13 

7 46 6 cl 2 

0.19 

165.68 

0,52 



“ 0.43 

160 .; S ‘0 

4.12 





0.28 

70,20 

3.94 





89.40 

“ 166.16 

0.52 



8019 

4570462.46 

2.94 

586586 .05 

6.73 

4386418.74 

3,10 


43 43 33,27 

0, 14 

7 17 57.49 

0.30 

394.73 

0.70 



0.13 

98.57 

6.77 





0,09 

8.57 

4,16 





80,85 

“115 ,60 

0.70 



8030 

4205628.97 

4.24 

163700.74 

7.80 

4776538.76 

3.70 


48 48 22.13 

0.18 

2 13 44,64 

0.38 

183.25 

0,44 



“ 0.02 

86 .56 

7.82 


• 



0.03 

176.56 

5.58 





89.97 

34.03 

0.44 



8031 

3593850.57 

9.46 

“ 202757.45 

6 .47 

5248061.70 

6.35 


55 44 0,93 

C .36 

356 46 15,33 

0.38 

301 . 3 ? 

0.53 



- 0.03 

160.96 

11.75 





0.09 

70,9 6 

5,79 





89.90 

“ 126.40 

0,53 



8032 

4213543.23 

1.62 

820770.77 

2.10 

4702811 .61 

1.47 


47 48 5.86 

0,07 

11 1 22 . 1 1 

0.10 

958.04 

0.43 



“ 0.37 

“ 137.43 

2.14 





“ 0.21 

132.57 

2.11 





89 .58 

“ 165.83 

0.43 
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Table 7 (Continued) 


P033 


80?4 


860‘» 


8701 


8702 


8703 


8705 




pOOB 


4041862.78 

4,62 

620629,85 

5.88 

487862f ,75 

>.93 

50 13 11.42 

<•.20 

F 43 46.60 

0,29 

IF.8.12 

C,53 


0.03 

145,85 

6.35 




0.11 

55.85 

5.54 




89, £9 

-110.84 

0.53 



3919688.91 

3,18 

298822.94 

3.68 

5005R9<^,5/; 

2,50 

52 ? 41.00 

0,13 

4 21 34.53 

0,19 

22.27 

0,69 


0.17 

-174.45 

3.98 




C.07 

95.5 5 

3.69 




89.81 

-17,21 

0,69 

- 


4578365.32 

2.10 

4579 46.52 

3.45 

4403145,28 

2.18 

43 55 55.69 

0.10 

5 4? 43,04 

0.15 

676.02 

0,69 


0,52 

94.79 

3.45 

. 



C.3? 

4.79 

2.94 




89.39 

-116.88 

0.69 


■* 


i 

i 

I 


4194427.35 3.64 

47 4 0.75 0.15 

-0.19 

-C.C4 

E-9.8C 


402791P.19 4.74 

50 47 50.85 0.21 

“O.OP 

-0.27 

69.72 


1152697.02 2.09 

152936,23 0.11 


-157.13 4.86 

102.87 2.12 

-178.89 0.53 

307000.05 7,58 

4 21 30.76 0,39 

120.93 8.17 

-149.07 5.54 

-165.00 0.53 


4547204. P7 3.28 

492.22 0.53 


4919436.45 4,09 

115.30 0,53 


3513633.03 6.B2 

55 44 20.34 0.25 

0,02 

-0.04 

69.95 


778935.33 6.26 

12 29 58.97 0.39 

-144.00 8,42 

125. 1C 5.76 

92.60 0.53 


5248193.93 4.?2 

52.41 0,53 


4530507.04 17,56 

44 50 2.73 0.82 

0.00 

-0.05 

69,95 


-4 17 34.52 33,58 


359 28 19.97 1.52 

54.78 39,95 
144.78 12.67 
146.77 0.5? 


4474381.19 17.92 

108.87 0,53 
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Table 7 (Continued) 


B 706 

A 587 fl ‘> l . P 2 

2.10 

419519.25 

3.45 

43964 ^ 6.79 

2.18 


A 3 51 8.55 

0.10 

5 13 28 .66 

0.15 

227.00 

0.69 



0.18 

94.46 

3.45 





0.21 

4.46 

2.94 





89.72 

- 135.38 

0.69 



8710 

3818503 . A 2 

4,80 

708050,48 

5.07 

5042640.36 

3.66 


52 35 3.07 

C.io 

10 30 17.37 

0.20 

81.53 

0.53 



- 0.17 

- 1,55 .45 

6.14 





- 0,02 

114.55 

5.01 





89.83 

- 160,67 

0.53 



8711 

4806387.24 

1.30 

1516 X 70.77 

2.17 

3856670.81 

1.62 


37 2 A 30.10 

0.06 

15 2 44.56 

0,09 

16,23 

0.67 



- 0.03 

107.57 

2,22 

* 




1.73 

17. 5 f 

1.91 





88.27 

- 163.39 

0.67 



8712 

433551 8.00 

17.72 

1063083.85 

30.46 

4540933.17 

19 . 5 i 


45 40 55.55 

0.90 

13 46 38.26 

1.34 

395.11 

0.53 



- 0.00 

132.27 

33.8 5 





0.00 

42 .27 

21.83 





90,00 

153 .73 

0.53 



8713 

4628610.14 

19.57 

1471055.77 

39.67 

4120 ^ 68.17 

27.57 


40 29 58.70 

1.15 

17 38 28.73 

1 .50 

196,50 

0.53 



- 0.00 

133.11 

44 . cc 





- 0.01 

43 .1 1 

26.30 





89,99 

69.32 

0,53 



871 A 

4885402.98 

23.85 

78 40 67.24 

29.82 

-^ 011527.38 

32.58 


39 13 18.51 

1.36 

9 7 3.90 

1.14 

146.52 

0.53 



0.00 

- 27.06 

46,12 





0,02 

62.94 

19.81 





89.98 

- 126.70 

0.53 



8715 

4806391.28 

1.30 

1316179 .00 

2.17 

3856662.03 

1.62 


37 26 38.66 

0.06 

15 2 44.34 

0.09 

16.23 

0,67 



- 0,03 

107,5 8 

2,22 





1.73 

17.58 

1.91 





88.27 

- 163 .39 

0.67 
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PT^PRODUCIBIUTY of the 
T able 7 (Continued) OPIGINAIj PAGE IS POOR 


6716 

A8 5(-67?.5r) 

7 .95 

“315«^( 8 .11 

7. 57 

9116626 ,23 

0.65 


<•0 2 7 1.36 

• •.91 

356 16 25.61 

0.31 

668,59 

0,53 



0.03 

22.12 

13.98 





~C.06 

112.12 

5,61 





89,9 3 

139.0 5 

0,53 



6717 

067 3. 61 

7.95 

^315910.81 

7.57 

9116625,95 

0.65 


90 27 1.25 

0,91 

256 16 ?5.9«5 

0.31 

6E8.59 

0.53 



C.03 

22ol2 

13.9 8 





~0,C6 

112.12 

5 ,6 jj. 





80,93 

124.02 

0,53 



8718 

9196522.13 

5.83 

613106.77 

6.92 

9791987,56 

3,08 


99 ff 39.79 

• '.25 

6 2^t 39«i>5 

0.32 

195.12 

0,53 



0.01 

176.8 5 

7.72 





o.oe 

t-6 .6. 5 

6,9i 





89.92 

“f6.89 

0,53 



8719 

9882057.33 

1.30 

1306098. 10 

2.17 

3879629.73 

1,62 


57 9x 35.59 

C.06 

19 58 28.82 

0.09 

1791 .79 

0.67 



-0.16 

107.53 

2,22 





1 ,95 

17.59 

1.91 





8 8.09 

“167.09 

0,67 



8 720 

3109203,00 

10.91 

998359 ,38 

8.62 • 

5963280.35 

5.11 


59 20 16,58 

C ,3? 

17 99 92.92 

C . 66 

91, *9 

0,53 



<:.02 

-132. <^5 

12.91 





-O.ljl 

127.05 

6.98 





£9,97 

66.16 

0.53 



8721 

3978930.80 

6.59 

1051053.13 

6.80 

9857553.91 

3.66 


99 55 16,19 

0.19 

19 97 59.99 

0.39 

519,73 

3.02 



-5. 89 

57.19 

8.66 





8.90 

1*6.33 

9.17 





79,75 

1 .69 

2.86 



8722 

2591995.20 

10.29 

•-107 89 99.51 

18.59 

5707662. OR 

9.80 


63 57 92,50 

C.36 

337 29 30.95 

1.37 

0.90 

0,53 



0.00 

78.3? 

18 .96 





-r.05 

1 68 . 3 3 

10.61 





89,95 

17i.no 

0.53 
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ORIGINAL PAGE IS POOR- 


872.3 

3800631.86 
52 22 4P.66 

6.?0 

0.14 

f810A4,lo 
1? 3 51.76 

12.31 

0.69 

5028??4.8o 

112.55 

4.31 

3,10 



-1.6? 

-2.13 

87.24 

77.12 

167.1? 

128,84 

13.40 

5.23 

3.07 




9004 5i05579.7o ?.J-4 “?‘5524f'.5 1 4,7? 

34 27 44,45 0.12 353 47 35,42 0.14 


1.60 58.22 
12.94 -32.15 
76.46 155.13 


5,1? 

3.10 

2 . 6 ? 


3760675,36 3.*^3 

47,65 2.66 


NORMAL TERMINATION 
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5. DISCUSSION 


Both solutions WEST 34 and WEST-ISAGEX 36 (W.I.36) are summarized in Table 8. 
The ISAGEX data set added three stations to the WEST 34 system. Station 9004 
(SAN FERN), though an addition to WEST data, is already included in the WNl4 
solution. Due to the small number of ISAGEX observations, only a minor improve- 
ment could be gained by the addition of the ISAGEX data. Their positive influence 
is reflected in the standard deviations, particularly at the stations with ISAGEX 
observations . hi Table 17 a three parameter transformation is shown between 
WEST 34 and WEST-ISAGEX 36. The residuals for some stations amount to the magni- 
tude.of their standard deviation, but there is no significant shift in the origin. A, glance 
at Table 8 further shows that a number of stations could not be determined accurately. 
The coordinates of stations 8705 (BRDUX), 8712 (OPICI), 8713 (ORIAA) and 8714 
(SRDIN) still exhibit extraordinarily large standard deviations. This is an imme- 
diate result of the increased variances in the observational data (see Table 4 
[Ehmsperger, 1974] for the standard errors stationwise, as obtained by the smoothing 
procedure.) The adjusted coordinates of 8711 (CATAN) could also be expected to 
exhibit a large standard deviation according to Table 4 [Ehrnsperger , 1974], but 
because this station has been connected to the nearby WN14 station, 6016 (CATANIA), 
by relative constraints, it does not exhibit a large variance. The large standard 
deviation of station 8722 (REKVK) is due to unfavorable geometric conditions and 
scarcity of observations. Similar reasoning would apply to stations 8716 (MADRD) 
and 8720 (LOVOA). A general characteristic of such stations is that the tie lines 
are ill -distributed (i.e. , they are not well spread around 360°). This problem, 
of course, is unavoidable in a non-global satellite triangulation. But it can also be 
seen that the standard deviations of coordinates of common stations have decreased 
in many cases by a factor of one-half as compared to the corresponding WN14 values. 

It is seen that the ISAGEX data add much strength to the station 8034 (DELFY), 
which is connected by relative constraints to the nearby WN14 station 8009 (DELFTH). 
Adjustments were performed applying no constraints to either of these stations (i.e., 
no station, height or relative constraints). Using only the WEST observations in the 


^8 


adjustment, it is found that station 8009 (DELFTH) shows no significant change in 
either the coordinates or the standard deviations; while the standard error in 8034 
(DELFY) rapidly increases to 35 m in each direction and its coordinates also change 
significantly. In contrast, when the combined data is used, neither of the two stations 
^ows aty significant change in coordinates or standard deviations . In fact, station 
8034 (DELFY) adjusted to the same position as obtained when using all constraints 
(especially in the case of relative constraints). 

The consistency of the stations 8721 (ONDREJOV) and 8723 (POTSDAM) can 
be estimated from comparisons with various solutions. Two solutions have been 
published based on the geometric and dynamic data of the ISAGEX campaign. A 
coordinate comparison is given in Table 9. If the origin of the coordinate system 
of the solution in [Marsh, et al., 1975] is taken as being at the geocenter, these 
coordinates can be transformed to the WN14 system by adding the shifts AU = - 21m, 
AV = - 5m and A W = 2 m. A similar comparison is possible with the coordinates 
given in [Gaposchkin, et al. , 1975]. But for computations here, the transformation 
parameters from Table 27 are used to transform the solution. A further attempt 
at comparison can be made using ED50 coordinates. The ED50 coordinates for 
8723 (POTSDAM) are taken from [Bracket, 1973]; while those of 8721 (ONDREJOV) 
from [Marsh, et al., 1975]. Marsh’s height presumably refers to an ellipsoid 

centered at the geocenter. But in the comparison here it is necessary to have the 
actual hei^t above the ED50 ellipsoid which is derived from MSL height as given 
in [Karsky, et al. , 1974], plus the undulation referring to the Potsdam geoid. The 
result after transforming to the W.I.36 solution, is given in Table 9, which seems to 
confirm the weak determination of the two stations involved. This is a consequence 
of both the lack of observations and their di-distribution. Marsh’s and Gaposchidn's 
values for station 8723 (POTSDAM) seem to be in relatively good agreement. 
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Table 8 


w 

o 


Summary of Cartesian Coordinates 
(Solution WEST 34 and W. I. 36) 


STATION 
NO I NAME 


6006 

6016 

6065 

8004 

8009 

8010 
8011 
8015 


8016 

8019 * 

8030 

8031 

8032 

8033 

8034 
8609 


8701 

8702 

8703 

8705 

8706 
8710 ' 

8711 

8712 


8713 

8714 

8715 

8716 

8717 

8718 

8719 

8720 


8721 

8722 

8723 
9004 


TROHSO 

CATANIA 

PEISEN 

BRNSG 

DELFTH 

ZIMLD 

MALVRN 

HAUTE PROVENCE 


STRBG 

NICEFR 

MUDDNI 

E0N5G 

HOPBG 

FRNFT 

DELFY 

ST.MICHEL FR 


GRAZA 

BRXOR 

COPHN 

BRDUX 

GOULT 

WSNDF 

CATAN 

OPICI 


ORI AA 
SRDIN 
TANIA 
KADRD 
MAORI 
KLSRH 
CATNA 
LOVOA 


ONDREJOV 
REKVK 
POTSDAM 
SAN FERNANDO 


o w 
> 




o 



SOLUTION 

WEST -34 




1 


SOLUTION 

WEST-1 SA6EX-36 



u 

V 

W 

Clu 

Civ 

cr 

1 

u 

V 

W 

CTu 

Q 

c 

cr w 

2102929 .0 

721666.4 

5958181,0 

2,2 

2.7 

2.4 

1 

1 

210292® .1 

721666.4 

5958182.4 

. 2.1 

2.6 

1.0 

4896308,5 

1316171.0 

3056668,1 

1.7 

2,3 

2,2 

I 

4896388,8 

1316171,2 

3856668.6 

1.3 

2.2 

1.6 

4213563 .9 

820824.8 

4702705,6 

l.o 

2,2 

2,1 

1 

4213564.0 

820824.1 

4702786.5 

1 «6 

2.1 

1.5 

3810505.9 

700050,1 

5042640.0 

5.0 

5,1 

4,7 

1 

3818505,2 

708049 ,2 

5042640.2 

4.9 

5.1 

3.7 

3923399,5 

299883,0 

5002971,2 

4.2 

5,2 

4,0 

I 

3923398 .0 

299876,5 

5002970.7 

3.2 

3.7 

2.5 

4331299,3 

567503.6 

4633113.2 

3.5 

4,9 

3,6 

1 

4331298.3 

567502 .3 

4633114.5 

2.4 

3.6 

2.2 

3920172.4 

” 134744,0 

5012722,2 

3.1 

2,6 

3.4 

1 

3920169.5 

- 134745,0 

5012719.5 

2.7 

2.5 

2.2 

4578317.8 

457962,2 

4403193,0 

3,4 

4,3 

3.6 

1 

4578322.3 

457962.3 

4403189,3 

2.1 

3.4 

2.2 

4108648,1 

571418,2 

4760143.7 

3.2 

4 . 1 

3.8 

1 

4188647.8 

571417.7 

4760144.8 

3.1 

4,0 

2.8 

4579462,5 

586587.1 

4386418.7 

3.6 

6,9 

3,8 

1 

4579462,5 

586586.0 

4306418.7 

2-9 

6.7 

3.1 

4205629.7 

163702.5 

4776538,9 

4.7 

7,9 

4.1 

! 

4205629,0 

163700.7 

4776538.8 

4.2 

7.8 

3.7 

3593859.0 

” 202757.0 

5248060,5 

9.4 

6,4 

7,0 

1 

3593859,6 

” 202757.4 

5248061.7 

9.5 

6,5 

6.4 

4213543.1 

820771,4 

4702810.6 

1-.9 

2,2 

2.1 

1 

4213543.2 

820770.8 

4702811.6 

1.6 

2.1 

1.5 

4041064.4 

620630,4 

4878629,4 

4.9 

6.0 

4.6 

1 

4041862.8 

6206 29,9 

4878628.7 

4.6 

5,9 

3.9 

3919690.4 

2.9 0829.4 

5005897,1 

4,2 

5.2 

4,1 

1 

3919688 ,9 

298822.9 

5005896.5 

3.2 

3.7 

2.5 







1 

4578365,3 

457946,5 

4403145.3 

2.1 

3,4 

2.2 

4194425,7 

1162696,8 

4647199.6 

3,9 

2.2 

4.2 

1 

4194427,3 

1162697.0 

4647204 ,9 

3,6 

2.1 

3.3 

4027918 .9 

307002,0 

4919436,3 

4.8 

7.7 

5.3 

1 

4027918,2 

307000.1 

4919436,5 

4.7 

7.7 

4.1 

3513633.9 

778936.2 

5248194.0 

7.0 

6.2 

5.0 

1 

3513633,0 

778935,3 

5248193.9 

6.8 

6.3 

4.3 

4530506.3 

” 41732,3 

44743 B 2.1 

17.4 

33.2 

18,0 

1 

4530507,0 

” 41734.5 

4474381,2 

17.6 

33.6 

17.9 

4587807.3 

419519,2 

4396450,5 

3.4 

4,3 

3.6 

1 

4507891,8 

419519,2 

4396446,8 

2.1 

3.4 

2.2 

3818504.1 

708051,4 

5042640.1 

5,0 

5.1 

4,7 

1 

3818503,4 

708050*5 

5042640.4 

4,9 

5,1 

3.7 

4096386.9 

1316170.5 

3856670.2 

1.8 

2.3 

2.2 

1 

4896387.2 

1316170,8 

3856670,8 

1,3 

2.2 

1.6 

4335518.6 

1063083,4 

4540932.6 

17.6 

30.1 

19-4 

1 

4335518,0 

1063083,8 

4540933.2 

17,7 

30.5 

19.5 

4628609 .9 

1471955.7 

4120468,5 

19,5 

39-2 

27,3 

1 

4628610,1 

1471955,8 

4120468.2 

19.6 

39,7 

27,6 

4805403.5 

784066,5 

4011526.6 

23.6 

29,5 

32.3 

1 

4885403.0 

784067,2 

4011527,4 

23.9 

29,8 

32.6 

4896390.9 

1316178.8 

3856662.4 

1.8 

2,3 

2,2 

f 

4896391 ,3 

1316179,0 

3856662.9 

1.3 

2.2 

1.6 

4850674.5 

- 315907,3 

4116627.0 

8,1 

7.5 

10,0 

! 

4850673,5 

” 315908.1 

4116626.2 

7.9 

7.6 

9,7 

4850674.6 

- 315910.0 

4116626.7 

8,1 

7.5 

10,0 

1 

4050673,6 

” 315910,8 

4116626,0 

7,9 

7.6 

9.7 

4146533.6 

613109,8 

4791487,6 

6.1 

6,5 

5.6 

1 

4146532.1 

613108,8 

4791487,6 

5.8 

6 i»A- 

5.1 

4683057,0 

1306097.9 

3879629,2 

1.8 

2.3 

2,2 

1 

4883057,3 

1306098,1 

3879629,7 

1.3 

2,2 

1.6 

3104204,5 

998359,4 

5463280.1 

10,5 

8.6 

5,6 

! 

3104203,0 

998359,4 

5463280.4 

10-4 

8-6 

5,1 







1 

3978430,8 

1051033.1 

4857553.9 

6.6 

6 • 8 

3.9 

2591994,3 

- 1078495.2 

5707863.4 

10.2 

18.3 

5.8 

1 

2591995,3 

- 1078494,5 

5707863.0 

10.2 

18,5 

4.9 







1 

3800631-9 

881944.2 

5028839-9 

6.8 

12.3 

4.3 







I 

5105579-8 

- 555246.5 

3769675.4 

2.8 

4.8 

3.5 


S § 


h---* 

a J f. 


Table 9 

Coordinate Comparison for ONDREJOV and POTSDAM 


COORDINATE DIFFERENCES 

Station 

Au(m) 

Av(m) 

Aw (m) 

8721 

-11.67 

16.30 

3.56 

8723 

-35.22 

37.25 

49.28 

[Marsh, et al. 

,1975] - W.I.36 

8721 

11.83 

-22. 33 

20.68 

8723 

-44.03 

29.41 

35.00 

[Gaposchkin, et al. , 1975] - 

W.I.36 

8721 

-42.14 

93, 69 

14.89 

8723 

-18.82 

39.96 

8.72 


ED 50 

- W.I.36 



A summary of transformation parameters for solution WEST-ISAGEX 36 with 
respect to various other systems is given in Table 10. The results of the individual 
transformations reflected in this summary table, are presented in Tables 17 and 21 
through 32. In all cases the standard Molodenskii and Bursa transformation models 
have been used. In Appendix III a detailed discussion on the differences between the 
two models is given. Tables 18 through 20 present the results of transformation 
between WEST 34 and other systems. 

The transformation with respect to the common stations of the W1SI14 system 
indicates no significant shift or rotation for the W. I. 36 solution as intended. The resi- 
duals in the transformation between WNl4 and W.I.36, however, show systematic 
corrections for the coordinates of various stations which are especially significant 
in the y coordinate (Table 29;). The newly derived translation parameters for the ED50 
(Table 22), with respect to the geocenter, are (Molodenskii transformation model): 

- 94.1 + 21 « 73m 

- 115.5 + 5 -110m 

- 125.3 - 2 fa 123m 

where 2lm, 5m and - 2m are the coordinates of the WN14 origin with respect 
to the geocenter. The translation parameters are in good agreement 
with the other solutions, giving appropriate consideration to 
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Table 10 


Transformation Parameters between Other Systems and WI36 


System 

Trans- 

formation 

Parameter 

ED50 

NGS^ 

SEf 

WN14 

GEM’lf 

WEST 34 

, 25 Stations 

3 Stations 

7 Stations 

10 Stations 

4 Stations 

24 Stations 

AU(mf 
AV (m)s 
AW(m)s 

-123. 9± 18.8 
-112. 2± 30.6 
-152. 7± 18.8 

Bursa T 

6.2± 20.0 
-93. 0± 77.9 
9.0± 10.3 

ranslation Par 
-14. 7 ± 18.0 
-31, 7± 27.6 
21. 1± 18.7 

ameters* 
2.4± 10.0 
13. 0± 21.8 
6.2± 8.9 

-40. 2± 8.8 
-33. 4± 15,2 
11. 6± 9,3 


AU(m) 

AV(m) 

AW(m) 

- 94. 1± 3.4 
-115. 5± 3.3 
-125. 3± 3.5 

Molodensk 
- 9.2± 2.7 
-lO.Oi 3.3 
16. 0± 3.5 

ii Translation 
- 8.2± 2.5 
-15. 3± 2.9 
9.7± 2,6 

Parameters 
0.2± 1.7 
0.7± 2.4 
0.1± 1.9 

-19,5± 2.0 
-30. 7± 2.2 
6.1± 2.0 

■ 

BwlH 

Hl9i 

6.12± 2.77 
0.70± 0.77 
-0.15± 0.65 
-0.17± 0.86 

0.90± 1.62 
~1.68± 1.58 
0.55± 0.53 
2.10± 2.13 

-0.32± 2.45 
-0.23± 0.68 
-0.38± 0.67 
0.53± 0,78 

-1. 19± 1.42 
0. 04± 0.50 
-0. 12± 0.31 
-0.47± 0,63 

1.56± 1,35 
-0.66± 0,37 
-0.73± 0.31 
-0.53± 0,45 


v("y 

,0.63± 0.58 
-0.18± 0.64 
-0.32± 1.00 




■ 


Factor® 

1.03 

9 

1.07 

2.10 

m 

1.00 

2.40 

1.04 

0.09 



’’Coordinates for common stations from [Schmid, 1974]. 

® Coordinates for common stations from [Gaposchkin, 1974], 

® Coordinates for common stations from [Lerch, et al., 1974], 

^See Appendix III- 

®If (geooenter-datum) is sought, add to the tabulated values of ^U, AV, AW the respective 
quantites 21m, 5m and -2m. 

i/), € when positive, represent counterclockwise rotations about the respective W,V, U 
axes, as ■viewed from the end of die positive axis. 

’The positive directions of the axes considered at the geodetic initial point are south, east , 
aiid along the ellipsoidal normal outwards with the respective rotations about these axes 
as T} , C • o • 

® Factor by which the variances of the systems had to be scaled in order to obtain the 
variance of unit weight Oo close to unity. 

® Variance of each coordinate in the European datum is I44m® . 
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the variances. However, these do not have geometrical significance as explained in 
Appendix in. Also see [Anderle, 1974 and Weightman, 1975]. A graphical repre- 
sentation of the residuals is given in Figures 4, 5 and 6. 

The transformation parameters with respect to the solutions NGS and GEM 6 
shoiild not be overemphasized since there is a lack of an adequate number of common 
stations. 

Height constraints have been used very successfully in the WN14 ^obal solution. 
The use of weighted height constraints provides a unique tool to select the scale to 
fit some criterion. In a non-global satellite triangulation, as in the present case, 
hei^t constraints are only used to strengthen stations at the edge of the triangulation 
area, or stations with insufficient observations. Solution WEST 34 has been repeated 
without the use of any height constraints (WEST 35) and is given in Appendices I and n 
along with the transformation between WEST 35 and WEST 34. There is no significant 
shift between the systems, although a rotation around the W and U axes seems to be 
present. The height residuals are large for all weaMy determined points. But also 
two supposedly well determined stations, 8031 (EDNBG.) and 8701 (GRAZA), exhibit 
a large residual in height. 

The scale factor for the European datum was sought with respect to the satellite 
system by using the original ED50 coordinates instead of the revised coordinates as 
explained in [Weightman, 1975]. He indicates the possibility that current publications 
give revised coordinates for the terminal station of the European, base lines, which 
are based on the base line adjustment. This seems to be partly confirmed in Table 
11, which shows the difference between die lengths, as computed from ED 50 coordi- 
nates (Table 1) and values for the base lines as obtained from traverses. 


Table 11 

Comparison of Base line Length and ED Chords 


STATION 

ED 50 minus traverse lengtl 

1 From 

To 

(m) 

ppm 

6006 

TROMSO 

6016 

CATANIA 

2.36 

0.67 

6006 

TROMSO 

8032 

HOPBG 

0. 39 

0.16 

8032 

HOPBG 

8711 

CATAN 

-10.44 

9.09 

8011 

MALVRN 

8032 

HOPBG 

-1. 04 

1.00 

8032 

HOPBG 

8701 

GRAZA 

0. 24 

0.69 
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Latitude Differences After Transforming EOSO to W.L36 ( in Meters ) 
7 Parameter Tronsfermotion (WL36>ED50) 



Fig, 4 


34 



Longitude_Oiffereflces After Tronsforming E050 to W,L36 ( in Meters) 
7 Parameter Trensfonmtion (VKZ.3S«ED50) 1 
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The conclusions drawn from Table 11 are, of course, not based on the best 
evidence. But due to lack of additional information, it is assumed that all but the 
coordinates of station 8711 (CAT AN) are revised coordinates. The originality of 
the CATAN coordinates seems to be confirmed when computing the distance between 
6016 (CATANIA), as given in [NASA, 1973] and 8711 (CATAN), as given in [Ehrn- 
sperger, 1974] . The computed distance is 14, 59m as opposed to the 2, 76 m given 
in Circular Letter No. 37 [WEST, 1966-1972]. The difference of -10. 44m in Table 
11 might also be interpreted such that the southern terminal point of the base line is 
8715 (TANIA), since this station is approximately 12m south of 8711 (CATAN) [WEST, 
1966-1972]. But comparing the adjusted chords and their measured values (Table 12), 
it can be seen that this interpretation is probably wrong. Notice that the last two 
chords which are actually parts of the chord 6006 (TROMSO) - 6016 (CATANIA), 
were not constrained in the W.I. 36 solution. 


Table 12 

Comparison of Adjusted Chords 


LINE 

Given Length 
(m) 

-Adjusted (W.I. 36) 
ppm 

6006-6016 

-1.32 

0.37 

8011-8032 

-3.68 

0.35 

8032-8701 

-.32 

0.92 

6006-8032 

0.75 1 

0.30 

8032-8711 

-0.48 

0.40 


Thus it was concluded that station 8711 (CATAN) is the only base line terminal 
station which can be used in subsequent scale factor investigations. In addition, those 
stations which are wealdy determined were also excluded for the same considerations. 
In smnmary, the following stations were excluded: 


6006 TROMSO 
8032 HOPBG 
soil MALVRN 
8701 GRAZA* 
8705 BRDUX* 
8712 OPICI 


8713 ORIAA 

8714 SRDIN 
8716 IVIADRD 
8720 LOVOA 
8222 REICVK* 

* See additional comments in Table 1. 
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The remainir^ station coordinates were used to perform a standard seven 
parameter Molodenskii transformation resulting in a scale factor of 1, 74± 3, 70 ppm. 
Compared to the previous value of 6.12± 2,77ppm, there seems to be a significant 
change, which gave reason to further investigate the individual stations. After com- 
paring chords it was found tliatthe chord differences starting from 8031 (EDNBG) 
were systematically different from others . Deleting this station and repeating the 
transformation, a scale factor of 5.92± 3. 54 ppm was obtained. This result clearly 
shows the weak determination of the scale factor and gives support to the high 
variance attached to it. Table 13 gives the chord differences between W.I.36 
and ED 50 systems. ' It is readily seen that there is a great- dispersion , in 
the scale factor. . .. • • 

A similar result was found by [Marsh, et al. , 1971]. He also investigated 
specific areas within Europe. According to [NASA, 1973], the adjustment of the 
Central European Network was carried out between 1945 and 1947. This triangulation 
network covers the region that lies from 47° to 56° North latitude and between 6° 
and 27° East longitude. The Central European Network was extended by the addition 
of two separate adjustments of large networks of triangulation known as the South- 
western Block and Northern Block. The Central Network was substantially held 
fixed and, with the addition of the two blocks, forms the European Triangulation 
based on v/hat is now designated as the European datum. Tables 14, 15 and 16 
are modest attempts to represent the scale factor within these three blocks. 

Differences between Satellite and Groimd Survey 
Chord Lengths (meters) 

(Satellite - Ground Survey) 

Table 14 
Northern Block 



8009 

8030 

8031 

8030 

1.5 



8031 

m 

-30.1 


8702 

- 7.3 

■B 

-45.1 


39 













Table 15 
Central Block 



8004 

8016 

8033 

8703 

8016 

7.9 


HH 


8033 

7.5 

0.5 

HH 


8703 

2.6 

10.7 

16.2 


8718 

12.3 

-0.5 

6.6 

15.1 


Table 16 

Southwestern Block 



8010 

8019 

8706 

8711 

8019 

- 3,3 




8706 

-0.5 

-14. 6 



8711 

18.2 

33.2 

21.0 


9004 

m 

-0.3 

7.4 

45.2 


From the above three tables the following conclusions might 
be in order: 

1) In both the Northern and the Southwestern Blocks there seems to be 

a lai^e variation in chord length differences. The differences tend to be generally 
negative in the Northern Block. 

2) in the Central Block the variations are relatively small and the differences 
are generally positive. 


In conclusion, the above analysis has yielded slightly improved values of 
coordinates of the European stations in the \\TSri4 solution. Two stations could 
be added to the WN14 solution. In the absence of any statistics, the use of the 
seven image data in analysis has been deferred until the reduced data for all the 
WEST plates is available. At that stage the whole data maj^ merit fresh attention. 
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Table 17 

Transformation 


WEST 34 -TO- ( 3 PARAMETER) 

43^ 444444c 4414 4444 4 444 444>4<4>ti4 4444 44 44 4444 


SOLUTION FOR 3 TRANSLATION PARAMETERS 
(UNITS - METERS) 


(USING variances only) 


1 

t 


i 

i 

I 

I 


DU 

METERS 

0.27492E92P-C1 
♦ 0«,30 


DV 

METERS 

- 0 . 508634770+00 
+ 0,35 


D\V 

.METERS 

0r21750O85D+0n 

± O.?o 


i 


VARIANCE - COVARIANCE MATRIX 


M02= 0.09 


0,876958800-01 

0.0 

o 

. 

o 

c 

* 

o 

0.12217624D+00 

0.0 

o 

. 

o 

c 

. 

o 

0.867643660-01 



COEFFICIENTS 0^ CORRELATION 


O.lOOOOCOOD+01 

0.0 

0.0 

0.0 

o.ioroooonn+oi 

0,0 

c 

. 

c 

C 

• 

C 

0, looononoD+ni 
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Table 17 (Continued) 


‘RESIDUALS V 



VII 1 

WEST 34 ) 



V 2 ( W . 

1,36 ) 


VI 

- V ? 


60 C 6 

0.0 

0,3 

0.3 

6006 

“ 0 . 0 

- 0.3 

- 0.1 

0.1 

0,6 

'^,4 

8 CC 4 

•> 0.3 

“ 0.2 

0.0 

8004 

0.3 

0.2 

“ 0.0 

- 0.7 

” 0,4 

0,0 

8 CC '> 

“ 1.0 

“ 3.9 

- 0.5 

8 009 

0.6 

2.0 

0.2 

“1 .5 

- 5.9 

“ 0.7 

8010 

- 0.7 

“ 0.4 

0,8 

8010 

. 0. 3 

0.2 

- 0.3 

- 1,0 

- 0.6 

1.1 

eon 

- 1.6 

“ 0.2 

- 2.1 

8011 

1.3 

0.2 

0.9 

“ 2.9 

- 0.4 

- 3.0 

8016 

- 0.2 

f'.O 

0,6 

8016 

0.2 

- 0.0 

- 0.3 

“ 0.3 

0.0 

0,9 

8019 

- 0.0 

“ 0.2 

“ 0.1 

6019 

0.0 

0.2 

0 , 1 

“ 0.0 

- 0,4 

- 0,1 

8030 

“ 0.4 

“ 0.6 

” 0.2 

8 C 30 

0 . 3 

0.6 

0.2 

“ 0.7 

“ 1.2 

-0. 4 

BC 31 

0.3 

O.I 

0.6 

8 031 

“ 0,3 

“ O.i 

“ 0.5 

0.6 

0.1 

1.0 

8032 

0.0 

- 0.0 

0.5 

8032 

- 0,0 

0.0 

- 0.3 

0.1 

- 0,0 

0,8 

8033 

“ 0.9 

0.0 

“ 0.5 

8033 

0.8 

-0 . 0 

0,4 

“1 .7 

0,1 

- 0.0 

8 701 

0,9 

0.4 

3.1 

8701 

- o,ft 

- 0.4 

” 1.9 

1 .6 

0.8 

5,0 

8702 

- 0.4 

“ 0.7 

- 0.1 

6 702 

0 ,^ 

0.7 

0,0 

- 0.7 

-1 .3 

- 0.1 

8703 

“ 0.5 

- 0.2 

“ 0.2 

8703 

C.4 

0.2 

0.1 

“ 0.9 

- 0,3 

- 0,3 

8705 

0.4 

“O.e 

“ 0.5 

8705 

- 0.4 

0 . 8 

0.5 

0.7 

” 1.7 

- 1,1 

8706 

3.2 

0.4 

“ 2.9 

8 706 

“ 1.3 

- 0,3 

1.0 

4.5 

0.7 

-^.9 

8711 

0.2 

0.4 

0.2 

871 1 

“ 0 . 1 

- 0.4 

- 0.1 

0,3 

0.8 

0.4 

8712 

“ 0.3 

0.5 

0.2 

8712 

0.3 

- 0,5 

“ 0 . 2 

- 0.7 

1.0 

0 . 3 

8713 

0 . 1 

0.3 

- 0.3 

8713 

- 0 . 1 

- 0.3 

0,3 

0.2 

0.6 

- 0 . 5 

8716 

- 0.3 

0.7 

0.3 

8714 

0.3 

“ 0.7 

- 0 , 3 

- 0.6 

1 .3 

0.5 

8716 

- 0.5 

“ 0.1 

“ 0.5 

8716 

0,5 

0 , 1 

0 . 5 

“ 1.0 

“ 0.2 

- 1.0 

8718 

“ 0.8 

“ 0.2 

” 0.1 

8 718 

C . 7 

0.2 

0 . 1 

“) .5 

“ 0.5 

- 0.3 

8720 

- 0.8 

0.3 

0.0 

8720 

0 . 7 

“ 0.3 

“ 0.0 

- 1,5 

0,6 

0.0 

8722 

0.5 

0.6 

- 0.4 

8722 

“ 0.5 

“ 0,6 

0,3 

1.0 

1.2 

- 0,6 


UNIT OF RESIDUALS ( METERS ) 


42 



Table 18 

T ransf ormation 


ED-50 -TO- WFST 34( 3 PARAMETER) 


SOLUTION FOR 3 TRANSLATION PARAMETERS 
(UNITS - METERS) 


(USING VARIANCES ONLY) 



DU 

DV 

DW 


METERS 

METERS 

METERS 

« 

-0.9352 87.68D+02 

-0, 11079R06D+03 

-0.12670269D+0? 


+ 2.P8 

+ 3.00 

+ 2,86 


VARIANCE - COVARIANCE MATRIX 


HO 2= 1.03 

0.83177967D+01 

0.0 

0.0 


0.0 

0. 897993710+01 

0.0 


0.0 

0.0 

0.81659603D+01 


COEFFICIENTS OF CORRELATION 


0. iooooooon +01 
0.0 
0.0 


0.0 

0. lOOOCOOOD+01 

0.0 


0,0 

0.0 

O.lOOOOOOOD+Ol 
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Table 18 (Continued) 


6006 

80 (vt 

P009 

8010 

8011 

8016 

8019 

8030 

8031 

8032 

8033 

8701 

8702 

8703 

8705 

8706 
8 711 

8712 

8713 

8714 
8716 
8718 
8720 
8722 


RESIDUALS V 


rOTODUCIBIUTY OF the 
OUIGINAE PAGE IS POO- 


VI ( 1 

ED “50 

) 


V 2 ( WEST 34 ) 


VI 

“ V 2 


“ IP . 3 

29.1 

6.0 

6006 

0.6 

“ 1,4 

“ 0.2 

“ 18,9 

30.5 

6.2 

“ 3.3 

0.8 

1.4 

8004 

0.6 

“ 0.1 

“ 0.2 

“ 3,9 

0,9 

1.7 

7.5 

“ 5.1 

3.5 

8009 

“ 0.9 

1.0 

“ 0.4 

8.4 

“ 6.1 

3.9 

1.3 

“ 14.7 

2.7 

8010 

- 0.1 

2.5 

“ 0.2 

1.4 

- 17.2 

3.0 

14.4 

“ 4.1 

“ 0.4 

8011 

“ 0.9 

0.2 

0.0 

15.3 

~4 • 3 

- 0.4 

2.1 

1.0 

“ 4.1 

8016 

- C .2 

- 0,1 

0.4 

2.3 

1.1 

“ 4.6 

“ 1.4 

“ 19.3 

5 .P 

8019 

0 . 1 

6,4 

“ 0.6 

“ 1,5 

“ 25.7 

6.4 

4.7 

“ 15.1 

4.3 

8030 

“ 0.7 

6.6 

“ 0.5 

5.4 

“21 .7 

4.8 

6.2 

16.2 

" 6.9 

8031 

“ 3.9 

“ 4.7 

2.3 

10.1 

20.9 

“ 9.3 

“ 9.5 

9.5 

13.7 

8032 

0.2 

“ 0.3 

“ C .4 

“ 9,7 

9.8 

14.1 

2 . 5 

“ 1.3 

“ 1.5 

8033 

“ 0.4 

C .3 

0.2 

2.9 

- 1.6 

- 1,7 

0.8 

5.9 

“ 9.7 

8701 

- 0. 1 

“ 0.2 

1.2 

0.8 

6.1 

“ 10,9 

“ 1.6 

” 18.3 

1.7 

8702 

0.2 

7.5 

“ 0,3 

“ 1.8 

“ 25.8 

2.0 

“ 4.3 

3.0 

2.2 

8703 

1.5 

“ 0.8 

“ 0.4 

“ 5.8 

3.9 

2 . 6 

“ 0.3 

2.1 

“ 0.2 

8705 

0 . 6 

“ 16 . 0 

0.4 

“ 0.9 

18 . 1 

“ 0,6 

“ 4.1 

” 10.9 

4.2 

8706 

0.3 

i .4 

“ 0.4 

- 4.4 

“ 12.3 

4.5 

“ 10.9 

17.4 

2.4 

8711 

0 . 2 

“ 0.6 

“ 0 . 1 

“ 11.1 

18.0 

2.4 

“ 0.5 

5.4 

“ 2.6 

8712 

1 . 1 

“ 33.7 

6.8 

“ 1.6 

39.1 

“9 . 4 

“ 4.5 

“ 1.7 

3.8 

8713 

1 1.9 

18.5 

“ 19.9 

“ 16,4 

“ 20,2 

23.7 

“ 8.2 

11.4 

3.7 

8714 

31.7 

“ 68.7 

“ 26.9 

“ 39.9 

80.1 

30,6 

27.8 

“ 21.1 

“ 29.3 

8716 

“ 12.6 

8,3 

20,2 

4 0.4 

“ 29.4 

- 49.6 

6.0 

6.3 

“ 6,3 

8718 

“ 1.6 

“ 1.9 

1.4 

7.6 

8.2 

“ 7.7 

“ 2.2 

“ 1.4 

1.4 

8720 

1.7 

0.7 

“ 0.3 

- 3.8 

“ 2,2 

1.7 

“ 4.2 

5,0 

4.2 

8722 

3.0 

“ 11.6 

“ 1,0 

“ 7,2 

16.6 

5,1 


UNIT OF RESIDUALS { METERS ) 
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Table 19 



Transformation 



ED-50 -TO- WF«^T 34( MOL MOnPL ) 

SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PARAMETERS 


(USING VARIANCES ONLY) 


DU 

DV 

OW 

delta 

OMEGA 

PSI 

EPSILON 

METERS 

METERS 

MFTFRS 

1X1, 0+6) 

SECONDS SECONOS 

SFCONOS 

-94.12 - 

114.17 -125.66 

6.17 

0.3? 

“0.16 

0r09 

+ 3.50 

± 3.28 

+ 3.5P 

*. 3.02 

± 0.8 6 

+ 0.73 

+ 0.9R 


VARIANCE - COVARIANCE MATRIX 


S02= 1.03 

C.123D+02 “C.234D + 00 0.212D + C1 -0.217D-05 0.725D-05 -0.3«>?D-05 -0.4650-05 

-0.23AD+00 0.1070+02 -0.398D + 00 0.3120-05 0.131D-05 0.22SDr06 0.22RD-05 

0.212D+01 “0.3980+00 0.129D+02 0.1770-05 0.4540-05 -0.369D-05 -0.9510-05 

-0.2170-05 0.3120-05 0.1770-C5 O.^lAD-ll -0.3840-12 0.160D-12 -0.48-»D-12 

,0.7250-05 0.131D-05 0.4540-05 -0.3840-12 0.1720-10 -0.3210-11 -0.1030-10 

-0.392D-05 0,228D-Cft -0.3690-05 C.16C0-12 -0.3210-11 0.125D-10 0.4120-ii 

-0.465D-05 0.228D-C5 -0.9B1D-C5 -G. 4870-12 -0,1030-10 0,4)20-11 0.7?Bn-10 


CnFFFTCIENTS OF CORRELATION 


O.lOCO+01 -0.2040-01 0.1680+00 -0,2040+00 0.44on+00 -0.3160+00 -0,27«0+00 

-0.2040-01 0. 1000 + 01 -0.3380-01 n,314D+00 0. 0620-01 0.1970-01 0.1460 + 00 

0.168D+00 -0.3380-01 O.lOOO+Ol 0,1630+00 0,3050+00 -O,?9l0+O0 -0,57?^+00 

-0,2040+00 0,3140+00 0,1630 + 00 0,1000+01 -0.3060-01 0.14OD-01 -0.3370-01 

0.499D+00 0 . 9620-01 0.3050 + 00 -0,3060-01 0.1000+01 -0.2190+00 -0,51'?D+OO 

-0.3160+00 0.1970-01 -0,2910+00 0.1490-01 -0.2190+00 O.lOOD+01 0,2450+00 

-0.278D+00 0.1460+00 -0.5720+00 -0. 3370-01 -0.519D+00 0.245D+C0 0. 1000+01 

NOTE : THE POSITIVE AXES ARE TOWAROS MEAN 

GREENWICH, EAST, AfJO CIO. THE ROTATIONS 

PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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Table 19 (Continued) 


reproducibiuty of the 

ORIGINAL PAGE IS POOR 

RESIDUALS 


VU ED-SO J V2< WEST 34) VI •> V2 


6006 

-8.1 

24.6 

-1 .9 

6006 

0.3 

-1.2 

0.1 

“8 . 4 

2r?.8 

-2.0 

P004 

-2.7 

“0,5 

0.4 

8004 

0.5 

0. 1 

“0,1 

“3.2 

-0.6 

0.5 

80C9 

8.2 

“4.1 

2.6 

80C9 

-l.C 

0.8 

“0. 3 

9.2 

"4 . 9 

2.0 

8010 

-0.5 

“14.4 

4.3 

8010 

C.O 

2.4 

“0.4 

-0.5 

-16.0 

4,7 

8 011 

15.7 

-0,4 

— 1 .6 

8011 

-1.0 

0.0 

0.1 

16.8 

-0.4 

-1.7 

8016 

1.1 

1.0 

-? ,4 

801 6 

-0.1 

“0,1 

0.3 

l.I 

1.1 

“3.8 

8019 

-4.3 

“18.8 

C.o 

£019 

0,4 

6.2 

-0.0 

“4.7 

-25.0 

0,8 

8030 

4.1 

“13.3 

4.P 

8030 

-0.6 

5.8 

“0.6 

4.8 

-10.1 

5.4 

8031 

8.3 

19.0 

“0.2 

PC31 

-5. 1 

“5.5 

3.1 

13.5 

24.5 

”12.3 

8032 

“11.1 

8.1 

14. c- 

8032 

C.3 

“0.3 

“C.5 

“1 1.4 

8.4 

15,4- 

8033 

2.2 

”1.7 

“1.5 

8033 

“0.4 

0.4 

0,2 

2.5 

-2.1 

-1.7 

9701 

“1.0 

2.4 

“8.2 

8701 

0.1 

“0.1 

l.Tf 

-1.1 

2.5 

-0.2 

8 702 

-1.4 

“17.4 

1.2 

F702 

0.? 

7.1 

-0.3 

“1.6 

-24.5 

1.6 

8703 

“2.6 

1.0 

-G.O 

8703 

0.9 

“0.3 

0.0 

-3.4 

1.3 

-0, 1 

8705 

-0.9 

2.6 

0.6 

8705 

2.0 

“2C.0 

“1.4 

-2.9 

22.6 

2.0 

8706 

“7.0 

“0.4 

7.2 

8706 

0.5 

1.2 

-0.7 

“7.5 

-10.6 

7,8 

8711 

-16.7 

14.4 

9.4 

8 711 

0.4 

“0.5 

-0.3 

“17,1 

15.0 

9.7 

8712 

-1.3 

5.0 

-1.9 

8712 

2.9 

-31.4 

5.0 

“4.2 

36.4 

-6.0 

8 713 

-5.8 

“2.1 

4.7 

8713 

15.3 

22.7 

-24.4 

“21 .? 

-24.8 

79.2 

8714 

-9.3 

11.4 

4.4 

8714 

35.8 

“68.8 

-32.2 

“45.1 

80,2 

36.6 

8716 

25.6 

-16.2 

-26.4 

8716 

“li,6 

6.4 

18.2 

37.1 

-22.6 


8 718 

5.3 

6.1 

-5.0 

8718 

“1.4 

“1.8 

1.3 

6 .6 

7.0 

“7. 1 

8720 

0,3 

“4.5 

-2.0 

8720 

”0.2 

2.3 

0. 4 

0,5 

-6,8 

-2.5 

8722 

2.0 

7.2 

“1.5 

8722 

“1.4 

“16.7 

0.4 

3.4 

23.0 

-1.9 


UNIT OF RESIDUALS (METERS) 


* Besiduals in the Cartesian coordinate system. 
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Table 19 (Continued) 


reproducibility of the 
original page is poor- 


VIC EO-5C ) V2C WEST 34J VI " V2 


6006 

- 0.98 

25 . 8 ® 

- 1.7 

6006 

0.15 

- 1.24 

0.0 

- 1.13 

27.14 

- 1.7 

R004 

2.43 

- 0.03 

- 1.3 

8004 

— 0.41 

O.Pl 

0.2 

2.84 

“ 0.04 

“ 1.6 

P00«) 

-A .5 6 

- 4.75 

6.9 

8 CC 9 

0.56 

0.35 

- 0.8 

“ 5.12 

- 5.6 0 

7.7 

iolo 

4.63 

“ 14.25 

1.5 

8010 

“ 0.52 

2.40 

“C.O 

5.15 

- 16.66 

1,6 

8011 

- 13,40 

0. 12 

8.4 

8011 

0.88 

-0.01 

- 0.5 

- 14 , 2 Q 

0.14 

8.0 

8016 

- 3.16 

0.83 

“ 1.8 

8016 

0.29 

“ 0.10 

i ).2 

- 3.45 

0,93 

“ 2 rO 

8019 

11.04 

“ 18.07 

1.3 

8010 

-1 .47 

6.11 

0.2 

12.52 

- 24,10 

1.1 

8030 

0.4 A . 

- 13.44 

6.0 

8030 

- 0.07 

5.84 

- 0,7 

0.51 

- 19 . 2 P 

6 . 7 

8031 

- 11.14 

10.48 

- 3.5 

8031 

5.74 

- 5.77 

“ 0.1 

-16,88 

25. 25 

“ 3.4 

8032 

16.94 

10.07 

4.6 

PO ?2 

- 0.47 

" C .32 

- 0.2 

17.41 

10.38 

5.0 

8033 

- 2,39 

- 2.02 

0,1 

80 33 

0,36 

0.48 

- 0.0 

- 2.75 

- 2.50 

0.1 

8701 

“ 5.31 

2 . 6 ? 

- 6.2 

8701 

0.62 

“ 0.11 

0.8 

“ 5.93 

2 . 73 

- 7.0 

8702 

2.95 

“ 17.20 

- 0.7 

8702 

— fi .76 

7.09 

0.3 

3.70 

- 24.29 

- 1. 0 

8703 

1.86 

1.58 

- 1.3 

87 03 

- 0.65 

“ 0,47 

0.5 

2.51 

2.04 

- 1 . c . 

8705 

1.11 

2.60 

“ 0.3 

8705 

- 2 . 5 C ; 

“ 19,94 

0.6 

3.62 

22.54 

- 0 . 8 

8706 

10.56 

“ 8,74 

- 0.6 

8706 

“ 0 . 0 ? 

1.13 

0.0 

11 , 4 P 

- 9.87 

“ 0,7 

8711 

14.97 

18.28 

- 4,1 

8711 

“ 0.37 

- 0.59 

- 0.0 

15.34 

18.87 

“A. 1 

8712 

- 1.26 

5.17 

- 1.4 

8712 

6.85 

“ 31,15 

0.3 

“ 8,11 

36 , 3 ? 

“! c 7 

8713 

7.61 

- 0.27 

- 1,6 

8713 

- 32.55 

16,98 

0,5 

40,16 

- 17.25 

- 2 , 1 

8714 

8.08 

12.73 

- 2.9 

87 14 

- 40.40 

“ 73.56 

-1 .4 

4 R .48 

86,20 

“ 1.5 

8716 

“ 37.34 

- 14.5 1 

3.1 

8716 

21 .64 

5.60 

2.7 

“ 58,98 

- 20.10 

0.3 

8718 

- 8.44 

5.24 

- 0.4 

8718 

2 .04 

- 1.56 

-n. 1 

“ 10.40 

6,«o 

“ 0 .’ 

87 20 

- 0.10 

- 4.38 

- 2.3 

8720 

- 0.18 

2.25 

0.6 

0.09 

“ 6.63 

“ 2.9 

8 722 

0.17 

7.38 

“ 1,8 

8722 

“ 4.45 

- 16.00 

2.6 

4.63 

23.33 

- 4.4 


UNIT DF RESIDUALS (METERS) 

♦The spherical residuals represent the residuals in the Cartesian coordinate system as 
transformed into the curvilinear (<p, X, h) system. 
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Table 20 

Transformation 


ED-50 -TO- WF*-T 5A{ VF.I5 MOnFL ) 

SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 R^TATTON OARAMFTFRS 


(USING variances ONLY) 


DU DV DW delta ALPHA KSI ETA 

METERS METERS METERS (Xl.D+6) SECONDS SECONDS SECONDS 


“ 

94.12 -114.17 -125.66 

6,17 

0.2P -0 

cie “0. 

16 

♦ 

3.50 + 3.28 + 3.59 

± 3.02 + 

0.63 + 0 

.69 + 1 r 

17 


VARIANCE 

- COVARIANCE MATRIX 



S02 = 

1.03 





0.123D+02 

“0.2340+00 0.212D+01 

-0.2170-05 

0.2440-05 

-0.2770-05 

-0.8710-0 5 

>0.2340+00 

0.107D+02 “0.39RD + 00 

0.3120-05 

0.243D-05 

-0,2940-06 

0 .lOOD-05 

0.212D+01 

“0 .3980+00 0.129D+G2 

0. 1770 -05 

-0.274D-05 

-0,1380-05 

-0.1 10 0-04 

0.217D-05 

0.3120-05 0.1770-05 

0.91^0-11 

-0 .5720-12 

0.2660-12 

-0.113 0-12 

0,2 44D-05 

0,2430-05 -0.274D-05 

-0.5720-12 

0 .9390-11 

0 ,2300-12 

0,5830-12 

-0.277D-05 

“0 .2940-06 -0.1380-05 

0.2660-12 

0.239D-1? 

0.1120-10 

0.1620-11 

■C.8 71D-C5 

0.1000-05 -0.110D-C4 

-r. 1130-12 

0 .583D-1? 

0.1620-11 

0.3196-10 


COEFFICIENTS OF CORRELATION 



0. 1000+01 

-0.204D-01 0.168D+00 

-0.204D+00 

0.227D+00 

-0.236D+00 

-0.4400+00 

•0.2040-01 

O.IOOD+Dl -0.33RD-0I 

0.314D+00 

0.241D+00 

-0,2680-01 

0.5420-01 

O.lflfiD+OO 

-0.33SD-01 0.1000+01 

0. 1630+00 

-0 .2490+00 

-0.1140+00 

-0. 5420+00 

■0 .2040+00 

0.314D+C0 0. 1630 + 00 

O.lCCD+01 

-0 .6180-01 

0,2630-01 

-0.659D-0? 

0.227D+00 

0.241D+00 -0.2490+00 

-0.618D-01 

O.lOOD+01 

0,2^30-07 

0.3370-01 

•0.2360+00 

^0.26PD-01 -0.114D+00 

0. 2630-01 

0 .2330-01 

O.looO+oi 

0,3590-01 

•0.4400+00 

0,5420-01 -0.5420 + 00 

-0.6590-02 

0. 3370-01 

0. "590-01 

0 , 1 00 D +0 1 
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Table 21 

Transformation 

ED-*>0 -TO- W.I.36 ( 3 P4RAMFETER) 

♦****i»#*******<‘X'=>=^<‘=F 


original page lb 


SPLUTinri FPR 3 TRANSLATION PARAMFTFRS 
(UNITS - METERS) 


(USING VARIANCES ONLY) 


DU 

METERS 

-0.9A066S310+02 
± 2.8A- 


DV 

METERS 

-0. 1191A556D+03 
+ 2.95 


DW 

METERS 


“0.126AE272P+03 

♦ 2.TP 


VARIANCE “ COVARIANCE MATRIX 


MO 2= 1.07 

0.804138910+01 

0.0 

0,0 


0,0 

0.e711£834D+01 

0.0 


0,0 

0,0 

0.773589100+01 


COEFFICIENTS OF CORRELATION 


0 . 


lCOOOOCOD+01 


0,0 


0.0 


0,0 


0. lOOOOCOOD+01 


0.0 


0,0 


0,0 


O.lOOOOOOOD+01 
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Table 21 (Continued) 


REPRODUCIBILITY OP T-r' 
ORIGINAL PAGE IS POOR 


Ru«;ir'iiAL? V 



Vl( FD-50 

1 


V2( W. 

1.36 ) 


VI 

- V2 


6006 

-17. F 

31.3 

6.5 

6006 

0,5 

”1 .5 

-0.0 

“18.3 

32.8 

A. 6 

6004 

-3.4 

2.0 

1.5 

8 004 

0,6 

-0.3 

“O, 1 

-4.0 

2.3 

2.6 

PC09 

7,0 

-9.4 

5.0 

8009 

”0,5 

0,9 

“0. 1 

7.5 

-10.2 

3.2 

8010 

0.9 

“14.7 


8C10 

40,0 

1.3 

“0.1 

0,9 

“16.1 

4.0 

8 011 

12.4 

-2.8 

“3.3 

eon 

“C.6 

O.i 

0. 1 

13. C 

“3.0 

“3.4 

8016 

2.3 

2.6 

“3.5 

8 016 

“0.2 

“0.3 

C.2 

2.5 

2.9 

”3.7 

8019 

-0.0 

-18.5 

5.8 

8019 

0.1 

5.8 

-0,4 

-1 ,0 

“24.3 

6,2 

8C20 

4, 6 

“14.8 

4.1 

8 030 

“0.6 

6.3 

“0.4 

5.2 

-21.1 

4.4 

8031 

6.9 

17.6 

“6.4 

8 C3 1 

“4.3 

“5.1 

l.C 

11.2 

22.8 

“P.3 

8032 

-9.0 

11.2 

14.6 

8032 

0.2 

-0.3 

- 0.2 

“9,1 

11.6 

14.9 

8033 

1.6 

0,2 

“2.3 

FC33 

“0. 2 

“C.O 

0,2 

1.8 

C.2 

“2.6 

8701 

2.8 

8.4 

“5.5 

8701 

-0.3 

-0.3 

0.4 

3.0 

8,6 

-'^.P 

8702 

“1.7 

-Ifi.n 

1.7 

870? 

0.3 

7.4 

“0,2 

-1 ,9 

-25 ,4 

1.9 

8703 

—4. 6 

4.2 

2.C 

8703 

1.5 

-1.1 

“0.3 

-6.1 

5.3 

7.3 

8705 

0.1 

2.1 

“0.5 

8 705 

“0.3 

-16.1 

1.2 

0,4 

18.2 

“1.7 

8706 

0.6 

“9,2 

0,6 

8706 

“0.0 

O.S 

-0.0 

0.6 

“ 0.0 

0. 6 

8711 

“10. i 

2C.0 

2.7 

6711 

0. 1 

“0.7 

“O.c 

-10.2 

20.6 

2.8 

8712 

-0.5 

5.6 

“2.5 

8712 

1.1 

“36.2 

6.6 

-1.7 

41 .8. 

-9,1 

8713 

-4,3 

-1.5 

3.7 

8713 

11,4 

16.4 

-19.4 

-15.7 

“17.9 

23.1 

8714 

-8.1 

11.6 

3.7 

8714 

31.9 

“71.6 

-27.4 

“30.0 

83.1 

51,1 

8716 

?7.8 

“19.9 

“30.7 

8 716 

“12.2 

7.9 

19.8 

40,0 

-27.P 

“50.6 

8718 

5.3 

7.4 

“6,8 

8718 

“1,3 

“2.1 

1.2 

6,6 

9.5 

-8,0 

8720 

“2.7 

0.1 

1.5 

8720 

2.G 

-0.1 

“0. 3 

-4. 8 

0.2 

1.8 

8722 

-3.3 

5.8 

3.9 

8722 

2.4 

-13.8 

“0. 6 

-5.6 

19.6 

4.5 

9004 

-6.0 

-21.1 

2.4 

9 004 

0, 3 

3.3 

“0,2 

-6.3 

“24.5 

2,6 


UNIT OF RESIDUALS (METERS) 
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Table 22 
Transformation 

en-50 -TO- W.K?6 ( VOL ) 

SOLUTION FOR ?. TRANSLATION, 1 SCALF AND ■; ROTATTn^? PA’?^M^T'=RS 

(USING VARIANCES ONLY) 

DU DV PW DELTA OMEGA OSI FPSTLON 

METERS METERS METERS (XI, 0+6) SECOVOS SECOfOS SECONDS 

-9A.10 -125.53 -175,?? €.1? C,70 -0,15 0.17 

♦. 3c45 1 ?.?fr 1 3.AB ± 2.77 + 0.77 ♦ ± 0.E6 

VARIANCE - COVARIANCE MATRIX 

S02= 1.0? 

0.119D-»f>2 -n,17«5n + on 0.17ir.+oi -0,225D-O5 0,6270-05 -0.3000-05 

-C.179D+0C 0.1060+02 -0.3PeO + 00 0,3010-05 0.161O-O5- 0.1720-^6 

0.1710+01 -O.BEEO+nc 0.1210+0? 0.1820-05 0.?56D-05 -0.25AD-05 

-0.225D-05 0.3010-05 o,lf-?n-n5 0.766D-11 -0.215D-1? O.?320-l? 

0.627D-05 C. 1610-05 0.3560-05 -0.215D-1? 0,1400-10 -0.71?0-I2 

-0.30Cr— C5 0,1720-06 -0.2540-05 0.232D-12 -''.7130-12 0.9E4D-11 

-0.3470-05 0.227D-O5 -O.P 170-05 -0.3760-12 -0.7690-11 0.1030-11 

COEFFICIENTS OF CORRELATION 

O.lOOO+Oi -0. 1590-01 0.142D+00 -0.2360+00 0.4860+00 -0,2780+00 -0,2400+00 

-0.159D-01 0.1000+01 -0.342D-01 0.334D+00 0,1320+00 0,1680-01 0.1660+on 

0,1420 + 00 -0,3420-01 O.ICOO+CI C.18'’D+C0 0.273D + 0C -0.233D+OC -0.561D+00 

-0.2360 + 00 0,3340 + 00 0.189D+00 O.IOOD+Ol -0.2080-01 0,2670-0] -0.3250-01 

0.4860+00 0.1320+00 0. 2730+00 -C.2080-01 0.1000+01 -0.607D-01 -0. 4^00+00 

-0.27PD+00 0.1680-01 -0.2330+00 0.2670-01 -0.6070-01 O.lOOD+01 0.7870-01 

-0.240D+00 0,1660+00 -0,5610 + 00 -0.3250-01 -0,4900 + 00 0,7970-01 0,lO0D+''l 

NOTE : THE POSITIVE AXES ARC TOWARDS KEAN 

GREENWICH, EAST, ANO CIO. THE SQTATIONS 

PRINTED ARE A LOUT ?.P.O, 2ND, ANO 1ST AXES RESPECTIVELY. 


-0,347r-Of 
0.227D-05- 
-0.817P-05 
-0.3760-1? 
-0.7 69 0-1 ] 
0.1030-11 
0.1 75 0-10 
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Table 22 (Continued) 



RFSIPUALS V * 



VH 1 

ED -50 

) 


V 2 ( K . 

1.36 ) 


VI 

- V 2 


6006 

- 7.8 

22.6 

“ 1.6 

6006 

0,2 

- 1.1 

0.0 

- P.O 

23.7 

- 1.6 

etioA 

- 3.0 

- 0.2 

0,2 

8004 

0.5 

0.0 

— O . 0 

- 3,5 

-0 . 2 

o . 3 

8009 

8.2 

“ 0.0 

1.7 

8 004 

- 0.6 

o . E 

- C.l 

8.7 


1.8 

8010 

“ 0,9 

- 14.4 

5.2 

8010 

0,0 

1.3 

- 0.2 

- 0.9 

- 15.6 

•*.4 

8011 

15.0 

0.1 

- 5.1 

8011 

- O.F 

- 0,0 

0.2 

15.8 

0,1 

- 5,2 

8016 

1.3 

2.4 

“ 3.0 

8 Cl 6 

- 0 . 1 

- 0.3 

0.2 

1 .4 

2.7 

- 3.2 

8019 

- 4.0 

- 17,5 

8.7 

8019 

0.2 

5.5 

- 0.6 

“ 4.2 

- 23.0 

9,3 

8 C 30 

4,8 

- 13.2 

4,1 

8030 

- 0.6 

5.6 

- 0.4 

5,4 

~ 1=.7 

4.5 

8031 

9.4 

1 * 5.3 

“ 4.2 

Ff-^l 

“ 6,1 

- 5.6 

2.6 

16.0 

24.4 

- 11 . P 

8032 

“ 11.0 

9.7 

15.7 

8032 

0.2 

“ 0.3 

“ 0.2 

“ 11.2 

10,0 

16.0 

8033 

1.2 

“ 0.7 

“ 2.6 

8 03 3 

- 0 . 7 

0.2 

0,3 

1 .4 

-O.p 

- 7.8 

8701 

- 0.0 

4.8 

“?.4 

8 701 

0 . 

- 0.1 

0.3 

“ 0.0 

4.0 

- 4.2 

8702 

- 1.1 

“ 17.5 

1.0 

8702 

0.2 

7.2 

" 0 . 1 

“1 .? 

- 24.6 

1 . 1 

8703 

- 3.1 

0.9 

“ 0,4 

8 703 

J.C 

- C .2 

0.1 

- 4,1 

1 .2 

- 0.5 

8705 

- 0.2 

2,6 

O.l 

8 705 

0.4 

“ 20.5 

“ 0,2 

“ 0.5 

23. 2 

0,3 

8 706 

- 2.1 

“ 6.9 

3.4 

8706 

0 . 1 

0.6 

“ 0.1 

- 2.2 

- 7.5 

3 . 5 

8711 

- 17.4 

1 8.4 

9.8 

8711 

0,2 

- 0.6 

“ 0.2 

- 17,6 

l°.o 

9,9 

8712 

- 1.6 

5.3 

" 1.8 

8 712 

3 . 5 

- 34.0 

4 . 4 

- 5.2 

3 °. 3 

- 6.7 

8713 

- 6,0 

“ 1.8 

4.6 

8 713 

16 , C 

10.7 

- 24 . 1 

~??,0 

-21 .5 

2 ^’. 6 

8714 

- 9.2 

1 l.F 

4.4 

8 714 

36.3 

- 72.8 

- 32.4 

-4 5,5 

84.6 

36 . S 

8716 

26.6 

- 14.2 

- 28,2 

871 6 

-11.7 

5.6 

18.2 

38.3 

- 10.8 

- 46.4 

8713 

4.5 

6.9 

- 6.5 

8718 

-1.1 

“ 2.0 

1.2 

5.6 

5 ,9 

- 7,6 

8 720 

“ 0,7 

“ 4.5 

- 2.1 

8 720 

0.5 

2.3 

0.4 

- 1.2 

- 6.9 

- 7.4 

8722 

4,6 

6.9 

“ 2 . P . 

8 722 

- 3.3 

- 16.4 

0,5 

8.0 

23.3 

- 3.3 

9004 

- 8.0 

“ 11.9 

8.2 

9004 

0.4 

1.9 

- 0.7 

- 8 -, 5 

- 13,8 

8.0 


UNIT OF RESIDUALS ( METERS ) 


♦Residuals in tiie Cartesian coordinate system. 
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Table 22 (Continued) 




RFSIPUALS V* 
SPHERICAL 



VI ( 

EH- 50 J 



V2 ( W. 

1 . 36 } 


VI - V7 


6006 

-0.54 

23.90 

-1.5 

6006 

0.12 

-1.08 

-0.0 

-0.66 

24.9? 


8004 

2.53 

0.3P 

-1 .6 

8004 

-0.41 

-0,06 

0,3 

2.94 

0,44 

- 1.9 

8009 

-4.. SO 

-0.60 

5.9 

80OQ 

0.35 

0.59 

-0.4 

-5.15 

-lO.AO 

6.3 

8010 

5.59 

-14.11 

1.9 

8010 

-0.26 

1.26 

0.0 

5.85 

-!5,3F 

1.9 

8011 

-14.97 

0.58 

5.2 

8011 

0 .70 

-0.03 

"0.3 

-15.67 

0,61 

5. 5 

8016 

-3.22 

2.24 

-1.2 

8016 

0 .20 

“0.25 

0.0 

-3.4? 

2.40 

-1.3 

8019 

10.58 

-16.83 

1.6 

8019 

-1 .07 

5.43 

0.3 

11.64 

-22.26 

1.3 

8030 

-0.48 

-13.35 

5.9 

8030 

0 .03 

5.58 

-0.5 

-0 . 51 

-18.94 

6.5 

8031 

-12.44 

19.82 

-2.7 

8031 

6 .27 

-5.94 

-1.1 

-18.70 

25.76 

-1. 5 

8032 

17.16 

1 1.62 

5.7 

8032 

-0.26 

-0.35 

-0.1 

17.42 

11.95 

5.7 

8033 

-2.45 

-0.84 

-1.3 

60 33 

0.29 

0.19 

0,1 

-2.74 

-1,0? 

-1 . 4 

8701 

-3.58 

4.62 

-2.0 

8.70 i 

0.23 

-C.14 

0.2 

-3.81 

4.76 

-2. 2 

8702 

2.52 

-17.35 

“O.R 

B70? 

-0,63 

7.12 

0.4 

3.14 

-24,47 

-1. 1 

8 703 

2.10 

1.57 

-2.0 

8703 

“0,72 

“0.46 

0.6 

2.8^ 

2.0? 

-2.5 

8705 

0.20 

2.62 

“0.1 

8705 

-0 .5A 

-20.54 

0.2 

0.74 

23.16 

-0.3 

8706 

4.34 

-6.72 

0.4 

8706 

-0,16 

0.56 

0.0 

4.51 

-7.2B 

0.4 

8711 

15.03 

2 2.2° 

-3.6 

8711 

-0.17 

-0.64 

-0.1 

15,19 

22.93 

-3. 5 

8712 

-1.06 

5.52 

-1.5 

8712 

6.7? 

-3 3, 89 

0.2 

-7.79 

39,41 

-1.7 

8713 

7.55 

0.11 

-1.8 

8713 

-32 .09 

13.94 

0.5 

39.64 

-13.84 

-2.3 

8714 

7.97 

13.10 

-2.8 

B714 

-40,50 

-77,67 

“1.7 

48.47 

90.78 

-1. 1 

8716 

-39.25 

-12.39 

2.6 

£716 

21 ,64 

4.86 

2.7 

-60 . 90 

-17,28 

-0.0 

8718 

-8.40 

6.16 

-1.3 

8718 

1 .78 

-1.79 

-0.0 

-10.18 

7.95 

-1.3 

8720 

0.69 

-4.12 

-2,8 

8720 

-0..S3 

2.08 

0.0 

1.5? 

—6.20 

-3.7 

8722 

-2,67 

8.13 

-1.8 

8.722 

- 2.69 

-16,47 

1.8 

0.0? 

24.60 

-3.7 

9004 

10.56 

-12.73 

-0.5 

9004 

-0 .71 

1.93 

-0.2 

11.28 

-14.66 

-0.7 


UMIT OF RESIDUALS IMETERS) 


♦ The spherical residuals represent the residuals in the Cartesian coordinate system as 
transformed into the curvilinear (<p, X , h ) system. 
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Table 23 

T ransf ormation 


j 


JG 


ED-50 -TO- W.I.36 ( VFT5 MOPEL ! 

^i^i^i^tf************ ******* **'^*** *^^^**** 

SOLUTION FOR 3 TRANSLATIONt . 1 SCALE AND ? ROTATION PAR^METPRS 


<USING VARIANCES ONLY) 


DU 

DV 

DW 

DE LTA 

ALPHA 

KSI 

ETA 

METERS 

METERS 

METERS 

(Xl.D+6) 

SECONDS 

SECONDS 

SECONDS 


-115.53 “125.32 

4.12 

0.63 

“O.IP 

-0.32 

± 3 .45 

♦ 3.24 

± 3.48 

± 2.77 

i 0.5S 

+ 0.44 

♦ l.O'' 


VARIANCE - COVARIANCE MATRIX 


S02= 1.03 


0.119D+02 

-0.17QD+0C 0.171D+01 

-0.225D-05 

0. 2490-05 

-Q.214D-05 

-0.7C'4D-0 5 

-0.179D+00 

0.106D+02 -0.3P8D+00 

0.301D-05 

0.2650-05 

- 0 . 3470-06 

0 .?.oin -06 

0.171D+01 

“0.388D+C0 0.1210+02 

0. 182D-05 

-0.239D-05 

-0 ,630D-^6 

~0 ,pO3D-0‘5 

-0.2 25D-0 5 

0.3C1P-05 0.1820-05 

0.766D-11 

-0.3620-12 

0.311D-12 

- 0 . 1170-12 

0.249D-05 

0.265D-05 -0.2390-05 

-0.362D-12 

0.796D-11 

0.358D-12 

-0.2270-1? 

-0.214D-05 

-0.3470-06 -0.6300-06 

0.311D-12 

0,3500-12 

0.O7F0-11 

-0.12^0-11 

“0.7040-05 

0.8010-06 -0.8950-05 

-0.1170-1? 

-0.2270-1? 

-0.1240-11 

0,237 0-10 


COEFFICIENTS OF CORRELATION 



O.lOOD+01 

-0.159D-01 0.1^20+00 

“0.236D+00 

0,2560+00 

-0.19OH+00 

-0.420D+00 

- 0 . 1590-01 

0.1000 + 01 -0,3420-01 

0,3340+00 

0,2880+00 

- 0 . 3400-01 

0.5 040-01 

0.1420+00 

-0.542D-01 O.lOOD+01 

0.189D+00 

-n, 2440 + 00 

-0.57OD-O1 

-0 .528D+00 

“0.236D+00 

0.334D+CC C.169D+C0 

O.lOOD+01 

-0 .4640-01 

0.360D-01 

-0.8730-0? 

0 . 2560+00 

0.280D+00 -0.2440+00 

-0.4640-01 

0.1 000 +01 

0.406P-01 

-0.1650-01 

-0.199D+00 

-0.340D-01 -0.579D-01 

('.3600-01 

0,4060-01 

0 .lOOP+Ol 

- 0 . 8190-01 

-0.420D+C0 

0.5040-01 -0. 5200+00 

-0.E73D-02 

-0.1650-01 

-0,8180-01 

O.io-JD+Ol 
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Table 24 
T ransf ormat ion 


)ly ! 


orige;^^^ 




ED-50 -TO- W.I*36 ( BURSA } 

41*41 * 411 ^ 4i4i4M^4»ti 4c 4i4t4!4i4r4!4t4< 

SOLUTION FOR 3 TRANSLATION, I SCALE AND 3 ROTATION PARAMETERS 


fUSING VARIANCES ONLY) 


DU OV DW DELTA OMEGA PS I EPSILON 

METERS METERS METERS (Xl»D+6) SECONDS SECONDS SECONDS 

-123.91 -112.22 -152.66 6.12 0.70 -0.15 0.17 

118.62 130.64 1 I 8 . 8 O 1 2.77 1 0.77 + 0.65 + 0,86 

VARIANCE - COVARIANCE MATRIX 


S02= 1.03 


0;354D+03 

-0.695D+02 

-0.372D + 02 -0.300D-04 

-0.8880-05 

0,4620-04 

0.994D-05 

-0.695D+02 

0 .939D+03 

-0.503D + 01 -0.267D-05 

0.9380-04 

-0.7940-05 

-O.1 15 0-03. 

-0.372D+02 

-0 .5030+01 

0.353D+03 -0,3790-04 

0. 5680-06 

-0.402D-04 

0.526D-05 

-0.300D-04 

-0.267D-05 

-0.379D-04 0,7660-11 

-0.2150-12 

0,2320-12 

-0.376D-I 2 

-0,8 88D~05 

0.938D-04 

0.5680-06 -0.2150-12 

0.1400-10 

-0.713D-12 

-0. 7690-11 

0.462D-04 

-0.794D-05 

-0.4020-04 0.232D-12 

-0.713D-12 

0.984D-11 

0.103D-11 

0,9940-05 

-0.115D-03 

0.526D-05 -0,3760-12 

-0.7690-11 

0.103D-11 

0. 1750-10 



COEFFICIENTS OF CORRELATION 



O.lOOD+01 

-0.121D+00 

-0.105D+00 -0.576D+00 

-0.1260+00 

0.783D+00 

0.126D+00 

-0.121D+00 

O.lOOD+01 

-0.873D-02 -0. 3150-01 

0. 8170+00 

-0.827D-01 

-0. 8950+00 

-0.105D+00 

-0 .8730-02 

O.IOOD+Ol -0.7290+00 

0.8060-02 

-0.682D+00 

0.668D-01 

-0.576D+00 

-0.315D-01 

-0. 7290+00 0. 1000+01 

-0.2080-01 

0.267D-01 

-0 .3250-01 

-0.126D+00 

0.817D + 00 

0. 8060-02 -0. 2080-01 

0,1000+01 

-0.607D-01 

-0 .490D+00 

0.783D400 

-0.827D-01 

-0.682D+00 0.267D-01 

-0. 6070-01 

O.lOOD+01 

0.707D-01 

0.126D+00 

-0.895D+0C 

0.668D-01 -0.325D-01 

-0.4900+00 

0.787D-01 

O.IOOD+Ol 


NOTE : THE POSITIVE AXES ARE TOWARDS MEAN 


GREENWICH, EAST, AND CIO. THE ROTATIONS 

PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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REPRODUC.IBiLI'r/ r;, ’ 
Table 25 ORIGINAL PAG^^ 

Transformation 

NGS -TO- W.1,36 ( Mnt MOtiFL ) 

SOLUTION FOR 3 TRANSLMinN, 1 SCALE AND 3 ROTATION PARAMETERS 

I USING VARIANCES ONLY) 


DU 

nv 

DW 

DELTA 

OMEGA 

PSI 

EPSILON 

METERS 

METERS 

METERS 

(Xl.D+6) 

SECONDS 

SECONDS 

SECONDS 

-9.17 

-10.04 

I5.9S . 

0.90 

“1 .68 

0.55 

2,10 

± 2.68 

± 3.3C 

♦ 3.4P 

± 1.62 

+ 1.58 

+ 0,53 

± 2.13 



VARIANCE - COVARIANCE MATRIX 



S02= 

“ 1.07 





0.717D+01 

0 .259D+01 

0.2P3D+-01 

0.100D-05 -0. 9780-05 

0,?54n-05 

0.1270-04 

0.259D+01 

0 ,109D+02 

0.477D+-C1 

-0.401D-06 -0.1640-04 

0,4140-05 

0,2 100-0^* 

O.283D+01 

0.477D + 01 

0. 1220 + 02 

-C.267D-06 -0.170D-04 

0.5510-05 

0.2370-04 

O.lOOD-05 

-0.401D-06 

-0.267D-06 

0.263D-11 -0.190D-1?. 

-0. 3260-12 

-0.7 97 0-1? 

-0.9 780-05 

-0 .164D-C4 

-0.17CD-04 

-0,190[)-13 0.585D-I0 

--n,14.6D-in 

-0.756D-in 

0.2 54D-05 

0 .414D-05 

0.551D-C5 

-0.326D-1? -0.I46D-10 

0,6690-1 1 

0. 1900-10 

0.127D-04 

0.21CD-C4 

0.2370-04 

-0.797D-13 -0.7560-10 

0. 199D-10 

0.1C6D-09 



CDEFFIC lENTS OF CORRELATION 



O.lOOD+01 

0.292D+CC 

0.3G2O+00 

0.230D+00 “0.478D+00 

0.3670+00 

0.459D+0C 

O.292D+O0 

O.lOOD + 01 

0.413D+00 

-0.7480-01 -0.6490+00 

0.4840+00 

0 ,6170+00 

0.302D400 

0 .413D+00 

O.lOOD+01 

-0.4720-01 “0.6360+00 

0.6100+00 

0.6 57 0+0 6 

0.230D+00 

-0 .7480-01 

-0.472P-01 

0,1000+01 -0.153D-02 

-0,7780-01 

-0.476 0-0 2 

-0.478D+00 

-0 .6490 + 00 

—0 .6360+ no 

-0,1530-02 O.IOOP+OI 

-0.740D+P0 

-0 .9590+00 

C.367D+00 

0 ,484D + 0(. 

0,6100+00 

-0.7780-01 -0.7400+00 

O.lOOD+01 

0.7450+06 

0.459D+OD 

0 ,6170+00 

0,6570+00 

-0,4760-02 -0.9590+00 

0.7450 +00 

0. 1000+01 

NOTE : the 

positive axes are towards mean 



GREENWICH, EAST 

, AND cin. 

THE ROTATIONS 




PRINTED ARE ABOUT 3RDt 2NDt AND 1ST AXES RESPECTIVELY. 
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Table 25 (Continued) 


REPRODUClGnjri 
ORIGINAL PAGE A 


RESIDUALS V 


Vl( NGS I V?( W.I.36 ) VI - '/2 


6006 

“0.5 

1.1 

“0.0 

6006 

0.2 

“0.5 

0.0 

“0.7 

. l^T 

“1.0 

6016 

“1.9 

l.A 

“0.4 

6016 

0.? 

“0.6 

0.1 

“2.2 

1.9 

“0. 4 

6065 

2. A 

“2.3 

1.2 

6065 

“0, 5 

0.9 

“0.2 

3.0 

“3.1 

1.3 


UNIT OF RFSIDUALS (METERS) 
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Table 26 

REPRODUCIBILITY C ” 
Transformation ORIGINAL PAGE L3 Lv> .• 


NGS -TO- W.r.36 t BURSA ) 

SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PARAMETERS 


fUSING VARIANCES ONLY) 


DU DV DW DELTA OMEGA PSI EPSILON 

METERS METERS METERS (Xl.D-t-6) SECONDS SECONDS SECONDS 

6.19 “92.95 8.96 0.90 -1.68 0.55 2.10 

♦19.97 +77.89 +10.30 ♦ 1.62 +1.58 + 0.53 ♦ 2.13 


- 

VARIANCE 

; - COVARIANCE MATRIX 



S02s 1.07 






0.399D+03 “0.128D+0A 

0.370D+02 

“0.116D-04 

“0.127D-03 

0.474D-04 

0.169D-03 

0.128D+04 0.607D+04 

‘'0.1B7D+03 

-0.226D-05 

0.586D-03 

-0,1510-03 

-0.798D-03 

0.370D+02 “0.187D+03 

0,1060+03 

-0.U3DH54 

-0.173D-04 

-0 .4840-05 

0.2770-04 

0.116D-04 -0.226D-05 

-0.113D-OA 

0.263D-11 

“0.190D-13 

-0.326D-12 

-0.797D-13 

O.l 27D-03 0 .586D-03 

-0.173D-04 

“0. 190D-13 

0.585D-10 

-0.146D-10 

-0.756D-10 

0.A74D-04 -0.151D-03 

“0.464D-05 

-0.326D-12 

“0,1460-10 

0.669D-11 

0.1990-10 

0.169D-03 “0 .7980-03 

0.277D-04 

“0.797D“13 

-0.756D-10 

0.1990-10 

0.106D-09 


COEFFICIENTS OF CORRELATION 


0.1000+01 -0. 8230+00 

0.180D+00 -0,3580+00 

-0. 8290+00 

0,9180+00 

0.818 0+00 

0.8230+00 

O.lOOD+01 

-0,2330+00 -0.1790-01 

0. 9830+00 

-0.7480+00 

-0.9930+00 

0.1800+00 -0.233D+00 

0,1000+01 -0.6780+00 

-0.2200+00 

-0,1820+00 

0. 2610+00 

0.3 580+00 -0.179D-01 

-0.6780+00 0.1000+01 

-0.1530-02 

-0.7780-01 

-0 .4760-02 

0.8290+00 

0.983D + 00 

-0.2200+00 -0.1530-02 

0.1000+01 

-0.7400+00 

-0.9590+00 

0.918D+00 -0.748D+00 

-0.1820+00 -0.7780-01 

-0.7 400+00 

0. 1000+01 

0.7450+00 

0,8180+00 -0,9930+00 

0.2610+00 -0,4760-02 

-0. 9590+00 

0.7450+00 

0. 1000+01 


NOTE I THE POSITIVE AXES ARE TOWARDS MEAN 

GREENWICH, EAST, AND CIO. THE ROTATIONS 


PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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Table 27 
Transformation 


REPR0DDX3IBILKY - ' 
ORIGINAL PAGE IS : 


SAD III “Tn~ W.1.3f ( MOL MGDfL ) 

SOLUTION FOP 3 TPANSLAIIONt 1 SCALE AND 3 ROTATION PARAMETERS 


(USING variances ONLY) 


DU DV DW DELTA OMEGA PS I EPSILON 

METERS METERS METERS (Xl.D+6) SECONDS SECONDS SECONDS 

“8.17 -16.30 R.66 -0.32 “0o23 -0^38 0,53 

+ 2,53 * 2.88 + 2.60 + 2.45 + 0.68 ♦ 0.67 + 0,79 




VARIANCE 

- COVARIANCE MAT»IX 



S02= 

: 0.P9 






0.6A0D+01 

“0 .236D + C0 

0.629D+C0 

-0.1820-05 

0.393D-05 

-o.i?in-o5 

-0.1680-05 

•0.236D+00 

0.E-2RD + 01 

-O.lRiPO+OO 

0.3190-05 

0.1 520-05 

-0, 3050-06 

0.2720-05 

0.629D+00 

-0.199D+00 

0.6740+01 

0. 1750-05 

0.1750-05 

-0.9320-06 

-0,5050-05 

0.182D-D5 

0.319D-C5 

0.1750-05 

0. 5980-11 

0.5750-12 

-0.1360-12 

0.4630-12 

0.3 93D-05 

0.152D-05 

0.1750-05 

0,5750-12 

O.llOD-lO 

0,4800-11 

-0.643 0-11 

C.121D-05 

-0,3050-06 

-0.9320-06 

-0.1360-12 

0.4(500-11 

0.I07D-10 

-0.5370-11 

0.16RD-05 

0.272D-05 

-0. 5050-05 

0.A630-12 

-0. 6430-11 

-0.5370-11 

0. 143 0-10 



COEFFICIENTS OF CORRELATION 



O.lfOD+01 

-C. 3250-01 

0.9500-01 

-0.2940 +00 

0 ,4690+00 

-0. 1460+00 

-0,1760+00 

•0.325D-01 

o.iooD-t-ni 

-0.2660-01 

0.4530+00 

0.1600+00 

-0,3250-01 

0,2 500+00 

0.9 59D-01 

-0.2660-01 

0.1000+01 

0.2760+00 

0.2030+00 

“0. 1100+00 

-0.5150+00 

■0.29AD400 

0,4530+00 

0,2760+00 

O.lOOP+Ol 

0.7100-01 

-0,1700-01 

0.5010-0! 

0.A69D+00 

0 ,1600 + 00 

0.2030+00 

0.710D-01 

0 .1000+01 

0.4440+00 

-0.5 14 0+00 

•0.146D+0C 

-0.3250-01 

-0,1100+00 

-0.17CD-01 

0 ,4440+00 

0.1000+01 

-0.4350+00 

0.176D+00 

0.2500+00 

-0.5150+00 

0. 5010-01 

-0.5140+00 

-0.4350+00 

0. 100 0+01 


NOTE ; THE POSITIVE AXES ARE TOVARDS MEAN 


GREENWICH, EAST, AND CIO. THE ROTATIONS 


PRINTED ARE AnOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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Table 27 (Continued) 


VII SAO III) 


6006 

“13.6 

1.6 

-7.5 

6016 

-14.. 0 

P.0 

“4.0 

8010 

l.I 

-1.8 

?.l 

8011 

6,7 

1.7 

-7.1 

8015 

1.2 

0.6 

-0.2 

8019 

-5.0 

-10.7 

-3.3 

9004 

-1.5 

-0.1 

O.A 


RESIDUALS V 


V2( W.I.36 ) 


6006 

C.5 

-0. 1 

0.1 

6016 

0.3 

“0.5 

0.2 

8010 

-C.6 

2.4. 

-1.0 

8011 

-C.4. 

-C.l 

0.3 

8015 

-1.0 

-2.4 

0.3 

8019 

0.6 

6.4 

0.4 

9004 

1.9 

0.2 

-0.7 


REPRODUCIBILOT 0’ 
original page is l 


VI 

- V2 


-14.1 

1.7 

-7.5 

-14.2 


-5.0 

1.7 

“4,2 

3.1 

7.1 

1.8 

-7,4 

3.0 

3.1 

-0.5 

-5.6 

-17.1 

-8.7 

-3.4 

-0.2 

1.0 


UNIT OF RESIDUALS IMETFRS) 
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DU 

DV 

OW 

DELTA 

OMEGA 

psr 

epsilon 

METERS 

METERS 

METERS 

(Xl.D+6) 

SECONDS 

SECONDS 

SECONDS 

-14.66 

-31.67 

21.12 

-0.32 

-0.23 

-0.38 

0.53 

♦18 .04 

±27.56 

±18.68 

± 2.45 

± 0.68 

± 0.67 

± 0.78 


VARIANCE - COVARIANCE MATRIX 



S02*= 0.99 


0.325D+03 

0.181D+03 

-0.914D+02 

“0.281D-04 

0.1 510-04 

0.4550-04 

—0.2360—04 

0.1810+03 

0.759D+03 

“0.230D+03 

-0.148D-O5 

0.7760-04 

0.4530-04 

-0,9190-04 

0.9140+02 

-0.230D+03 

0.349D+03 

-0.254D-04 

-0.2650-04 

-0.4960-04 

0,2710-04 

0.281D-04 

-0.148D-05 

-0.254D-04 

0.598D-U 

0. 5750-12 

-0.1360-12 

0. 4630-12 1 

0.151D-04 

0. 7760 -04 

-0.265D-04 

0. 5750-12 

0.1100-10 

0.4800-11 

-0.6430-11 

0.455D-04 

0.453D-04 

-0.496D-04 

-0.135D-12 

0.4800-11 

0.1070-10 

-0.5370-11 

0.236D-04 

-0.919D-04 

0.271D-04 

0.463D-12 

-0.6430-11 

-0.5370-n 

0. 1430-10 



COEFFICIENTS OF CORRELATION 



O.lOOD+01 

0 .3650+00 

-0.2710+00 

-0,6380+00 

0.2 520+00 

0.7730+00 

-0. 3460+00 

0.365D+00 

O.lOOD+01 

-0.446D+00 

-0.219D-01 

0.8490+00 

0.5030+00 

-0.8830+00 

0.271D+00 

-0 .446D + 00 

O.lOOD+01 

-0.556D+00 

-0 .4280+00 

-0. 8130+00 

0.3840+00 

•0.63CD+00 

-0.219D-01 

-0.556D+00 

O.lOOD+01 

0.7100-01 

-0.1700-01 

0. 5010-01 

0.2 52D+00 

0.849D+00 

-0.428D+00 

0.710D-01 

0. 1000+01 

0.4440+00 

-0. 5140+00 

0.773D+00 

0 .503D+00 

-0.813D+00 

-0.170D-01 

0.4440+00 

0.1000+01 

-0.4350+00 

■0.346D+00 

-0.883D+00 

0.384D+00 

0.501D-01 

-0.5140+00 

-0.4350+00 

O.lOOD+01 


NOTE : THE POSITIVE AXES ARE TOWARDS MEAN 

GREFNWICHt EAST, AND CIO. THE ROTATIONS 

PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 



61 


Table 29 

Transformation 

WN-14 -TO- W.I.36 ( MOL MODEL ) 

SOLUTIOrj FOR 3 TRAMSLATI ON, 1 SCALE ANO 3 ROTATION PARAMETERS 

(USING VARIANCES ONLY) 

DU PV PW DELTA OMEGA PSI EPSILON 

METERS METERS METERS (Xl.D+6) SECONDS SECONDS SECONDS 



0.19 0, 

68 0.12 

-1.19 0.C4 -0 

.12 -Or- 

47 

± 1.73 i 2. 

41 1 1.93 

+1.42 ±0.50 ± 0 

.31 + 0,! 

S3 



VARIANCE 

- COVARIANCE MATRIX 



S02 

= 1.00 





0.301D+01 

-0.887D-OX 

0.107D-01 

0.9910-07 0.3970-06 

-0,3770-06 

-0.317 0-0 6 

-0.BP7P-01 

0 ,5800 + 01 

-0.156D+00 ■ 

“0. 1250-06 -0. 1210-05 

0.74B0-07 

0.131 D-05 

0.107D-01 

“0,156D+0G 

0.374D+01 

0.1390-06 0.5150-06 

0.1570-06 

“0,iri6D-n5 

0.991D-07 

“0.125D-06 

0.13PP-06 

0.2020-11 0,8590-1? 

-0. 2330-12 

0.7340-13 

0.397D-06 

“0.121P-C5 

0.515D-06 

C. 6590-13 0. 5790-11 

-0.281D-12 

-0.4590-11 

-0.377D-06 

0,7480-07 

0.157D-06 

-0. 2330-12 -0,2810-12 

0.2270-11 

0,2660-12 

-0.317D-06 

0.131D-05 

-0.1 060 -05 

0.7340-13 -0,4590-11 

0.2660-12 

0. *>340-11 



COEFFICIENTS OF CORRELATION 



O.lOOD+01 

-0.212D-01 

0. 3200-02 

0. 4020-01 0.9500-01 

-0,1440+00 

-0.5*>5D-01 

-0.212D-C1 

O.IOOD+Ol 

-0.3 3 50-01 

-a, 3660-01 -0.2100+00 

0.2060-01 

0 ,178 0+00 

0.320D-02 

-0.335D-01 

O.lOOD+01 

0. 5070-01 0. 1110 + 00 

0.5390-01 

-0.1790+00 

O.AOZP-Ol 

-0.366D-01 

0.5070-01 

O.lOOD+01 0.2510-01 

-0.1090+00 

0.1690-01 

0.95CD-01 

-0.21CD+CC 

o.iiio+ro 

0.2510-01 0.1000+01 

•«0. 7750-01 

-0.6240+00 

-0.144D+00 

0.206D-01 

0.53PD-01 

“0,1090+00 -0.7750-01 

0.1000+01 

0.5770-01 

-0.598D-01 

0.1780+00 

-0.1790+00 

0,1690-01 -0.6240+00 

0.5770-01 

0 .1000+01 

NOTE : THE 

POSITIVE AXES ARE TOWARDS MEAN 



GREENWICH, EAST, 

AND CIO. THE POTATIONS 



PRINTED ARE APOUT 3RD, 2ND, 

AND 1ST AXES RESPECTIVELY. 
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Table 29 (Continued) 



VI ( 

WN-14 

) 

AOOA 

-1.3 

-0.3 

1.9 

6016 

1.3 

-2.0 

-1.2 

6065 

-0.6 

-5.0 

1.7 

8009 

“G. 1 

6.1 

-3.3 

8010 

-8.C 

9.6 

6.2 

8011 

14. 

57.6 

-12.5 

8 015 

0.6 

22.3 

“5.0 

8019 

-0.2 

P.4 

-0.2 

8030 

1.6 

12.6 

-0.3 

9004 

-0.0 

1 8.9 

1.3 


RFSIDUALS V 



V2( W. 

1.36 ) 


6006 

0.4 

0.1 

-0.1 

6016 

-0.3 

o.e 

0.3 

606 5 

0.2 

1.5 

“0.3 

8000 

0.0 

-0.3 

0.2 

8010 

C.6 

-o.e 

“0.4 

8011 

-0.6 

-0.7 

0.5 

8015 

-0.1 

“1.7 

0.5 

9019 

0.0 

-2.5 

0.0 

8 030 

-o. 3 

“3.4 

0.1 

9004 

0.0 

“1.8 

“0.4 


UNIT OF RESIDUALS (METERS) 


reproducibility op > 
ORIGINAL PAGE IS POOr. 


VI 

- V2 


“1.7 

-0.5 

2.0 

1.5 

-2.fi 

-1.4 

“0.8 

-6.5 

2.0 

“0.1 

6.4 

-3.5 

-8.6 

10.4 

6,6 

15.5 

58.3 

-13.1 

0,6 

24.1 

-5,5 

-0.2 

10.9 

-0.2 

1.9 

16.1 

-0.4 

-0.0 

20.7 

l.C 


63 



Table 30 

Transformation 

WN-X4 -TO- W.1.36 I BURSA )t 

*«««*«« »«*****« 4 >«***«««*:^« ******** 

SOLUTION FOR 3 TRANSLATION, I SCALE AND 3 ROTATION PARAMETERS 

lUSING VARIANCES ONLY! 

DU DV OW DELTA OMEGA PSI EPSILON 

METERS METERS METERS (Xl.0+6) SECONDS SECONDS SECONDS 

2.41 12.98 6.24 -1.19 0.04 -0.12 -0.47 

♦10,00 ♦21.82 ♦8.89 ♦ 1.42 ♦ 0.50 ±0.31 +0.63 

VARIANCE - COVARIANCE MATRIX 


EEPROD0CIBIIJXY r- 

original 



S02* 1,00 


O.lOOD+03 -0.364D+02 

-0.704D-01 

-0.9590-05 

-0,6040-05 

0.1150-04 

0. 4390-05 

-0.364D+02 0.476D+03 

-0.266D+02 

-0. 1770-05 

0,4470-04 

-0.2170-05 

-0 .6200-04 

-0.704D-01 -0.266D+02 

0.791D+02 

-0,8330-05 

-0.247D-05 

-0,8100-05 

0.5140-05 

-0.959D-05 -0,1770-05 

-0.833D-05 

0.2020-U 

0.859D-13 

-0.2330-12 

0.7340-13 

-0.6040-05 0.447D-04 

-0.247D-05 

0.B59D-I3 

0.579D-U 

-0.2810-12 

-0,4590-11 

0.115D-04 -0.217D-05 

-0.810D-05 

-0.233D-12 

-0.281D-12 

0. 2270-11 

0.2660-12 

0.439D-05 -0.620D-04 

0.514D-05 

0.734D-13 

-Oc459D-ll 

0.2660-12 

0.9340-11 


COEFFICIENTS OF CORRELATION 



O.lOOD+01 -0.167D+00 

-0.7920-03 

“0.674D+00 

-0.251D+00 

0.7640+00 

0.1440+00 

-0.167D+00 O.IOOD+Ol 

-0.137D+00 

-0.570D-01 

0.8510+00 

-0.659D-01 

-0 .9300+00 

-0.792D-03 -0.137D+00 

O.lOOD+01 

-0.659D+00 

-0.115D+00 

-0.6040+00 

0.1890+00 

-0.674D+00 -0.570D-01 

-0.659D+00 

O.IOOD+Ol 

0,2510-01 

-0. 1090+00 

0. 1690-01 

-0.251D+00 0.851D+00 

-0,1150+00 

0.251D-01 

O.IOOD+Ol 

-0.7750-01 

-0.6240+00 

0.764D+00 -0.659D-01 

-0.604D+00 

-0.109D+00 

-0,7750-01 

0.1000+01 

0.5770-01 

0.144D+00 -0.930D+00 

0,1890+00 

0.169D-01 

-0,6240+00 

0.5770-01 

0.1000+01 


NOTE : THE POSITIVE AXES ARE TOWARDS MEAN 

GREENWICH, EAST, AND CIO. THE ROTATIONS 

PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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Table 31 

T ransf ormat ion 

GFM 6 -TO- W.I.36 { MOL 'inDFL ) 

SOLUTID.N for 3 TAAr;5LATin.‘J, 1 SCALE AND ? ROTATION PARAMETERS 

I USING VARIANCES ONLY) 


ru 

DV 

DW DELTA 

OMEGA 

PSI 

=P$IL=N 

METERS 

METERS 

METERS (X 1,0+6) 

S ECO NOS SEC ON OS 

SECOMOS 

-19.53 

-30.67 

6.00 1.56 

—0 .66 

-0.73 

-0,5? 

+ 1 .98 

i 2.25 

+2.00 +1.35 

+ 0.37 

± 0,31 

♦ 0.45 


VARIANCE - COVARIANCE MATRIX 


502= l.OA- 

0.390D+01 -C,901D-C1 0,2 = 00-01 “C.A5FD-06 0.73«>n-C6 -C.PAf»D~Of. 

-0.901D-f>l O.507D+01 -0. 1730 + 00 0. 18^-0-06 -0.312D-06 -0. 1310-06 0,l?4r— 

0,2900-01 -0.1730+00 0.39=0+01 0.6010-05 0.4090-06 -0.145D-06 -0, =010-06 

-0.458D-06 0.1540-06 0.6C10-C6 0,1520-11 0.2090-12 0.5050-13 -f),TA4P-lA 

0.739D-06 -0.3120-06 0,40=0-06 0,20=0-12 0. 3190-11 0.4730-12 -0,1=90-11 

-0«868n-06 -0.1310-06 -0.1450-06 0.5050-13 0,4730-1? 0.2230-11 -0.5=50-12 

-0.356D-06 0.124D-05 -O.901D-06 -0.744p-l4 -0.1980-11 -0.595r?-I2 0,4830-11 

COEFFICIENTS OF CORRELATION 

0,1000+01 -0.2030-01 0.735D-O2 -0.172D+00 0.209D+OO -0,2=40+00 -O,9?oo-oi 

-0,2030-01 O.lOOD+01 -0. 3850-01 C.6C6D-01 -0.7750-01 -0.3=10-01 o. =510+0= 

0.735D-02 -0.3850-01 O.lOOD+01 0,2230+00 0.1150+00 -0.4R5D-01 -0.2050+00 

-0.172D+00 0.6060-01 0.223D+00 O.IOOD+Ol 0.868D-01 0.251D-01 -0.251D-02 

0.20=0+00 -0.7750-01 0.115D+00 0,8680-01 O.lCOD+Oi 0.177D+00 -0.5040+00 

-0.2940+00 -0.3910-01 -0.4850-01 O. 2510-01 0.1770+00 0.1000+01 -O.lPlC+00 

-0.8200-01 0.2510+00 -0.2050+00 -0.251D-02 -0.5040+00 -0.1810+00 0. 1000+01 

NOTE : the positive AXES ARE TOWARDS MEAN 

GREFr^WICHt FAST, AND CIO. THE ROTATIONS 

PRINTED ARE AB=UT 3RD, 2ND, AND 1ST AXES RES PECT IVELY . 
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Table 31 (Continued) 


RESIDUAL?: V 


n ' 

i- ^ ^r\Av<! 




ml 


VI ( GEM 6 J 


V2{ W.I.36 > 


VI - 


V? 


6C06 

6016 

6065 

9004 

-o.e 

-1.6 

0,3 

0.2 

4,6 

-O.P 

-3.2 

0.1 

-0.7 
C .1 
2.6 
-0,1 

6006 
601 6 
606 5 
9004 

C.2 

0. A 
-0. 1 
-1,0 

— 1 ® 6 
G.5 
1 . 1 
-3.1 

O.C 

-0.0 

-0.4 

2.3 

-1.0 

-2.0 

0.4 

2.1 

6.2 
-1 .4 
-4,3 
3.2 

-0.7 

0.1 

3.1 

-2.4 





UNIT OF 

RESIDUALS 

(METERS ) 





66 


Table 32 -REPRODUCIBILITY OF T 

Transformation OEIGINAL PAGE IS POO.' 

GEH 6 -TO* W.I.36 i BURSA } 

SOLUTION FOR 3 TRANSLATION, I SCALE AND 3 ROTATION PARAMETERS 

fUSING VARIANCES ONLY) 

DU DV DW DELTA OMEGA PSI EPSILON 

METERS METERS METERS 1X1. D+6) SECONDS SECONDS SECONDS 

“AO. 22 -33.38 11.63 1.56 -0.66 -0.73 -0.53 

♦8.82 *15.24 *9.34 +1.35 +0.37 *0.31 +0.45 

VARIANCE - COVARIANCE MATRIX 


S02« 1.04 


0.7780*02 

0.102D•^02 

-0.4810*01 

-0.8050-05 

-0.5380-06 

0. 9000-05 

-0.150D-05 

0.1020*02 

0.2320*03 

-0.368D+02 

-0.392D-06 

0.2230-04 

0*4620-05 

-0.298 D-04 

•0.4810*01 

-0.3680+02 

0.8720+02 

-0.8160-05 

-0.4 190-05 

-0. 1020-04 

0.5 600-05 

0.8050-05 

-0.3920-06 

-0.8160-05 

0. 1820-11 

0. 2090-12 

0.505D-13 

-0,7440-14 

0.5 380-06 

0.2230-04 

-0.4190-05 

0.2090-12 

0.3190-U 

0.4730-12 

-0.198D-11 

0.900D-05 

0.4620-05 

-0.1020-04 

0.5050-13 

0.4730-12 

0.2230-11 

-0.5950-12 

•0.1 500-05 

-0 .2980-04 

0.5600-05 

-0.7440-14 

-0.1980-11 

-0.5950-12 

0.4830-11 



COEFFICIENTS OF CORRELATION 


• 

O.l 000*01 

0.7550-01 

-0.5840-01 

-0.6770+00 

-0.342D-0I 

0,6840*00 

-0.772 D-01 

0.7 550-01 

0. 1000*01' 

-0.2590+00 

-0,1910-01 

O.810D+OO 

0,2030*00 

-0. 8900+00 

•0.5840-01 

-0.2590+00 

0.1000*01 

-0.6490*00 

-0.2510*00 

-0.7360*00 

0.2730*00 

■0.677D+00 

-0.1910-01 

-0.6490*00 

O.lOOD+01 

0,8680-01 

0.2510-01 

-0 .2510-02 

•0.342D-01 

0.618D + 00 

-0. 2510+00 

0.868D-01 

O.lOOD+01 

0. 1770+00 

-0.504D+00 

0.6 840*00 

0.2030*00 

-0,7360*00 

0.2510-01 

0,1770*00 

0,1000*01 

-0.181D+00 

-0.772D-01 

-0. 8900+00 

0.2730*00 

-0.251D-02 

-0.5040*00 

-0.1810*00 

0.1000*01 


NOTE : THE POSITIVE AXES ARE TOWARDS KEAN 

GREENWICH, EAST, AMD CIO. THE ROTATIONS 

PRINTED ARE ABOUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY. 
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APPENDIX I 


CARTESIAN AND GEODETIC COORDINATES 
(Solution WEST35) 


Sta.No 

u 



Or 

W 

Or 


<0 

— 

mSSSm 


H 

Om 


a. A. r. 




u,v,w Cartesian coordinates in meters (Orientation: u h the Greemvich 
meridian as defined by the B.I. H. ; v = X = 90“ (E); w s Conven- 
tional International Origin). 

Geodetic latitude and longitude in angular units (degrees, minutes 
and seconds of arc) computed from the Cartesian coordinates and 
referred to a rotational ellipsoid of a = 6378155. 00 m and 
b = 6356769.70m. 

H Geodetic (ellipsoidal) height in meters referred to the same 

ellipsoid. 

Ou,0v.<T¥ Standard deviations of the Cartesian coordinates in meters. 

a^,o^ Standard deviations of the geodetic coordinates in seconds of arc. 

Cm Standard deviations of ^he geodetic height in meters. 


a*,A,,r, Altitude (elevation angle), azimuth and magnitude of the major 

semi axis of the error ellipsoid, respectively. Angles in degrees, 
magnitude in meters. Altitude is positive above the horizon. 
Azimuth is positive east reckoned from the north (see section 
4.74). 

at,A),,r^ Same as above for the mean axis of the error ellipsoid, 

Ac, Tc Same as above for the minor axis of tlie error ellipsoid. 
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6006 

2102928 

.90 

2.30 

721666.28 

2.68 

5958181 .55 

2.86 



xa 

0,07 

18 56 26.82 

0.25 

112.26 

2.92 


{{■■I 

87.35 

21.30 

2.92 





“ 0.55 

99,26 

2.68 






2.59 

“170 .76 

2.22 



6016 

4896388 

.39 

1.84 

1316170.88 

2.28 

3856667.89 

2.28 


37 26 39 

.09 

0.07 

15 2 44.55 

0.09 

15*36 

2.12 



48.56 

“ 7.14 

2.30 

% 





12.47 

97.37 

2.28 





38.73 

“162 .42 

1.81 



6065 

4213563 

.66 

2.05 

820824.47 

2.20 

4702785.16 

2.46 


47 48 4 

.56 

0.07 

11 1 24.46 

0.11 

958.82 

2 !. 32 



51.48 

14.04 

2.47 






“ 1.24 

102 .48 

2.19 





38.50 

“ 168.51 

2.05 



800 ^ 

3818505 

.10 

5,83 

708051,35 

5.36 

5042642.31 

9,52 


52 35 3 

.06 

0.22 

10 30 17.40 

0.29 

84.18 

8.89 



59.94 

20.97 

9.72 





29.93 

- 153.30 

5.78 






2.48 

115.27 

5.05 



8009 

3923401 

.40 

4.81 

299883.04 

5.23 

5002973.63 

5.35 


52 0 6 

.34 

0.16 

4 22 15.14 

0.27 

45.73 

5.21 




24.70 

117.19 

5.46 






51.98 

“ 8.85 

5.35 






26.90 

- 139.32 

4..54 



8010 

4331299 

.42 

4,19 

567503.65 

4.95 

4633112.83 

4.90 


46 52 37 

.21 

0.15 

7 27 52.54 

0.23 

920.22 

4.65 




20.94 

123.95 

5.06 





47,95 

59,06 

4.81 






34.50 

“ 161.30 

4.16 
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ttoil 

3920174.14 

3.62 

- 134743.75 

2.69 

5012725.51 

4.77 


52 8 3 5.61 

0.12 

358 1 53.08 

0.14 

142.97 

4.72 



62.16 

31.94 

4.99 





27.41 

“ 158.96 

3.63 





“ 4.50 

113.37 

2.26 



6015 

4578318.14 

3.79 

457961 .82 

4.31 

4403193.46 

4.15 


43 55 57.84 

0.13 

5 42 43.93 

0.19 

676.73 

3.82 



15.79 

111.05 

4.38 





26.72 

12.8 7 

4.21 





58.30 

“ 131.71 

3.64 



8016 

4188649.22 

3.89 

571417.93 

4.19 

4760143.75 

5.55 


48 34 58.70 

0.15 

7 46 6.13 

0.20 

166.04 

5.04 



53.17 

11.96 

5.58 





11.15 

117.22 

4.21 





34.57 

-144 .98 

3.83 



8019 

4579462.47 

4.25 

586586.81 

6.97 

4386418.77 

4,47 


43 43 33.26 

0.14 

7 17 57.52 

0.31 

394.83 

4.50 



8.27 

98.17 

6.98 





61.04 

“ 7.07 

4.50 


. 



27.54 

“ 167.49 

4.20 



8030 

4205632.55 

6.26 

16 3702.89 

8.02 

4776541 .56 

5.77 


48 48 22.10 

0.19 

2 13 44.74 

0.39 

187.76 

6.26 



“ 9.34 

85.27 

8.09 





65.15 

154 .46 

6.30 





22.81 

“ 0.76 

5.63 



8031 

3593853.84 

10.02 

“2027 58.8 8 

6.51 

5248053.62 

12.13 


55 44 0.93 

0.38 

356 46 15.23 

0.39 

291.46 

10.20 



27,44 

“ 11.68 

12.82 





57.64 

133.35 

9.67 




“ 15.80 

69.87 

5.65 



8032 

4213542.86 

2.06 

820771.07 

2.20 

4702810.23 

2.47 


47 48 5.84 

0.07 

11 1 22-13 

0.11 

956.82 

2 . 3 ? 



51.51 

13 .9 8 

2.48 





“ 1.28 

102.37 

2.19 





38.46 

“168 .65 

2.06 
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8033 

4041868.23 

6*26 

620631.85 

6.27 

4878635.20 

7,25 


50 13 11.41 

0.20 

8 43 46.66 

0.30 

196.71 

7.51 



69.61 

74.59 

7.74 





- 6,36 

147.13 

6.44 





19,28 

- 125.11 

5.45 



8034 

3919692.35 

4.81 

298829.50 

5.23 

5005899.45 

5.35 


52 2 40,96 

0,16 

4 21 34.86 

0.27 

26.98 

5.22 



24,69 

117.20 

5.46 





52.01 

- 8.87 

5.35 





26.87 

- 139.33 

4.54 



8701 

4194418.40 

5.16 

1162693.39 

2.53 

4647182.33 

7.61 


47 4 0,48 

0.17 

15 29 36.17 

0.12 

469,19 

7.63 



68.28 

10.34 

8.00 





21.69 

- 166.55 

4.69 





1.07 

103,02 

2.20 



8702 

4027920.76 

5,37 

307002.41 

7.73 

4919441.96 

10.97 


50 47 50.89 

0,27 

4 21 30.87 

0.39 

121.29 

9.09 



46.93 

- 15,98 

11.07 





22.82 

100.76 

7.81 





34,20 

- 152.62 

5.03 



8703 

3513637.60 

9.58 

778937.20 

6.33 

5248201.26 

10.75 


55 44 20.34 

0.25 

12 29 59.02 

0.40 

61.20 

11,76 



82.06 

- 93.78 

11.83 





- 4.82 

- 146.56 

8.59 





6,30 

123.97 

5.80 



8 705 

4530504,22 

21.34 

- 41723.13 

33.91 

4474399.67 

45.54 


44 50 3.23 

1.06 

359 28 20.48 

1,54 

119,77 , 

38,37 



47.38 

15.76 

46.13 




- 26.71 

72.62 

37.81 





30.43 

145.4 3 

21.02 



8706 

4587887.69 

3.79 

419518.77 

4,31 

4396450.97 

4.15 

‘ 

43 51 8.74 

0,14 

5 13 28,65 

0,19 

226,89 

3.82 



15.45 

111.03 

4.38 





26,67 

13.06 

4,21 





58.54 

-132 .12 

3.64 
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8710 

3818503.32 

5.82 

708052.67 

5.36 

5042642 .47 

9.52 


5? 35 3.10 

0.22 

10 30 17.49 

0.29 

83.39 

8.89 



59.95 

20,93 

9,72 





29.93 

-153.24 

5.77 





2.52 

115,30 

5.05 



8711 

4896386.79 

1.85 

13161 70. A 4 

2.28 

3856670.06 

2.29 


37 26 39.18 

0.07 

15 2 44.55 

0.09 

15.36 

2.13 



48.35 

-8.09 

2.31 





13.01 

96 .97 

2.28 





38,70 

-162,37 

1.82 



8712 

4335514,42 

29.88 

1063081 .54 

30.79 

4540929,32 

34.88 


45 40 55.85 

0.90 

13 46 38.20 

1.35 

389.55 

37.76 



66.09 

-46.11 

38.59 





23.88 

130.93 

33.19 





1.10 

-138.59 

21.61 



8713 

462862 0.82 

61.20 

1471951.97 

40.91 

4120480.49 

56.71 


40 29 58.82 

1.17 

17 38 28.44 

1.98 

211.37 

71.77 



61.48 

-75.02 

79.11 





23.40 

142.20 

41 .65 





15.37 

45.3? 

25 .30 



8714 

4885390,86 

30.59 

784057.56 

31.96 

4011490.51 

59,59 


39 13 17.86 

1.58 

9 7 3.58 

1.19 

112.75 

48.13 



44.19 

-0.69 

59.75 





34.53 

131 .30 

40.23 





26.04 

-119.06 

17.88 



8715 

4896390.83 

1.86 

1316178.67 

2.29 

3856662.17 

2.29 


37 26 3 8,85 

0.07 

15 2 44.83 

0.09 

15.36 

2.14 



48.23 

-8 .6 4 

2.31 





13.32 

96 .74 

2.29 





38.68 

-162.33 

1.83 



8716 

4850690,58 

11.76 

-315908.83 

7-56 

4116651 .39 

16.72 


40 27 1,63 

0.4 A 

356 16 25.63 

0.31 

717,89 

15.53 



53.53 

11.78 

17,16 





35.36 

-151 .99 

12.23 





7.79 

112 .44 

5.58 




75 


8717 

4850600.69 

11.77 

- 315911.54 

7.56 

4116651.11 

16.72 


40 27 1 .62 

0.44 

356 16 25.51 

0.31 

717.88 

15.53 



53.54 

11 .78 

17.16 





35.35 

- 152.00 

12.23 





7.78 

112.44 

5.58 



8718 

4146546.32 

10.32 

6131 15.18 

7.59 

4791501 .83 

11.51 


49 0 3 9.73 

0.25 

8 24 39.26 

0.34 

165.70 

13.71 



76.05 

91.42 

14.04 





“ 1.46 

175.51 

7.75 




“ 13.87 

85.15 

6.29 



8719 

4883056.88 

1.86 

1306097.78 

2.29 

3879628.98 

2.29 


37 41 35.53 

0.07 

14 58 28.82 

0.09 

1740.87 

2.14 



48.55 

“ 8.39 

2.31 





13.08 

96.87 

2.29 





38.47 

- 162 .49 

1.83 



8720 

3104216.61 

14.17 

9983 56.71 

8.62 

5463292 .92 

12,40 


59 20 16.45 

0.33 

17 49 42.49 

0.68 

58.50 

14.38 



52.52 

- 130.78 

15.67 





37.41 

45.12 

11.88 





1.98 

136.6 4 

6.49 



8722 

2591994.43 

10.90 

- 1078496.29 

18.45 

5707868.71 

30.88 


63 57 42.58 

0.47 

337 24 30.31 

1.36 

5.50 

29.38 

- 


70.01 

0.15 

31.04 





- 3.59 

80.22 

18.72 





19.64 

168.94 

9.96 




NORMAL TERMINATION 
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Appendix II 
Transformation 

WEST 35 -TC- WEST 5A( MOL vnOFL ) 

SOLUTinrj EDR 3 TRANSLATinN', 1 SCALE A^p ?, POTATirni PARAMET»=RS 

(USING VARIANCES ONLY) 


DU 

DV 

DW 

DELTA 

OMEGA 

PSI 

EPSILON 

METERS 

METERS 

M FT ER S 

txi ,r+6) 

SECONDS SECONDS 

•seconds 

in 

• 

o 

0.23 

-0.23 

“0,01 

0.5? 

-0,11 

-0,70 

1 0.5.5 

± 0.f>2 

+ O.AP 

i 0.52 

+ 0.20 

+ 0,13 

♦ 0.2^* 


VARIANCE - COVARIANCE MATRIX 


SO?= 0.23 


0. 3350+00 

-0.654D-02 

0.393D-02 

0. A? CD -07 

0,8630-07 

-0.1860-07 

-0 ,7910-07 

-0.(SS4D-02 

0,3870 + 00 

-0.530D-02 

0.3270-07 

-0 . 12^0-06 

0.1800-07 

0,6^90-07 

0.393D-02 

-0 ,53on-o? 

0.467D+O0 

-0.1820-07 

0 .6310-07 

0.452D-07 

-0,135 0-06 

0.420D-07 

0.327D-07 

-0.1P2D-07 

0.2600-12 

-0. 1710-13 

-0.3950-13 

0.203D-13 

0.863D-07 

-0,1200-06 

0,6510-07 

-C.lTin-13 

0.9520-1? 

-0. 1960-12 

-0,991.0-12 

-0.186 D-07 

O.lROD-07 

0.4520-07 

-0.3950-13 

-0,1960-12 

0.392D-12 

0.236D-I2 

-0.791 D-07 

0. 6700-07 

-0.1350-06 

0.203D-13 

-0.8910-12 

0.2360-12 

0.162D-11 



COEFFICIENTS OF CORRELATION 

• 


O.lOOD+01 

-0.1900-01 

0.9O3D-02 

0.1400+00 

0.1530+00 

-0.513D-01 

-0.1 070+00 

-0.1 90D-01 

0.1000+01 

-0.1250-01 

0.1010+00 

-0.2130+00 

0.4630-01 

0.8500-01 

0.9«>3D-02 

-0.125D-01 

0.1000+01 

-C.514D-01 

0.9460-01 

0.1060+00 

-0. 1560 + 00 

0.140D+00 

O.IOID + OO 

-0.5140-01 

ionn +01 

-0.3390-01 

-0.1220+00 

O .3070-01 

0.153D+C0 

-0.2130+00 

0.Q46D-01 

-0. 3390-01 

0,1000+01 

-0,3200+00 

-0.7180+00 

-0.513D-01 

0 ,4630-01 

0.1060+00 

-0,1220+00 

-0,3 200+00 

C.lOOD+01 

0.2O7D+00 

-0.1 07D+00 

0 .8590-01 

-0,1560+00 

0.3070-01 

-0.7180+00 

0.2 970 +00 

o.ioon +01 


NOTE : THE POSITIVE AXES ARE TOWARDS MEAN 


GREENWICH. EAST. AND CIO. THE ROTATIONS 


PRINTED ARE AF.OUT 3RD, 2ND, AND 1ST AXES RESPECTIVELY, 


77 




RESinUALS V * 


VI ( WEST 351 


V?t WEST 3^J 



6C06 

“C. 1 

“0,4 

“C.l 

6006 

C. 1 

0.4 

0.1 

“0.2 

-O.P 

“0.2 

6004 

0.6 

-0.6 

“1.5 

8004 

-0,4 

C.6 

0.4 

1.0 

“1 .2 

“1.9 

R0C9 

“0.4 

0.0 

“0.2 

8 009 

0.3 

“0.0 

0.1 

“0.6 

o, 1 

-r.,4 

8010 

0.4 

”0.1 

1.0 

8010 

“0,3 

0. 1 

“0.6 

0.7 

-0,3 

1.6 

8011 

0.4 

“0.0 

0.2 

801 1 

-0.3 

n.o 

-n, 1 

0.6 

-0.1 

0,3 

8016 

“0.2 

0. 1 

O.P 

8016 

0 . 1 

“0. 1 

“0. 4 

-0,4 

0.2 

l.l 

8019 

C.5 

“0.1 

0.7 

8 019 

“0, A 

C.l 

“0.5 

0.9 

“0.3 

1.2 

8030 

“0.7 

“0.2 

0.0 

8C?Ci 

0.4 

0.2 

“0.0 

“l.l 

“C.3 

0.0 

8031 

4,0 

1.1 

7.9 

8 031 

“3.6 

“1 .0 

“2.6 

7,6 

2.1 

ir.5 

803? 

0.2 

0.0 

0.4 

8032 

“0,2 

-0.0 

“0,3 

0.4 

0.1 

0,6 

8023 

“2.0 

“0.8 

“3.5 

8 033 

1.2 

0.7 

1.4 

“3.2 

“1.5 


8701 

4.2 

1.7 

12.5 

8.701 

“2.4 

“1.3 

“3.8 

6.6 

3.0 

16.3 

8702 

“C.3 

“0.2 

-2.9 

8702 

0.2 

G.2 

0. 7 

“0,6 

“0.3 

“3.6 

8703 

“2.4 

“0.5 

-5.9 

870? 

1.3 

0.5 

1.3 

“3.7 

“0,0 

-7,2 

8 705 

2.7 

-4,8 

-12.1 

8 705 

“1.8 

4.6 

1.9 

• 4.5 

-0,4 

“14.0 

8706 

0.5 

0.0 

0.8 

8 706 

“0.4 

“0.0 

“0.6 

I .0 

0.0 

1.4 

8711 

“0.3 

“O.e 

“0.6 

8711 

0.3 

0.8 

0,5 

“0,6 

“1.5 

-1,1 

871? 

2.P 

0.7 

2.1 

8712 

“1.0 

-0. 7 

“0,6 

3.8 

1.4 

2.7 

8713 

“11. C 

1,3 

“11.3 

8713 

1.1 

“1.2 

2.6 

-12.1 

2.4 

-13.9 

8714 

6.3 

4.2 

28.^ 

E7J4 

-4,0 

“3.6 

“8.4 

13.3 

7.9 

35.8 

8716 

-7.9 

0,4 

-11.2 

8716 

4,8 

“0. 4 

8.4 

“12.7 

0,7 

-19, 7 

8 718 

“9.5 

“3.1 

“11.9 

8718 

2.6 

2.3 

1.3 

“12.1 

“5.4 

-13,3 

8720 

“8.3 

1.2 

“11.5 

8720 

4.6 

“1.2 

2.4 

“12.8 

2.5 

-13.8 

8 722 

2.3 

C.3 

1.1 

8722 

“2.0 

“0.3 

“0.0 

4.3 

0.5 

1.1 


UNIT OF RESIDUALS (METERS) 
* Residuals in the Cartesian coordinate system. 
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REPRODUCIBILITY OF T. 
ORIGINAL PAGE IS POO 


RFSIPIIALS V 
SPHEkICAL * 


VI { WEST 35) 


V2( WEST 34) 


VI - VP 


6006 

0.17 

“ 0,34 

- O '. 2 

6006 

- 0,17 

0.34 

0.1 

0.34 

- 0.67 

- 0,3 

eoc4 

- 1,28 

“ 0.73 

- G.o 

8004 

0.48 

0,65 

C.i 

- 1.76 

“ l.?r 

- l.G 

so 09 

0.13 

0,07 

- 0.4 

8009 

“0 .12 

“ 0.06 

0.3 

0.25 

0 . 1 ? 

- 0.7 

8010 

n .41 

-<^.20 

1.0 

8^10 

- 0.18 

0.18 

- 0.6 

0 . 5 ® 

- 0.38 

1,6 

8011 

- C .16 

“ 0.02 

0.4 

soil 

0.14 

0,02 

- 0.3 

- 0.29 

“ 0,05 

o . 7 

8016 

0.66 

0.13 

C .4 

8016 

“ 0 . 3 ? 

- 0.11 

- 0.2 

0.99 

0.24 

0,6 

8019 

0.14 

“ 0.18 

0.8 

8019 

“ 0.11 

If . 1 6 

“ 0.6 

0,25 

“ 0,34 

1.4 

8 030 

0.55 

- 0,14 

- 0.5 

8030 

- 0 . 3 1 

0,15 

0.3 

0.87 

- 0 . 2 ® 

-C,7 

8031 

1.19 

1,29 

E .P 

8031 

1 .42 

- 1.24 

-4 .1 

“ 0,23 

2.54 

12.9 

8032 

0,07 

- 0,00 

0.4 

8032 

- 0 . 0 ? 

- 0,00 

“ 0.3 

0.10 

- 0.00 

O.c 

80 33 

- 0.62 

“ 0.46 

“ 4.1 

F 033 

- 0.13 

G .51 

1 

- 0 . 49 

- 0.®7 

- 6,0 

8701 

5,17 

0,48 

12.2 

8701 

“ 0.65 

“ 0.61 

“4 .6 

5 . P 3 

1 .OQ 

16.8 

8702 

- 1 , 5 ® 

- 0.14 

- 2.5 

8702 

0.24 

0.15 

0.7 

- 1,83 

- 0 . 2 ® 

- 3.2 

8703 

- 1.29 

0.06 

- 6.3 

8703 

- 0.43 

0.17 

1.8 

- 0.86 

- 0.12 

-8.1 

8705 

- 10.54 

- 4.78 

- 6.6 

870 5 

2 .65 

4.58 

0.0 

“ 13,19 

- 9.37 

- 6.6 

8706 

0,21 

- 0.05 

0.9 

8706 

- 0,16 

0.04 

- 0,7 

0.37 

“ 0.09 

1.7 

8 711 

- 0.15 

- 0.66 

- 0.7 

8711 

0.13 

0.66 

0.7 

- 0.28 

- 1.33 

- 1.4 

8712 

- C .65 

0.01 

3.5 

8712 

0.35 

- 0.42 

- 1.2 

“ 1 . 00 

0 . 4 ? 

4 . F 

8713 

- 2.05 

4.55 

- 15.0 

871 ? 

1.53 

“ 1.45 

2.2 

“ 3 . 5 ® 

6.00 

- 17,3 

8714 

16.42 

2.87 

24,9 

8714 

-3 .04 

“ 2.78 

- 9.5 

19.46 

5.65 

34,4 

8716 

- 3.43 

- 0,15 

“ 15.3 

8716 

3 .28 

- 0.04 

9.2 

“ 6,71 

“ 0,11 

- 72.5 

8718 

- 0.41 

“ 1.70 

- 15.5 

£718 

-1 .31 

1.91 

2.9 

0.00 

“ 3.61 

- 18.3 

8720 

0.60 

3.71 

“ 13.7 

8720 

- 2.22 

- 2.56 

4.1 

2.82 

6,27 

- 17.7 

8722 

- 1.34 

1.11 

1.8 

8722 

1.54 

- l.OO 

- 0.8 

“ 2.87 

2.11 

2.6 


UNIT OF RESIDUALS (METERS) 


* The spherical residuals ropresent the residuals in the Cartesian coordinate system as 
transfoi-med into the curvilinear (<o , X , h ) system . 
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PlffiCEDING PAGE BEANK NOT riEME 
Appendix ni (Referenced in Section 5) 

Part 1 

A DISCUSSION ON TRANSFORMATION MODELS 
1 . Introduction 

The availability of geocentric station coordinates obtained by satellite geodesy 
made it possible to compute the relative positions of the ellipsoids which are used in 
classical triangulations with respect to an earth centered ellipsoid. The coordinate 
system of the latter is the Average Terrestrial system (AT), where the z axis is 
directed toward the average north terrestrial pole as defined by the Internationa] 
Polar Motion Service (IPMS) and is commonly known as the Conventional Interna- 
tional Origin (CIO). The zx plane is parallel to the mean Greenwich astronomic 
meridian as defined by the Bureau International de I’Heure (BIH). The classical 
triangulation is computed in the geodetic coordinate system (u,v,w), with the w 
axis directed toward the North Pole and the wu plane coinciding with the geodetic 
Greenwich meridian. 

Several procedures for transformations have been published in recent years. 

A summary of the three most commonly used models is given in [Badekas, 1969]. 
Some clarifying remarks about these models are given here. One additional trans- 
formation model, although intended for a different purpose, has been suggested by 
[Vanicek, 1974 and 1975], 

2. Transformation Mbdel 1 (Bursa) 

This model was introduced in [Bursa, 1962 and Wolf, 1963] . It treats the 
Cartesian coordinates as observations and the seven transformation parameters, 
(three shifts, one scale and three rotation angles), as quantities to be solved for 
in a least squares adjustment. 
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Bursa Model 


The notation is as follows: 

X,Y,Z ; Averj^e Terrestrial coordinate system (satellite system) 
with origin at 0. 

U, V, W : Geodetic coordinate system with origin at 0 
Pj : Terrain point. 

pQ % Shift vector between the origins 0 and O'in the satellite system, 

rj : Position vector of Pj in the geodetic system, 

p^ ; Position vector of Pj in the satellite system. 

The transformation of point Pj is described in the satellite sj ihm b3’- the equation 

Fj = Po + (l + k)R?j - = 0 (1) 

where R is the product of three orthogonal rotation matrices and k denotes the 
scale factor. Because only differential rotations are being considered, the sequence 
of rotations is irrelevant. Deleting small terms of the second order and mixed 
teims, one arrives at a rotation matrix of the form 
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10 0 


' 

1 

3 
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1 0) -0 

R s I + Q = 

0 10 

+ 

-w 0 c 

= 

-to 1 e 


0 0 1 


Ip -€ 0 


0 -e 1 


The rotation angles 03, 0 , c , when positive, represent counterclockwise rotations 
about the respective W, V, U axes, as viewed from the end of the positive axis . If 
the term kQrj in Eq.(l) is ignored, one obtains the expression 

Fj = po+ Tj - pj +krj + Qrj = 0 (3) 


which forms the mathematical model for a least squares solution. Here, the 
model [Uotila, 1967] used is of the form 

F(La, Xa) = 0 (4) 

where La denotes the adjusted observations and Xa denotes the estimates of 
parameters . The linearized form gives 

BV + AX + W = 0 (5) 


in which B = 9 F/S La , A = B F/S Xa and W = F (L^ , Xo ) . Lb denotes the 
observations, e.g. , the Cartesian coordinates in both the satellite and the 
geodetic coordinate system (pj and r^ ) and Xo denotes the approximate 
values for the parameters (po k, 60 , 0 , c ). That part of Eq. (5) which pertains 
to point Pj can be written as 


100-100 
0100-10 
0 0 1 0 0 -1 


Vv 

V„ 

Vx 

Vy 

V, 


A, 

1 0 0 u V -W 0 
0 1 0 V-U 0 w 
0 0 1 W 0 U -V 



X 

4* 

Wj 


DX 


u - X 

• 

DY 

4* 

V - Y 

j 

DZ 


w- z 

k 




to 




0 




€ 




= 0 


( 6 ) 
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3. Transformation Model 2 (Molodenskii) 


The model described in this section is attributed to Molodenskii. At first 
glance it looks very much like Bursa’s model as expressed in Eq. (1). One may 
also be inclined to use a figure similar to Figure 1 to interpret or even to setup 
the transformation equation. The purpose of the discussion here is to clarify 
this lype of misinterpretation. The Molodenskii transformation equation is 
commonly stated as 


Fj = Po+ Tk + <l+k)R^^j ~ Pi = 0 


where Po 


-> 

ric 


R 

rkj = rj-r;, 

Pi 


; is usually interpreted as the shift vector between the 
origins of the two ellipsoids in the satellite system. 

% is the vector of the initial point in the geodetic system. 

It is assumed that the components refer to a coordinate 
system that is already parallel to the satellite system, 
hi practical computation, however, it is simply the 
vector of the initial point expressed in the geodetic system. 
Thus the transformation of from geodetic to satellite 
system is ignored. 

: is the same as in Eq. (2). 

; is the relative position vector of the point Pj with respect 
to the initial point in the geodetic system. 

: is the position vector of Pj in the satellite system. 


If the term kQrj^ j is deleted from Eq. (7), then one obtains the new expression 
Fj = Po + - Pi + = 0 . 


The analog to Eq, (6) reads 




4. Transformation Model 3 (Veis) 

This model (usually referred to as the Veis Model), is identical to Model 
No. 2 (Molodenskii model) except for the axes about which the rotations occur. 

The point of rotation is again the initial point of triangulation [Veis, I960]. The 
first axis is tangent to the geodetic meridian with positive direction toward the South; 
the second axis is perpendicular to the meridian plane and is positive eastward and 
the z axis. is along the geodetic normal with its positive direction upward forming a 
right-handed system with axes one and two. Thus, the only difference between the 
Molodenskii and Veis models is that the rotations in the Veis model are familiar 
quantities; for example, a rotation about the third axis corresponds to a rotation 
in azimuth . 


5, Comparison between Models 1 and 2 (or 3) 


rT represents the vector of "translation” as obtained by the various 
models. A simple comparison between Eq. (1) and (7) indicates that 


ri = Po 

(bursa ) 

Tt = Po- IV + (1+k) Rrk . 

(KOUOOENsKl I ) 
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Eq. (10) is very important as far as the geometric interpretation is concerned. 
Clearly, r^ is not the vector between the origins of the two coor- 

(MOLOOENsK II ) 

dinate systems, but rather a function of the scale, the rotations and the choice 
of the point of rotation rjt. Eq. (10) is based on the knowledge that both 
models give identical values for the scale and rotation parameters as shown 
in the following paragraphs. In both cases the same adjustment model is 
used. The only difference is the design matrix A, as can be seen from 
Eq. (6) and (9). If, respectively, Af, and A denote the design matrices of the 
Molodenskii and Bursa models and they are partitioned into (3 x 3) and 
(3 X 4) submatrices, then one has 



r 

I 

Af/|x 


I 

Ax " 

D 


Am ~ [ I 1 Am ] 

I 

Am 2 

= 

I 

A2 " 

p 

(11) 


i 

Am r 


i 

A,- 

D 



where r is the number of stations to be transformed. Both design matrices differ 
only by the submatix D, which is the same for every station. If Eq. (5) is partitioned in a 
similar manner, and assuming there is a block diagonal variance -covariance matrix 
(no correlation between the coordinates of different stations), the normal equations 
can be conveniently solved using the general formulas for the inverse of partitioned 
matrices [Uotila, 1967]. Generally, if 
Nxx Nx 2 

” " Ni's Ns3 

has full rank and Nxx Nsa ai'e square and non-singular, then 


Nix Nx2 
N22 


Ni^+ (N 2a“Ni2 Nii_N^)>i^Nu_;' . 

- (N22- Nx'sNii'^Nxs)'^ N4 Nu" 


,1 (N32 

1 (N2 2-N^Nu-"Nx2)-' 


(12) 

If one minimizes the norm V'PV subject to the condition (5), one obtains the 
well known equations 
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X 


( 13 ) 




V = -P'^B'M'^(AX + W) 


where 


M = 


( 14 ) 


( 15 ) 


and P'^ is the variance-covariance matrix of observations which is taken to be 
block diagonal . Since B has the form indicated in Eq.(6), M has the form 

Pi 0 


M = 


0 


P'r 


where Pj = *V[ S g + Lg ] 

Oh * 


Ob 


Qo being the a priori variance of unit weight and Eqj and Egj being the variance- 
covariance matrix of geodetic and satellite coordinates, respectively, at point i. 
With the notation 

S = [SA/PiAj -(SPi Ai)'(LP/(SPt Ai)]\ 

one can express the normal matrix of the Molodenskii model with the help of Eq. 
(12) as 


( 16 ) 


(A^tM•U)^ H 


I ~l 

(SPi )•' + (LP,)-"(SPi A,- SPiD) S(EPj A,-SPiD) '(SP*)'' ■ -(5 :p,)\LP,Ai- EPjD) S 


-S (SPA-SPiD)^(SPi)-^ 


( 17 ) 


all summations being taken over r, where r is the number of points to be transformed. 
Continuing the solution of Eq. (13), one obtains 

Xmi = [(SPi)‘^+(SPj)-'(SA/Pi)' S (SA/ Pi)(SPi)-'](SPiWO 
-(SPi)-^(SA/Pi)'S (LA/P^Wi) 

+ DS [ (SA/Pj Wi) - (LA: Pi) (LPi)’' (LPiWi ) ] . (18) 

Xa = S [S(A/PiWi) -(SPiAi) '(SPi)'"(SPiWi)] . (19) 


In Eq, (18) the term within the last brackets [] is not eqiial to zero. The 


87 


parameter vector X is partitioned into Xi and Xg . The subscript M distinguishes 
the solution in the Molodenskii model from the solution in the Bursa model. According 
to partitioning, X^i contains the translation parameters and Xf^s contains the scale 
and the rotation parameters. It is readily seen from the above equations thatXMi 
does depend on the submatrix D and X^s does not. In fact, the first two terms of 
Eq.(18) represent exactly the solution which is obtained when using the Bursa 
model. 

Next it is proven that the product AX is also independent of the submatrix D. 

AX is invariant in the two systems: 


=> 


- a!A.-D)Pj 

( 20 ) 


VXg/ 

VXg / 


=> 

Xi+AjXs 

- Xmx+AiXs-DXs 


=> 

Xx 

= X,i-DXg 

( 21 ) 


However, from Eq. (18), we have 
X„i = Xi+DS[], 

where the brackets [] denote the last term on R.H.S. of Eq.(18). 

X„i -DXg = Xi+ DS[ J-DXg . (22) 

If XgOf Eq. (19) is substituted in Eq.(22), it is then seen that the last two terms 
of Eq.(22) cancel, thus proving the identity (20). Therefore, the products 

AmXm and AX 

are invariant, and as a result both transformation models yield the same residuals. 
Note also that Eq. (21) is equivalent to the second equation in (10) — considering 
that the mixed (small) terms were neglected in the mathematical model of the 
adjustment. 

The final step is proving that the direct transformation from one system to 
the other yields identical results when either of the two solution vectors and its 
correspondiag variance-covariance matrix is used. It has already been proven that 
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( 23 ) 



holds where the subscript T denotes the point which is to be transformed. 
Also, r|. is the position vector of point T in the new system; consists of 
the componenets of the station vector in the old system and it is just a matter of 
matrix multiplication that shows the variance-covariance matrix is indeed 
invariant. One has 



(IAmt) 



(24) 


Since 

Sx^ ~ , (25) 

one can use the identity (17) to carry out the operations necessary in Eq. (24). 

The result is 

Ex = (SP,)-^+(i:Pi)-^(EPiAi)S(EPiAi)'(EPi)-" 

T 

-2(EP0^SPiAi)SA/ + AySA/ . (26) 

Once again, Eq. (26) is independent of the submatrix D. It is identical to the 
expression obtained when substituting Ay (instead of Amj) in Eq. (24). Even if 
variances are attached to the coordinates to be transformed, the complete error 
propagation yields a variance-covariance matrix which is independent of the 
submatrbc D. 

Conclus ions: 


1) It has been shown tliat the so-called Molodenskii and Bursa ti-ansfor- 
mation models give the same scale factor, rotation angles and residuals. Never- 
theless, they give different translation parameters. It has been shown that 
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while the Bursa model immediately gives geometrically meaningful shifts between 
the origins of the two Cartes i^ coordinate systems, the translation parameters 
Computed from Model 2 need to be suitably modified to give the same shifts. For 
example, the translation parameters of Table 10 are related by Eq, 10 or Eq. 

21, respectively. 

2) It is true that Model 2 gives si^ificantly smaller correlation coefficients 
between its parameters. However, the correlation coefficients in Models 1 and 2 
cannot be compared since the meaning of the respective parameters is different, 

3) in the previous derivations no use was made of any specific pattern of 
the submatrix D. Therefore all the above conclusions are valid for any rotation 
point (intersection of axes about which the rotations are defined) . Model 2 (which 
specifies the initial point as the rotation point), therefore, is just one special 
case of an infinite number of possible models, differing only in the location of the 
rotation point. 

4) If one is interested in obtaining transformation parameters only to 
transform an arbitrary set of points from one system to another, no geometrical 
interpretation needs to be attached to the transformation parameters. Either 
model can be used since both give identical results. 

5) At this point the question that should be raised is whether or not the 
so-called Molodenskii model as treated in [Badekas, 1969] is, indeed, treated 
identically in [Molodenskii, 1962]. 

6) The first four conclusions listed here are also valid for the Veis model. 

7) When transforming a system of non-global coverage, it can be shown 
that a three parameter transformation yields shift vectors with magnitudes that are 
close to the ones obtained by a seven parameter transformation (Molodenskii). 

A three parameter transformation, therefore, should be used with great discretion, 
particularly when the systems have a large difference in scale. 
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Appendix HI 
Part 2 

RECOVERY OF SCALE FACTOR 
AND ROTATION ANGLES FROM CHORD COMPARISONS 

It is well toiown that the chord between two points expressed in terms of three- 
dimensional Cartesian coordinates is invariant with respect to a shift or rotation 
of the coordinate system. This makes it possible to independently determine the 
scale factor from chord comparisons and the rotation angles from comparisons of 
the directions of the chords. However, care has to be taken that all the correlations 
between various chords are included. In addition, for scale factor computations 
in a network of n points, one should use only those (3n - 6) chords which determine 
the three-dimensional configuration completely. Any additional chords do not con- 
tribute scale information, but instead are dependent on the other chords. This 
becomes clear if it is remembered that the given set of station coordinates merely 
determines a polyhedron uniquely. Since the chords are derived directly from the 
station coordinates, one can at most use only as many chords as are needed to 
uniquely determine the polyhedron. 

To appreciate how only (3n - 6) chords can be chosen to determine the scale 
factor, consider the case of 5 stations with coordinates known in the two systems. 

The station coordinates in each system are treated as observations with their 
variances and covariances. In obtaining a scale factor from chords, if one treats 
the chords as derived observations from the station coordinates of 5 stations, one 
can generate 10 chords. But if all these 10 chords are taken as observations, the 
variance -covariance matrix of observations becomes singular. This can be 
proven as follows: 

Let Xj , Yi , Zi (i = 1, 5) be the given coordinates of the 5 stations. 

Let Lk (k = 1,10) be the dei-ived chords. Then 

Lk = 7(Xj -X,)" + (Yj - Y,)" + (Zj - Zi)’^ i j , j > i . (27) 


91 



If Ll represents the variance-covariance matrix of the 10 cliords» then 


Sl = GExG' 


( 28 ) 


where Sx the 15 X 15 variance-covariance matrix of the station coordinates and 


BLv dLk aLk BLk bl^ 


G = 


axj aYj BZj axi aY^ azi 


has the size. 10 x 15 and takes the form 
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where 

AXk = Xj - Xj j > i 

AYk = Yj - Yi 
AZk = Zj - Zi 


(29) 


( 30 ) 


The order can easily be ascertained from a glance at the G matrix, e.g. * 
AXs = Xb - Xs . Although the matrix G is of the size 10 x 15, it has only a ranlv 
of 9 because the rows considered as vectors form a dependent set. Specifically, 
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when the chords are considered as vectors it can be seen, for example, that 




X5-X4 


Xs-Xi + Xi-X4 


AXv 


~AXe' 

AYio 

= 

Y5-Y4 


YS-Y1 + Y1-Y4 

= 

<1 

- 

AYb 

A Zio_ 


/Zb-'Z4 


_Zs-Zi + Zi-Z4_ 


AZ 7 _ 


AZs 


( 31 ) 


Thus, row 10 (chord 10) is a linear combination of row 7 (chord 7) and row 6 (chord 6). 
This was numerically tested and is also evident from the fact that 6 chords of 4 stations 
determine a three-dimensional figure uniquely and three more chords are adequate 
to fix the fifth point. 

Since the rank of G is 9, that of 'Ll can at most be 9, causing a ranlc deficiency 
of 1 in the 10 X 10 matrix Sl» irrespective of whether Ex is a full or diagonal matrix. 
Thus, in general, for n stations one can at most use (3n - 6) chords to compute the 
scale. Of course there are many sets of (3n - 6) chords which could be chosen to 
determine the scale. They are all equivalent as far as scale factor computation is 
concerned. 

The rotation angles can be computed separately in a similar manner since all 
the rotational information of the coordinates is also present in the chords. When 
giving the direction of a chord by two independent angles, (e.g. , Greenwich Hour 
Angle/ Declination), one should again use only as many directional angles as are 
needed to determine the shape of the polyhedron completely. The same computa- 
tional procedures are obtained when one eliminates in the transformation equation, 
the translations and rotations or the translations and scale, respectively. 

The subject matter of this Appendix is treated in greater detail in a separate 
report [Leick and van Gelder, in press]. 
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■ receding page BrANK NOT FIEMED 

Appendix IV 

FORTRAN IV PROGRAM MTH SUBROUTINES 
PROBLEM CODE DEFINIT-IONS 

COLUMN MEANING 

1. OVERALL PROBLEM CODE 

PC0DE(1)=1 MEANS OPTICAL PROGRAM,GEOMETR IC MODE, OSU FORMAT 

2 MEANS RANGE, geometric MODE 

3 MEANS SOLUTION ONLY RUM 

A MEANS ORBITAL MODE, OPTICAL OBSERVATIONS 

5 MEANS ORBITAL MODE, RANGE OBSERVATIONS 

6 MEANS ORBITAL MODE, MIXED OBSERVATIONS. 

PC0DE{1)=7 MEANS OPTICAL PROGRAM,GEOMETR IC MODE, GEOS FORMAT 

2« PERFORM SOLUTION? 

PC0DE(2)=1 MEANS YES 

0 MEANS NO 

PCODEm = 3 IMPLIES PCODE(2) = l 
3 . MAXIMUM NUMBER OF ITERATIONS? • 

PCODE(l) MUST EQUAL 1, 2, OR 7 
PC0DE(2) MUST EQUAL 1, 

PCODEI5) MUST EQUAL 1, FOR ONE OR MORE COMPLETE ITERATIONS 

5 . FORM NORMALS? 

PROCESSING CODES 

1 MEANS YES, 0 MEANS NO 

6. SIMULATE GUIDE MATRIX? 

7 . PRINT NORMALS? 

8. PERFORM SUMMARY BY OBSERVED LINES? 

9. PUNCH NORMALS IN ASD FORMAT? 

10. SUMMARIZE RESULTS 

PC0DE(10)=0 DO NOT PRINT SUMMARY 

=1 PRINT THE DX'S AND STANDARD DEVIATIONS 
=2 PRINTS THE X,Y,2*S AND STANDARD DEVIATIONS 
=3 PRINTS THE LAT I TUDE, LONGITUDE AND HEIGHT 
PRINTS BOTH X,Y,Z C LAT., LONG, & H 

11 . PRINT SATELLITE POSITION FOR EACH EVENT? 

0 MEANS NO 

1 MEANS PRINT XYZ AND GEODETIC COORDINATES 

2 MEANS PRINT XYZ ONLY 

3 MEANS PRINT GEODETIC COORDINATES ONLY 

12 . THIS PARAMETER DESCRIBES WHERE THE STANDARD DEVIATIONS OF THE 

INDIVIDUAL OBSERVATIONS (USED TO FORM THE WEIGHTS) ARE TO BE FOUND 
PC0DE(12)=0 MEANS TO READ THE OBSERVATIONAL STANDARD DEVIATION 
FROM THE CARD CONTAINING THE OBSERVATION. 

PCOOE(12)=1 MEANS TO ASSOCIATE A SINGLE STANDARD DEVIATION WITH 
ALL OBSERVATIONS FROM A GIVEN STATION.*# THE STANDARD DEVIATIONS 
TO BE ASSOCIATID WITH EACH STATION ARE GIVEN IN COLUMNS 73-79 OF 
THE CARD CONTAINING THE INPUT COORDINATES OF THE STATION. 
PCODE(12)=2 MEANS TO ASSOCIATE A SINGLE STANDARD DEVIATION WITH 
ALL OBSERVATIONS.** THIS NUMBER IS FOUND IN COLS. 21-30 OF THE 
CARD CONTAINING THE TEST DISTANCE (OPTICAL) OR TEST VARIANCE 
(RANGE ) . 

** IN THE CASE OF OPTICAL DBS ERVART IONS , THIS NUMBER IS INTERPRETS 
AS THE STANDARD DEVIAION OF THE DECLINATION AND OF THE RIGHT 
ASCENSION TIMES THE COSINE OF THE DECLINATION, AND THE 
COVARIANCE IS SET TO ZERO. 

PCODE(12)=3 MEANS TO READ ONLY THE DIAGONAL ELEMENTS OF THE 
VARIANCE-COVARIANCE MATRIX (CPGS CORRELATED DATA ONLY) 

13. COMPUTE AMD PRINT CORRELATION MATRIX FOR EACH STATION 
( C£GS CORRELATED DATA ONLY). 
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CODES WI4ICK APPLY TO ORBITAL MODE PROCESSING ONLY 
lA. TREAT COORDINATES OF CENTER' OF MASS AS' UNKNOWNS? (ORBITAL MODE ONLY) 
15. PUNCH UPDATED ORBIT ELEMENTS? (ORBITAL MODE ONLY) ' 


SOLUTION CODES 

16. WRITE NORMALS AND INVERSE DURING SOLUTION PROCESSING? 

0 MEANS PRINT NOTHING 

1 MEANS PRINT PIVOT ELEMENTS 

2 MEANS ALSO PRINT NORMALS AND INVERSE 

3 MEANS ALSO PRINT REARRANGED NROMALS AND INVERSE 

17. PUNCH ADJUSTED STATION XYZ AND VARIANCES FOR INPUT TO BADEKAS' 
DATUM TRANSFORMATION PROGRAM? 

18. PUNCH ADJUSTED STATION POSITIONS? 

19. COMPUTE EIGENVECTORS OF VARIANCE-COVARIANCE MATRIX 

20. COMPUTE CORRELATION COEFFICIENTS 

COMMON/NSTA/NSTA 

INTEGER*? ENDSIG/IHE/,C0NTIN • 

INTEGER*? PC0DE(?0) 

COMMON/PCODES/PCODE 
REAL*8 TITLE! 10) 

3 CONTINUE 
WRITE! 6.6001 } 

6001 FORMAT! lHlr?0( /) ) 

4 READ(5,5001) title, CONTIN 
5001 F0RMAT(9A8,A7,AJ ) 

IFICONTIN.EO.ENDSIG) GO TO 5 
WRITE(6,6012) TITLE 
601? F0RMAT(30X,9A8,A7) 

GO TO 4 

5 CONTINUE 
C 

REA0(5,5050) PCODE 
5050 FORMAT! 8011 ) 

WRITE! 6,6050) PCODE 

6050 FORMAT(////10X, 'PROBLEM COOES lOX , ?01 1 ) 

CALL STAIN 
CALL REAOIN 
CALL AS0360 
CALL FORMRN 
STOP 
END 

DOUBLE PRECISION FUNCTION DPDOT!X,Y,N) 

DOUBLE PRECISION X{N),Y!N) 

DPD0T=0.0 
DO 10 I=1,N 

10 DPDOT=DPOOT+X(I)*Y(I) 

RETURN 

reproducibility of the 
ORIGINAL PAGE IS POOR 
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SUBROUTINE STAIN 
IMPLICIT REAL*8( A-H,0-Z) 

COMMON/PCOD6S/PCOD6 
INTEGER ENDSIG/1HE/,C0NTIN 
COMMON/NSTA/NSTA 
COMMON/STAORD/KORDERI 150) 

INTEGER STANAM*lD5-2 
INTEGER*2 PLUS/1H+/ 

INTEGER={=2 I SGNP , I PH I D» I PH I h, LONGD, LONGM, I SGNL 
COMMON/STALOC/STAUVW { 3» 150 ) , DATPRMI 2» 15 ) ,OATNAM (4, 15 ) , 
1STANAM(5,150)»I0S( 150) 

COMMON/STAPLH/STAPLH( 2t 150 ) 

COMHON/OBSO/OBSD( 150) ,OVOBSD 

MAXSTA=150 

WRITE(6,6000) 

6000 FORMAT! IHl ) 

6001 FORMAT! 1H1,20(/)) 

WRITE(6,6002) 

6002 FORMAT! /////4X,29HDATUMS INVOLVED IN ADJUSTMENT,//) 

C INPUT DATUMS 

10 READ(5,5002) IDD, AE,BE,C0NT1N 

5002 FORMAT! 12, 2F12. 3, 53X,A1 ) • • 

IFICONTIN.EO.ENTSIG) GO TO 30 ' 

DATPRM! 1, IDD)=AE 
DATPRM!2, IOO)=BE 

READ!5,5003) !DATNAM! I, IDD), I=1,A) 

5003 FORMAT (4A8) 

WRITE!6,6003) IOO,!DATNAM| I, IDD ) , I =1, A ) , ! DATPRM !I , IDD) , 1=1 , 21 

6003 FORMAT!6HODATUM,I3,3X,AA8,3HA= ,F10.2,12H METERS B= ,F10-2, 

17H METERS) 

GO TO 10 
C 

30 CONTINUE 

C STATION INPUT 

WRITE!6,6005) ' 

6005 FORMAT! 1H1///40X,29HINPUT COORDINATES OF STATIONS) 

KSTA=0 

35 KSTA=KSTA+1 

READI5,5005)IDDrl0TS,!STANAM! I , KSTA ) , I = 1 , 5 ) , ISGNP, I PHI 0, I PHI M, PHIS 
1,LONGO, LONGM,FLONGS,H,CONTIN 
5005 FORMAT! lA, I2,4A4,A2,A1,2! 213, F8.4), FIO. 2 , 16X , A1 ) 
IF!CONTIN.EO'.ENOSIG) GO TO 50 
PHI=ANRADD! I SGNP, I PHI D, I PHI M, PHI S ) 

1SGNL=PLUS 

FLONG=ANRADD! I SGNL , LONGD, LONGM, FLONGS ) 

KORDER! KSTA)=IDD 
I0S!KSTA )=IDTS 
STAPLH!1,KSTA)=PHI 
STAPLH!2,KSTA)=FL0NG 

CALL UVV'!D!OATPRM! 1 ,! OTS ), DATPRM! 2, I DTS ) , PHI , FLONG ,H, STAUVW 1 1 , KSTA) 
ItSTAUVW! 2rKSTA J , STAUVW! 3, KSTA ) ) 

WRITE! 6, 6006 HDD, ! STANAM! I , KST A ) , I =1 , 5 ) , IDTS , ! DATNAM ! I , I OTS ) ,1 = 1,4 
1), ISGNP, IPH ID, I PHIMtPHIS, 1 SGNL , LONGD, LONGM, FLONGS , H 

6006 FORMAT! iHo, 14, BX, 4A4, A2, lOX, 5HDATUM, I4,4X,4A8/10X,20HGE00ETIC COOR 
1D1NATES.2! 6X, A1,2I3,F8.4) ,F12.4) 

WRITE !6,6007) ! STAUVW! I , KST A ) , I =1 , 3 ) 

6007 FORMAT! lOX, 21HCARTESIAN CliORD I N ATES, 3F16. 3 ) 

GO TO 35 

50 CONTINUE 
NSTA=KSTA-1 
NSTAUN=3«NSTA 
RETURN 
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SUBROUTINE READIN 
IMPLICIT REAL'-!=8( A-HtO-Z) 

INTEGER>:--2 PC00E(20) 
common/ PCODES/PCOOF 

INTEGER*4 ENOS I G/ IHE/ 1 CQNT I N, DELCOD ( 2 ) / IH » IH*/ , ECOOE 
INTEGER«2 PLUS/JH+/ 

1NTEGER=!=2 ISGNP, IPHIO, IPHIM, LONGO,LONGM, ISGNL 

INTEGER-2 I DI 50 ) » KEY ( 50 ) , I HR ( 50 ) » Ml N ( 50 ) , I DAY ( 50 ) , I YR ( 50 ) » IGHHf 50) 
i,IGHM(50)» ISGND(50), I0eCD(50), lnECM(50)» IOAT(50,ll) 

COMMON/DEDITC/ GHA ( 50 ) , DEC( 50 ) , Ut 3, 50 ) , S ( 3) r D( 50) , 

1SDC(3,50).EVSUM,STAXYZ(3,50),G01, 

2TD,KSTATE(50) » IPASS(50),NSTE,NSUSED,EC0DE 
COMMON/NSTA/NSTA 
COMMON/STAORD/KORDERC 150) 

INTEGER STANAM, I0S=i=2 

C0MM0N/STAL0C/STAUVW(3»150),0ATPRM{2»15),DATNAM(4,155 , 
iSTANAM(5t 150), IDS ( 150) 

DIMENSION M0NTH(50) 

EQUIVALENCE! ID(1 ), IDAT( 1 , 1 ) ) , ( KEY ( 1 ) , IDAT ( 1 , 2 ) ) , C IHR { 1 ) , lOAT { 1 , 3 ) ) 
1, (MIM( 1 ) , IDAT( 1,4) ) , (I DAY! 1 ), lOAT! 1,5) ),! IYR( 1 ), I DAT 1 1,6) ) , I IGHH! 1 
2),IDAT!l,7)),!IGHM!l),IDAT!l,8)),!ISG.NiD!l),IDATIl,9) ) , ! IDECDI 1 ) , ID 
?AT!1,10) ),! IDECM 1 ) , I DAT ( 1 , 1 1 ) ) 

DIMENSION DAT(50,6), DECS! 50 ) , VARGHA ! 50 ) , VAR06C ! 50) , COVGHD! 5 

10),GHAS! 50),SEC(50) 

EQUIVALENCE! SEC! 1 ) , DAT ( 1 . 1 ) ) , ! GHAS ! 1 ) , DAT! 1 , 2 ) ) , ! DECS ! 1 ) , DAT ( 1 , 3) ) 
X, (VARGHA! 1 ) ,DATM,4) ), ! VARDEC! 1 ),DAT! 1,5) ) , ! CQVGHD! 1 ) , OAT ! 1, 6) ) 
COMMON/OBSD/OBSD! 150),0V0BSD 
3 MAX5TE=50 

PI=3. 141592653589793D0 
RH0=180.D0/PI 
SPR=RHO^=36O0.DO 
PI2=2.D0*PI 
WPWSP=0.D0 
C 

REA0!5,5004) TO 
WR1TE!6,6004) TO 
5004 FORMAT! F20«2,F10. 2) 

6004 FORMAT! //20X, ‘TEST DISTANCE =',F20.2,' SECONDS OF ARC') 

WRITE!3) TO 
C START DATA INPUT 

IEVENT=0 
KEVENT=o 
EPR=0.0 
IS=0 


ENTER HERE FOR A NEW OBSERVATION 

200 IS=IS+1 
C READ DATA CARD 
211 CONTINUE 

READ!5, 5000,END=901 ). ID! IS), lYR! IS) ,MONTH! IS) , IDAY! IS) ,IHR! IS) / 
1M1N!IS), 

2SEC! IS ) , IGHH! IS), IGHM! IS), GHAS! IS) , ISGND! IS) , 1 DECD! IS) , IDECM (IS) , 
3DECS! IS), VARGHA ! IS ), VARDEC! IS) ,COVGHD! IS) ,CONTIN 

5000 FORMAT! IX, 1 4, 1 X , 5 1 2 , F3. 0, 5X , 2 1 2 , F7. 4 , 5X , A 1 , 2 1 2 , F7.4, 1 X , 2F6. 3 , 
CF7.3,7X,A1) 

C 

IF!CONTIN.EO.ENOSIG) GU TO 250 
DDT=DFLOAT!MJD! IDAY! IS), MONTH! IS),IYR!IS))) 

DDT = DDT + lOFLOAT! ! IHR! IS)*60 + MIN(IS))*60) + SEC I I S ) ) /86400. 

IF! IS.LE.l ) GO TO 212 
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FOR THIS EVENT. BEGIN PROCESSING 


C THIS TEST SHOULD BE TRUE ONLY FOR THE FIRST CARO OF THE FIRST EVENT. 
C 

CHECK FOR END OF EVENT, ALLOWING 0.5 MILLISECOND DISCREPANCY. 
IF(OA8S(DDT-EPR).GT.0.58D-8) GO TO 250 

C 

C ENTER HERE TO BEGIN A NEW EVENT 

C THE FIRST ENTRY OF THE EVENT SHOULD ALWAYS BE MADE WIHT IS=1 

212 CONTINUE 
IDD=ID(IS) 

KSTA=KSTAIO( IDD) 

IF(KSTA.GT.O) GO TO 220 

WRITE (6, 6042) I D( IS ) , IHR (I S ) , MINI IS),SEC( IS),IDAY(IS) , MONTH! IS) , 
IIYR(IS) 

GO TO 211 
220 CONTINUE 

reproducibility of the 

GO TO 200 ORIGINAL PAGE IS POOR 

c 

C END OF INPUT 
250 CONTINUE 

NSTE = IS - 1 
lEVENT = lEVENT + 1 
6042 FORMAT (5 X, 'STATION' number NOT 
1F8.4,3X, I 3, A3, 12, •• OBSERVATION 
00 270 IS =1,NS1E 
ISGNL = PLUS 

GHAnS)=ANRADD( I SGML , I GHH (I S ) , I GHM{ I S ) , GHAS ( I S ) ) ^-'15 .DO 
DEC! 1S) = ANRADD( ISGND! I S ) , I DECO ( I S ) , lOECM! IS),DECS!IS) ) 

270 CONTINUE 

WRITEC3) lEVENT, NSTE, EPR, 

11! IDAT! IS, J ), J=l,ll ),MONTH( IS), IDATl IS, J), J=l,6) ,GHA ! IS) ,DEC! IS) 
2KSTATE! IS), IS=1,NSTE),C0NTIN 
C TEST FOR END OF INPUT 

IFICONTIN.EO.ENDSIG) GO TO 700 
C PREPARE FOR NEXT EVENT 
DO 610 1=1,6 

610 DAT!1,I)=DAT(NS1E+1,I) 

MONTH! 1 ) = MONTH! NSTE+1 ) 

DO 611 1=1, 11 

611 IDATd, I ) = I0AT(NSTE+1, I ) 

C RETURN TO START A NEW EVENT 
IS=1 

GO TO 212 


FOUND IN INPUT 
IGNORED* ) 


LIST* ,15, 3X, 213, 


700 RETURN 


C ERROR EXITS 
901 CONTINUE 
C ENTER HERE IF END 
CONTIN= ENDSIG 
GO TO 250 
END 


SIGNAL CARD IS MISSING FROM INPUT DATA DECK 
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SUBROUTINE ASD360 

C S/360 VERSION OF ASO FOR OPTICAL SATELLITE DIRECTIONS 
IMPLICIT REAL=!=8(A“H,0-ZJ 
INTEGER*2 PC0DE(20J 
COMMON/PCODES/PCOOE 


1NTEGER=!'4 ENDSIG/lHE/,CONTIN,DELCOO< 2 )/lH »1H*/,EC0DE 
1NTEGER*2 PLUS/IH+/ 

INTEGER*2 ISGNP, IPHIO, IPHl M, LONGO, LONGM, ISGNL 

INTEGER*2 I D( 50 ) » KEY { 50 ) , IHR ( 50 ) , MIN ( 50 ) , I DAY( 50) , I YR ( 50) , IGHH ( 50l 
1« IGHM(50)» ISGNDf 50), IDECDI 50), IOECM( 50), IDAT(50,11) 

COMMON/DEDITC/ GHA ( 50 ) , DEC ( 50 ) , U( 3 , 50 ) ,S ( 3 ) , 0( 501 , 

1SDC(3, 50) ,EVSUM, STAX YZ I 3, 50 ) , GQI , 

2TD,KSTATE( 50) , I PASS ( 50 ) ,NS1 E,NSUSED, ECODE 
eOMMON/NSTA/NSTA 
INTEGER STANAM,IDS-2 
DIMENSION MONTH! 50) 

COMMON/STALOC/Sl AUVWI 3, 150) , DATPRMl 2, 15 ) , DATNAM! A-, 1 5 ) , 

1STANAM(5, 150), IDS! 150) 

COMMON/STAORO/KORDERI 150) 

EQUIVALENCE! ID! 1 ), I DAT! 1, 1 ) ) , (KEY! 1 ), IDAT! 1,2) ) , ! IHR f 1 ) , I DAT 1 1 , 3 > ) 
l,(MlN{n,IDATn,A)),!IDAYm,IDAT(l,5)),! lYR! 1 ),IDAT! 1 , 6 ) ) ,( IGHHIl 
1 ), IDAT! 1,7) ), ! IGHM! 1 ) , I DAT ! 1 , 8 ) ) , ! ISGND! 1) , IDAT! 1,9) ), ( IDECDC 1) , ID 
3AT!1,10) ),! lOErXDflOATl ],11)) 

DIMENSION DAT!50,6), DECS ! 50 ) , VARGHA ! 50 ) , VAROEC { 50) ,C0VGHD( 5 

10),GHAS( 50) ,SEC!50) 

> VALENCE! SEC! 1 ) ,OAT! 1,1 ) ) , ( GHAS 1 1 ) , DAT ! 1 , 2 ) ) , ( DECS ( 1 ) , DAT ( 1 , 3 ) I 
1.. V'RGHA! 1 ) ,DAT! 1,A) ) , ( VARDEC ! 1 ) , DAT ! 1 , 5 ) ) ,!COVGHD! 1) , OAT! 1,6) ) 
i>i',cNSION SSDCf3,4),AW(3),A! 3,3),AT!3,3) ,DDN!3,3) ,W!3,3) , 

I N0BSTA(150),A1(3,3),A1T!3,3), DDK! 3 ) , D1 ( 3 , 3) , AH ! 3 , 3 ) , 

2BT(3,3),TEMP1!3, 3),TEMP2(3,3),BN!3,3,6),TEMP3(3) ,DN(3,3,50) , 
3DK!3,50),AK1!3),TA(?) 

COMMOM/WPW/WPW,XPU, IDEGF,NFSTA 
DIMENSION VPVSTA!150) 

MAXSTE=50 

PI =3.141592653589797300 

PI2=2.*PI 

RPD=180.0/PI 

SPR = { 180.’!=3600. )/PI 


WPWSP=0.0 

REWIND 2 
REWIND '3 
WRITE!6,4397) 

4397 FORMAT! IHl) 

READI3) TD 
WR1TE!6,6004) TD 
6004 FORMAT! //20X, 'TEST 




DISTANCE =',F20.2, • 


•SECONDS OF ARC) 


START DATA INPUT 
KEVENT=0 
EPR=0.0 

DO 70 KSTA=1,NSTA 
MOBSTA!KSTA)=0 
VPVSTA!KSTA)=0.0 
DO 70 1=1,3 
DK! I,KSTA)=0.0 
■ DO 70 J=l,3 
DN! 1, J,KSTA)=0.0 
70 CONTINUE 
DO 80 1=1,3 
DO 80 J=l,3 
A1!I,J)=0.D0 
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A1T(I,J)=C.D0 
80 CONTINUE 

00 31A J=l,3 
A1(J,J)=-1. 

31A A1T(J.J)=-1. 

210 CONTINUE 

READI3) IEVENT,NSTE,EPR, 

1(( lOATI IS, = ID.MONTHI I S ) , I DAT ( IS, J ) , J=1 ,6 ) ,6HA IIS),OECnS>, 
2KSTATE( IS), IS=1,NSTEI,C0N11N 
DO 272 IS=1,NSTE 
KSTA=KSTATE( IS) 

DO 272 M=l,3 

272 STAXYZ{M,I5)=STAUVW(H,KSTA) 

WRITE(6,6008) lEVENT 
6008 FORMAT!/ IX , • EVENT • , 16 ) 

C 

CALL DEDIT 
C 

DO 280 1S=1,NSTE 

WRITE (6, 6010) I0( IS), IHR( IS),MIN( IS), SEC { IS), IDAY( IS) ,MONTH( IS) , 
HYR( IS) , IGHH( IS ), IGHM! IS) ,GHAS( IS),IS6ND( IS) ,IDECD( I S ) , I DECM ( I S ) , 
2DECS( IS), VARGHAI IS ), VARDECI IS ) ,COVGHD( IS ) ,0( IS) ,DELCOO( IPASSI IS) ) 

6010 FORMAT! 1 7 , 2 1 3 , F9„ 5, 3X , 1 3, A3, 12, 2 1 3, F8. A, 3X , A1 , 12 , 1 3 ,F8. 4, 5X, 3F6.2 , 
1F10.0,2X,A1) 

280 CONTINUE 

IF(ECOOE.GT.l) GO TO 630 
IFIPCODEdD) 290,630,610 

610 IF! PCODE! 11 )-3) 611,612,611 

611 WRITE(6,6024) S 

602A FORMAT! • SATELLITE POS I T I ON ' , 3F15. 3 ) 

IF(PC0DE!ll)-2) 612,630,612 

612 IDTS= IDS! KSTATE! 1 ) ) 

CALL UVWTG2(S,DATPRM! 1, I01S),PHI ,FLAM,H) 

PHI = PHI=(=RPD 
FLAM=FLAM*RPD 
WRITE! 6,6023 ) PHI, FLAM, H 

6023 FORMAT! • GEOD. COORD. OF S ATELL I TE • , 2F1A. 6, F14, 1 ) 

630 CONTINUE 

WRITE!6,6012)GOI 
6012 F0RMAT!10X,«G0I=',F10.5) 

IF! ECODE.GT.l ) GO TO 290 
IFINSUSED.EO.O) GO TO 290 
RMSMC=DSORT! EVSUM/DFLOAT(NSUSEO) ) 

WRITE!6,6011) RMSMC 

6011 FORMAT! 1H+,27X, «RMS MISCLDSURE IN ME! ERS= » , FIO. 1 ) 

GO TO 300 

290 WRITE! 6,6015) ECOOE 

6015 FORMAT! 1H+,27X, 'ENTIRE EVENT DELETED, K0DE=',I4) 

GO TO 600 

C 

C SET UP GENERALIZED LEAST SQUARES EQUATIONS FOR THIS EVEN AND COMPUTE 
C CONTRIBUTIONS TO THE NORMAL EQUATIONS. 

300 CONTINUE 

IF! ECODE.GT. I ) GO TO 600 
KEVENT=KEVENT+1 

DO 310 1=1,3 . ■ 

DDK! 1 )=0.0 
DO 310 J = l, 3 

DDN!I,J)=0.0 • ‘ 

310 CONTINUE 

C 

JS=0 


„EPRODUCKlLn'V OF T^ 
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DO 390 IS=ltNSTE 

IF( IPASS( ISJ.GT. 1 ) GO TO 390 


C 


JS=JS+1 

JS IS THE COUNTER FOR NOW-DELETED STATIONS IM THE EVENT 
GHA ( IS)=PI2*-GHA (IS) 


RSOCSD=SDC( 1, rS)*^^2+SDC(2» IS)=«=«2 
RSO=RSCCSD+SDC( ?,IS)*=i=2 
RC0=0SORT(RS0CSD) 
RANGE=DSORT(RSO) 

SG=DSIN(GHA (IS)) 

CG=DCOS(GHA (IS)) 

SD=DSIN(DEC( IS) ) 

CD=DCOS(OEC( IS) ) 
A{1.1)*=SD^CG*RANGE 
A( 1,2) = SG=!:CD^RANGE 
. A(3,3)=-SD*RANGE 
A(2,1) = SD=!'SG*RANGE 





A( 2, 2 ) =-CG^CD=i=RANGE 

A(2,3)=-CD*SG*RANGE 

A(3,1)=-CD=!=RANGE 

A(3,2)=0.0 

A( l,3)=-CD*CG=f=RANGE 

AW(1)= SOC( 1»IS)-RANGE*CG*CD 

AW(2)= S0C(2»IS)-RANGE=J=SG=«'CD 

AW(3)= S0C(3» IS)-RANGE=i=SD 


C 

COMPUTE WEIGHTS 

VARGHA( IS)=fVARGHA( IS)/SPR)**2 
VARDEC( rS) = (VARDEC( IS)/SPR )=5==i=2 
COVGHD( IS )=COVGHD( IS )/SPR*=i'2 
OET=VARGHA( IS )*VARDEC ( I S ) -COVGHO( IS)#*2 
W(1,1)=VARDEC( IS)/DET 
W(1,2)=-C0VGHD( IS)/DET 


W(l,3)=0. 
W(2fl)=W(l»2) 
W(2,2)=VARGHA{ IS)/DET 
W(2f3)=0. 

H(3»l )=0. 

W(3t2)=0. 


W(3,3)=0. 

C 

KSTA=KSTATE( IS) 

C ELIMINATE DELETED STATIONS FROM THE LIST OF STATIONS INVOLVED IN 
C THE EVENT 

KSTATE( JS)=KSTA1E( IS) 

C 

CALL VERS0L(A,BT,3,3) 

DO 9A0 1=1,3 
DO 940 J=l,3 
940 A( I, J)=BT( 1, J) 

DO 821 1=1,3 
DO 821 J=l,3 
821 AT(J,1 )=A(1,J) 

CALL DGMPRD(AT,W,TEMP1,3,3,3) 

CALL DGMPRD(TEMP1,A,TEMP2,3,3,3) 

CALL 0GMPRD(A1T,TEMP2,BM( 1,1,JS),3,3,3) 

CALL DGMPRD(TEMP2,AW,TEMP3,3,3, 1 ) 

DO 915 .1 = 1,3 

DO 915 J=1 ,3 . 

915 BN( I, J, JS )=-BN( I, J,JS) 

DO 916 1=1,3 

DDK( I )=DOK( I )+TEHP3( I ) 
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6910 FORMATdH ,4D?2.1A) 

DO 916 J=l,3 

916 DDN( E, J)=DDN( 1 . J)+TEMP2( I ♦ J 1 
00 330 1=1,3 
DO 325 J=l,3 
TERM=0.0 
DO 320 n=l,3 
00 320 JJ=1,3 

320 TERM=iTERM+Al( II,n*1EMP2< 1 1 , J J l=J=A i IJ J, J ) 

DN( I, J,KSTA)=DN( I , J,KSTA )+TERM 
325 CONTINUE 
TERM=0.0 
DO 328 11=1,3 
00 328 JJ=1,3 

328 TERM=TERM+A1( 1 1 , n=«=TEMP2 ( JJ) 

DK{ I,KSTA)=DK( I,KSTA)-TERM 
330 CONTINUE 

CALL DGMPR0(S0C( 1, IS), TEMP2,AK1, 1,3,3) 

CALL DGMPRD ( AK 1, SDC ( 1, IS ),VPVT0, 1,3,1) 

WRITE<6,938) VPVTO 

938 FORMATdH , "NEW VPV‘,D20.J2) 

KN0=K0RDER(KSTA) 

VPVSTAtKSTA ) = VPVSTA ( KSTA ) +VPVTO 
NOBSTA ( KSTA ) =NnBSTA{ KSTA )+2 

390 CONTINUE 

FORK REDUCED NORMAL EQUATIONS 
INVERT DDN 

CALL VERSOLIDDN, 8T,3,3) 

CALL DGMPRD{DDK,BT,TA,1,3,3) 

CALL DGMPRD(TA, DDK, TB, 1,3,1) 

WRITE(6,939) TB 

939 FORMATdH , ‘WPW CONTRIBUTION FROM SATELLITE POS ITION » , 020, 12) 
WPWSP=WPWSP+TB 

NSUSEO=JS 

WRITE! 2) NSUSED,BT,DOK, ( ( ( BN (1 , J , JS ) , I =1 , 3 ) ,J=1,3) ,KSTATE( JS) , 
1JS=1,NSUSED),C0NTIN 
600 CONTINUE 

TEST FOR END OF INPUT 
IF(CONTIN.EO.ENOSIG) GO TO 700 
GO TO 210 


700 CONTINUE 

CHECK TO SEE IF END SIGNAL HAS BEEN WRITTEN ON DATA SET FT02 
1F( ECOOE. EO, 1 ) GO TO 710 
BACKSPACE 2 

C READ AND REWRITE LAST RECORD FROM LAST GOOD EVENT 

READ! 2) MSUSED,BT,DDK, ( ( < BN ( I , J , JS ) , 1 =1 , 3 ) , J = I ,3 ) , KST ATE !JS) , 
1JS = 1, NSUSEO) 

BACKSPACE 2 

WRITE! 2) NSUSED, BT,ODK, ! ! f BN ! I , J , JS ) , I =1, 3 ) , J=1 , 3 ) , KSTATE C JS ) , 
ljS=l,NSUSED),CaMTIN 
710 CONTINUE 

WRITE(2) !!! DN! I , J , KSTA ) , I = 1 , 3 ) , OK ! J , KST A 1 , J = ] , 3 ) , 

1KSTA=1,NSTA) 

WPW=0.0 

N0BS=0 


oeiginal page is 
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WR1TE(6,6019) 

6019 FORMAK /), lOX, 'ANALYSIS OF MISCLOSURES BY STATION'// 

ITIO, 'STATION', T?Or 'NUMBER OF 08SERVA1 IONS ', T50 ,» RMS MISCLOSURE'J 

DO 750 KSTA=1,NSTA 
N0BS=N0BS+NG8STA(KS7A ) 

WPW=WPW+VPVSTA(KSTA) 

RMSMC=0.0 

IF(NOBSTA{KSTA).GT,0) RMSMC=DSORT ( VPVSTA ( KSTA ) /DFLOAT (NOBSTA ( KSTAI 
1 ) ) 

WRITE (6, 6020) KORDER ( KSTA ) , NOBSTA ( KSTA ) , RMSMC 

6020 FORMATITIO, I7,T?5,I7,T50,F14.2) 

750 CONTINUE 

IDEGF=N0BS~3-KEVENT 

RMSMC = DSORT( WPW/OFLOAT( IDEGF) ) 

WRITE! 6,6021 ) NOBS, KBVENT , IDEGF, WPW,RMSMC 

6021 F0RMAT(////10X, 'TOTAL NUMBER OF GOOD OBSERVATIONS* ,T60, 18// 

IlOX, 'TOTAL NUMBER OF GOOD EVENTS ' ,T60, 18// 

IlOX, 'CORRESPONDING DEGREES OF FREEDOM ' , T60 , 1 8// 

IlOX, 'TOTAL SUM OF SQUARES OF MI SCLOSURES • , T60, FI 1 . 2// 

IlOX, 'CORRESPONDING STANDARD DEVIATION OF UNIT WE IGHT ' ,T60, FI 1 . 2 ) 
WPW=WPW-WPWSP 
WRITE! 6, 6022)WPW 

6022 FORMAT! 1H0,SX, 'kPW INCLUDING CONTRIBUTION FROM SATELLITE POSITION' 
1/15X, '! 1. E.,VPV+UX) ',T60,F11.2) 

RETURN 

END 
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SUBROUTINE FORMRN 
implicit REAL-a(A-H,0-Z) 

common/nsta/nsta 

INTEGER=^2 PC0DE(20) 

COMMON/PCODES/PCOOE 
COMMON/WPW/WPW, XPU, IDEGF, IFSTA 

DIMENSION DDN{3, 3), DDK( 3 ) , LI { 3 ) , L2 ( 3 ) , BNDDNI ( 3,03 ) , TN ( 3, 3) , 

.JTK(3) ,DDL(3),DN(3,3) 

INTEGER^=2 L,LS0LVE 
INTEGER C0NTIN,ENDSIG/1HE/ 

COMMON/STAORD/KORDER( 150) • 

COMHON/NORMEQ/L SOLVE 

DIMENSION REDN(3,3,1275),U(3,50),L( 1275) 

DIMENSION 5N{3,?,50),LG(50) 

C FORM REDUCED NORMAL EQUATIONS FOR UP TO 50 STATIONS 
DIMENSION KSTATE(50) - 

L0C(K)=(K#(K+l))/2 
MAXSTA=50 

IF(NSTA.GT.MAXSTA) GO TO 901 

THE REDUCED NORMAL EQUATIONS ARE STORED AS 3 X 3 BLOCKS IN THE ARRAY REDN. 
ONLY THE UPPER TRIANGULAR PART OF THE REDUCED NORMAL EQUATIONS IS STORED. 
THE BLOCKS OF THE REDUCED NORMAL EQUATIONS ARE NUMBERED 
ACCORDING TO THE FOLLOWING SCHEME: 

12 4 7 11 • 

3 5 8 12 

6 9 13 

10 14 

15 ET CETERA 

L(1275) IS THE GUIDE MATRIX 
L=1 SIGNIFIES A NON ZERO BLOCK 

L=0 SIGNIFIES A ZERO BLOCK 

IB=LOC(NSTA) 

DO 100 JB=1,IB 
DO 99 1=1,3 
DO 99 J=l,3 
99 REON( I, J, JB)=0.0 
100 L(JB)=0 

BACKSPACE 2 

READt2) n ( 8N( I, JtKSTA), 1=1,3),U(J,KSTA),J=1,3) , 

XKSTA=1 ,MSTA ) 

REWIND 2 

STASH DIAGONAL BLOCKS 
00 110 KSTA=1,NSTA 
IB =L0C(KSTA) 

DO 108 1=1,3 
DO 108 J=l, 3 • 

108 REDN( 1,J,IB)=PN( I,J,KSTA) 

110 CONTINUE 

FOEGF=IDEGF 

IF(PC0DE(9).EQ.l ) WR I TE ( 7 , 701 0 ) FDEGF,WPW 
010 FORMAT) 16X,2F16. 6) 

READ BLOCKS FROM EACH EVENT AND REDUCE NORMAL EQUATIONS 

150 READ(2) NST6,DDN,DDK, ( ( {BN( !, J, IS ), 1=1,3) , J=l,3) , 

IKSTATE) IS),]S=1,MST6),C0NTIN 
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00 100 IS=1,NSTE 
ISTA=KSTATE( IS) 

IB=1STA 

CALL DGMPRD(BNn,l,IS),00N,BNDDNl,3.3.3) 

CALL DGMPRD(BNDDNI,DDK»TK,3»3,1) - 

00 155 1 = 1.3 ■ 

155 un, ISTA)=U( I, ISTA )-TK{ I) A! i ^ 

00 100 JS=1,NSTE 
JSTA=KSTATE( JS) 

JB=JSTA 

SKIP IF ( ISTA.GT. JSTA), SINCE ONLY THE UPPER TRIANGULAR PART OF THE 
REDUCED NORMAL EQUATIONS IS BEING COMPUTED AND SAVED. 

IF( ISTA.GT.JSTA) GO TO 180 

<IB»JB) GIVES THE ROW AND COLUMN NUMBER OF THE BLOCK IN THE REDUCED 
NORMAL EQUATIONS CURRENTLY BEING PROCESSED. 


SET INDICATOR 
NB=L0C(JB-1) 

NB=IB+NB 
UNB)=L(NB)+1 
PERFORM REDUCTION 
00 156 1=1,3 
DO 156 J=l,3 
156 ONU, n-BM( I, J, JS) 

CALL DGMPRDf 6NDDN I , DN, TN, 3 ,3 , 3 ) 

5910 FORMATdH ,3020.12) 

DO 130 1=1,3 
DO 130 J=l,3 

130 REDN( I, J,NB)=REDNC I,J,NB)-TN( I,J) 

180 CONTINUE 

IF END OF DATA, GO OUT OF LOOP 
IF(CONTIN.EQ.ENOSIG) GO TO 400 
IF NOT, RETURN TO PROCESS ANOTHER EVENT 
GO TO 150 

ENTER HERE WHEN ALL EVENTS HAVE BEEN ’ PROCESSED. 
400 CONTINUE 

SIMULATE KRAKIWSKI'S GUIDE MATRIX 
IF(PC0DE(6).NE.l) GO TO 441 

WRITE(6,6001) 

>001 FORMAT! 1H1,10( /) ,20X, 'GUIDE MATRIX • ) 

DO 440 ISTA=1,MSTA 
IB=0 

LG(1)=1000 

DO 435 JSTA=ISTA,NSTA 
JB=LOC( JSTA-1 )+ISTA 
IF(L{ JB). EO.O) GO TO 435 
IB=IB+1 

LG(IB)=KORDER( JSTA) 

435 CONTINUE 

IB=IB+1 

IF(IB.GT.l) LG(IB)=999 

439 WRITE(6,6002) KORDER ( I ST A ) , { LG( I ) , I = 1 , I B ) 
i002 FORMAT! 20X, 15, 5X, 1815,200! /30X, 1RI5)) 

440 CONTINUE 

441 CONTINUE 


PRINT NORMALS IN ASD FORMAT, AND PUNCH 


DESIRED. 
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WR1TE(6,6003) 

6003 FORMAK IHl//' NORMAL EOUATIONS (SEE GUIDE MATRIX)'//) 

DO A50 1STA=IvNSTA 

00 AA2 1=1,3 
4A2 DDK 1 ) = -!)( I, ISTA) 

IB=0 

JB=L0C( ISTA) 

IF(L( J0).GT.O) 18 = 1 
C PUNCH normals 

IF(PC0DE(9).NE. 1 ) GO TO A43 
WRITE(7,7001 ) KCIRDER( ISTA) 

7001 FORMAT! IAI5) 

WRITE(7,7006) DDL 
7006 FORMAT! 3025. 16) 

WRITE(7,7008) ( (REDNU , J, JB ) , J=l, 3) , 1 = 1,3) 

7008 F0RMAT(3D25.16/3D25. 16/3025.16) 

C 

443 CONTINUE 

• C PRINT DIAGONAL BLOCK 

IF(PC0DE(7). NE. 1 ) GO TO 444 
WRITE! 6,6004) KORDER(ISTA) 

6004 FORMAT!// 15) • 

WRITE!6,6006) DDL 

6006 F0RMAT!3O25.16,//) 

WRITE!6,6008) ( ! REDN ( I , J , JB ) , J=1 , 3 ) , 1 = 1 , 3) 

444 CONTINUE 

C PRINT OFF-DIAGONAL. BLOCKS ' ' 

KSTA=ISTA+1 

IF! ISTA.EO.NSTA ) GO TO 448 
DO 445 JSTA=KSTA,NSTA 
JB=L0C! JSTA-l )+ISTA 
IF!L! JB).EO.O) GO TO 445 
IB=IB+1 

IF(PC0DE!9).N£.l ) GO TO 7445 
WRITE(7,7001) KDRDER(JSTA) 

WRITE(7,7008) ( ( REDN ( I , J , JB ) , J=1 , 3 ) , I =1 , 3 ) 

7445 CONTINUE 

IF(PC0DE!7).NE.l ) GO TO 445 
WRITE!6,6004) KORDER(JSTA) 

WRITE!6,6008) { ! REDN ( I , J, JB ) , J = 1 , 3 ) , I = 1 , 3 ) 

445 CONTINUE 
448 1=1000 

IF! I8.GT.0) 1=999 

IF(PC0DE!7).E0,1 ) WR I TE ( 6, 6004 ) I 
IF(PC0DE!9).EC.] ) WRITE(7,7001 ) I 
450 CONTINUE 

IF!PC0DE!8) .NE. ] ) GO TO 478 
WRITE!6,60I0) 

6010 FORMAT! 10( /) ,20X, 'OBSERVAl IONS ON EACH LINE') 

IB=MSTA-1 

DO 475 ISTA=1,IB 

KSTA=ISTA+1 

DO 475 JSTA=KSTA,NSTA 

WRITE! 6,6011 ) KORDER ! I STA ) , KOROER ( JSTA ) , L ! LOC ! JSTA-1 ) + I STA ) 

6011 FORMAT! 8110) 

475 CONTINUE 
478 CONTINUE 

RETURN 

901 CONTINUE 

WRITE! 6,9001 ) MAXSTA,NSTA 

9001 FORMAT!' FORMRN IS PRESENTLY DIMENSIONED TO HANDLE ONLY', 15, 

1' UNKNOWN STATIONS. »/20X, » this PROBLEM HAS',15,' UNKNOWN STATI 
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2I0NS. '/lOX, 'EXECUTION IS TERMINATED BY PROGRAM.*) 

! 


6008 F0RMAT(3025.16) 
STOP 
END 
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SUBROUTINE DEDIT 
implicit REAL=i=R(A-H,0-Z) 

COMMON/DEDITC/ ALES t50 » » DEC ( 50 ) » U( 3 » 50 ) ,S { 3 1 ,D(50) , 

lSnC(3,50),SUM,STAXYZ(3t50),G0I, 

2Tf),KSTATE(50 , I PASS ( 50 ) , NSTE, NSUSEO, KODE 
EDIT DATA BASED ON PRELIMINARY STATION POSITIONS AND DELETE BAD 
OBSERVATIONS AMD BAD EVENTS, BASED DM THE DISTANCE CRITERION TD 
this SUBROUTINE IS DIMENSION FOR A MAXIMUM OF MAXSTE=50 STATIONS 
PARTICIPATING IN ANY ONE EVENT, ALL AFFECTED ARRAYS ARE IN 
COMMON BLOCK /DEDITC/. 

THE NUMBER OF STATIONS PARTICIPATING IN THE EVENT IS NSTE, 

THE NUMBER OF STATIONS NOT DELETED IS NSUSED. 


COMMON/STALOC/STAUVWI 3, 150) 

DIMENSION 0(3,3) ,RHS(3),0I(3t3).VI(3) 


PI=3, 141592653589793D0 
RH0=180. DO/PI 
SPR=RHO=!=3600.DO 
TPI=2,D0=s=PI 
MAXSTE=50 
INITIALIZE 
K0DE=1 

DO 110 IS=1,NSTE 
110 IPASS( IS)=1 

IPASS=1 MEANS 


THIS DIRECTION OK 




IPASS=2 MEANS THIS DIRECTION DELETED FROM EVENT 


FORM UNIT VECTORS FOR ALL DIRECTIONS IN THIS EVENT 
DO 125 IS=1,NSTE 
STS=TPI“ALFS( IS) 

CA=DCOS(STS) 

SA=DSIN(STS) 

CD=DCOS(OEC( IS) ) 

SO=DSIN(OEC( IS) ) 

U( 1, IS)=CA*CD 
U(2, IS)=SA*CD 
U(3,IS)=SD 
125 CONTINUE 

INITIALIZE ARRAYS FOR THIS ITERATION 
130 CONTINUE 
NSUSE0=0 
DO 140 1=1,3 
RHS( I )=0.0 
S(I)=0.0 
DO 140 J=l,3 
0(I,J)=0.0 
140 CONTINUE 

ACCUMULATE EQUATIONS 
DO 190 IS=1,NSTE 
IF( IPASSC IS).EO. 2) GO TO 190 
NSUSED=NSUSED+1 
DO 170 1=1,3 
DO 169‘ J = l,3 

169 0I(I,J)=U(I,IS)*U(J,IS) 

170 01 ( I , I )=0I( I, I )-] .0 
DO 175 1=1,3 

DO 175 J=l,3 

0( 1, J)=0( I,J)+QI { I,J) 
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RHS( I )=RHS( n+OI ( 1, J)«STAXYZ( J, IS) 
175 CONTINUE 
190 CONTINUE 

TEST FOR DELETION OF WHOLE EVENT 
IF(NSUSED,LT.2) GO TO 420 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


INVERT AND SOLVE 

THE SATELLITE POSITION S IS SELECTED IN SUCH A WAY THAT THE SUM OF 
THE SOIJARES OF THE DISTANCES FROM S OF THE NON-OELETED RAYS IS MINIMIZED. 
DET=1.0 

CALLOMIMVIOf 3,DETtOI( lTl)tOI(1.2)) 

GOI=OABS(DET/DFLOAT(NSUSED)) 

IFIGOI.LT. 1.00-4) GO TO 430 
CALL 0GMPRD(0,RHS,S,3,3,1) 

COMPUTE DISTANCE FROM 5 FOR EACH RAY 
ISMAX=0 
DMAX=0.0 
SUM=0.0 

DO 280 IS=1»NSTE 
DO 270 r=l»3 
DO 269 J=l,3 

269 OintJ)=UU,IS)*U(J,IS) 

OI(I,I)=OI( I, D-l.O 
Vim=S( I )-STAXYZ( I, IS) 

270 CONTINUE 
0DI=DPD0T(VI,U(1, IS),3) 

DDI=DABS(DDI ) 

01=0.0 f 

DO 275 1=1 » 3 ■; 

OI=DI+(VI { I)-DDI*U{ U IS ) )**2 
SDCn, IS)=VI( I) 

275 CONTINUE • . N 

0( IS)=DS0RT(D1 )/D01*5PR , 

IF( IPASS( IS).E0.2) GO TO 280 ‘ 

SUM=SUM+DI * 

TEST D AGAINST TO AND DELETE IF NECESSARY 
IF(D( ISI.LT.DMAX) GO TO 280 
DMAX=D(IS) 

ISMAX=IS 
280 CONTINUE 

IFIDMAX.LT.TD) RETURN 
IPASSI ISMAX)=2 

GO BACK AND MAKE ANOTHER PASS THROUGH THE DATA 
GO TO 130 
400 CONTINUE 
C DELETE WHOLE EVENT 
DO 410 IS=1,NSTE 
410 IPASSf IS )=2 
NSUSE0=0 

RETURN * I 

420 CONTINUE I 

C DELETE FOR INSUFFICIENT GOOD OBSERVATIONS | 

K0DE=2 ' 

GO TO- 400 

C DELETE FOR INSUFFICIENT GEOMETRICAL SEPARATION BETWEEN OBSERVATIONS 
430 CONTINUE 
K0DE=3 
GO TO 400 
END 
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INTEGER FUNCTION KSTAID(ID) 

■ eOHMON/STAORD/KOROER( 150) 

COMMON/NSTA/NST/i 
KSTAI0=0 

C SEARCH TABLE OF STATION IDENTIFIERS FOR THE INTERNAL NUMBER OF THIS STATION 
DO 10 I=1,NSTA 
IF(KORDER( I I.NE.ID) GO TO 10 
KSTAID=I 
RETURN 
10 CONTINUE 
RETURN 
END 


reproducibility of inE 
I PAGE IS POOR 
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SUBROUTINE VERSOL ( ORGMAT, VERMAT , M) 

IMPLICIT REAL=«=8<A-HtO-Z) 

COMMENT I IS THE NUMBER OF ROWS IN ORGMAT» AND M IS I PLUS THE NUMBER OF 
C OF UNKNOWN COLUMNS- THE ORIGINAL VALUES OF ORGMAT. ARE RETAINED- 

DIMENSION ORGMAT {I,M), VERMAT (I,M) 

DIMENSION P(72) 

N=I-1 
M1=M~1 
DO 1 J=1tI 
DO 1 K=1,I 

1 VERMAT«J,K)=0RGMAT(J,K) 

DO 5 K=1,I 

00 2 J=1,MI 

2 P (J) « VERMAT « 1»J+1)/VEPHAT (Xvl) 

P(M)=1.0DO/VERMAT( 1, 1 ) 

DO 4 L=1,N 
DO 3 J=1»MI 

3 VERMAT (L.J5 = VERMAT (L+1»J+1) - VERMAT (L+ltl) ♦ P(J) 

4 VERMAT {L»M) = - VERMAT (L+l»l) * PtM) 

DO 5 J=1»M 

5 VERMAT n»J) = P( J) 

RETURN 

END 


SUBROUTINE UVWD( A , B , PHI » L AMDArH, U, V, W ) 

DOUBLE PRECISION PHI ^ LAMOA, N, E2» FAC,U, V, W,SP 

REAL=S=a A,B,H 

E2=1.0-(B/AJ=i*’S=2 

SP=DSIM{PHI) 

N=A/DSORT( l-0-E2-SP*SP) 

FAC=(N+H)*DCOS( PHIS 
U=FAC*DCOS(LAMDA ) 

V=FAC^DSIN( LAMOA J 

W=(N*(1.0-E2J+H)*SP 

RETURN 

END 

SUBROUT I NE UVWTG2 \ UVW , DATUM, PHI , L AM, H ) 

C CONVERT RECTANGULAR TO GEODETIC COORDINATES 
C ALIAS FOR UVWTG 

IMPLICIT REAL*8{A“Z) 

DIMENSION UVWI3),DATUM!2) 

LAM=DATAN2(UVW( 2) ,UVW( 1 ) ) 

IF(LAM.LT-O.O) LAM=LAM+6. 283 185 307 1795800 
OME2=(DATUM(2)/OATUM(1) )**2 
E2=1.0“0ME2 

P=DSORT(UVW( 1)*=!^2+UVW{2)*-2) 

WP=UVW(33/P 
TP1^=WP/0ME2 
PHI1=DATAN{TP1 ) 

5 TTP=TP1*TP1 

SECP=DSORT{ 1.0+lTP) 

N=OATUM( I )--;=SECP/DSORT( 1 . 0 + OME2-TTP J 
H = P+-SECP-N 

TP2=WP/{ 1.0-E2*N/(N+H) ) 

PHI=DATAN{TP2 ) 

IF(DABS(PHI-PHI1 ).LT.1.D-12) RETURN 

PHI1=PHI 

TP1=TP2 
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GO TO 5 
END 

DOUBLE PRECISION FUNCTION ANRADDt ISGM, IDEG,MIN,SEC) 
INTEGER=^2 MINUS/IH-/, PLUS/1H + /, AMPSAN/lH6/t ISGN, IDEG,M!N 
DOUBLE PRECISION SEC 
IF( lOEG.GE.O) GO TO 10 
1SGN=MINUS 
IDEG=-IDEG 
10 CONTINUE 

ANRADD=(DFLOAT( ( 1 0EG=!=60+MI N )*60 ) +SEC > /206264.80625DO 

IF{ ISGN.E0.HINUS)ANRADD=-ANRADD 

IF(ISGN.EO.AMPSAN) 1SGN=PLUS 

RETURN 

END 


FUNCT I ONMJDI DATE, MONTH, YEAR) 

COMPUTATION OF MODIFIED JULIAN DAY 
INTEGER^2 DATE, YEAR 
DIMENSIONMOMTHS( 2, 12) 

DATAM0NTHS/3HJAN,0,2HFE8,21,3HMAR,59,3HAPR,90,3HNAY,120,3HJUM, 15i, 
13HJUL, 181, 3HAUG, 212, 3HSEP,2A3, 3H0CT,273, 3HNQV,30A.,3HDEC,33^/ 
ID=365=<=( YEAR~50) + (YEAR~49 )/4 
00201 = 1,12 

I F(MONTHeEO. MONTHS! 1,1 ) )GCT025 

20 CONTINUE 

1F(M0NTH.LE.12) GO TO 21 
WRITE! 6, 6001 ) MONTH 

6001 FORMAT! 3X,22HM0NTH NAME MISPELLEO , A3 ) 

STOP 

21 I=MONTH 
MONTH=MONTHS( 1, n 

25 CONTINUE 

ID=ID+M0NTHS!2, I ) 

IFCMOD! YEAR^^1,4).E0.0.AND.KGT,2) ID=ID+1 

MJD=ID+DATE+33281 

RETURN 

END 


SUBROUTINE DMSTR ! A , R , N , MS A , MSR ) 
IMPLICIT REAL-B! A-H,0-Z ) 
DIMENSION A! 1 ) ,RI 1 ) 

DO 20 1=1, N 
DO 20 J=1,N 
IF(MSR) 5,10,5 
5 IF! 1-J) 10, 10,20 
10 CALL LOG! I, J, IR,N,N,MSR) 

IFIIR) 20,20,15 
i5 R!IR)=0.0 

CALL LOC! I, J, IA,N,N,MSA) 

IF(IA) 20,20,18 
IR R! IR)=A! lA ) 

20 CONTINUE 
RETURN 
END 
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